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DIGITAL COMPUTER PROGRAMS FOR GENERATING' 

OBLIQUE ORTHOGRAPHIC PROJECTIONS 
AND CONTOUR PLOTS 

By Gary L. Giles 
Langley Research Center 

SUMMARY 

User and programer documentation is presented for two programs for automatic 
plotting of digital data. One of the programs generates oblique orthographic projections 
of three-dimensional numerical models and the other program generates contour plots of 
data distributed in an arbitrary planar region. The user documentation gives a general 
description of the computational algorithms, user instructions, and complete listings of 
the programs. Several plots are included to illustrate various program options, and a 
single example is described in detail to facilitate learning the use of the programs. 

INTRODUCTION 

Studies based on mathematical models of structural configurations using a digital 
computer yield large amounts of output data which must be comprehensively analyzed. 
Graphical presentation of this data is often needed for effective evaluation by an analyst. 
The present paper contains a description of computational algorithms as well as user and 
programer documentation for two computer programs which generate graphical displays 
of data. These programs generate (1) oblique orthographic projections of three- 
dimensional models and (2) contour plots of data distributed in an arbitrary planar region. 
The programs were developed initially for the graphical presentation of cumbersome out- 
put data resulting from structural analyses by finite- element methods. However, the 
programs have been extended for display of a general class of data. Several plots are 
included in this paper to illustrate various program options, and a single example is 
described in detail to facilitate learning the use of the programs. 

SYMBOLS 

AjA^jAjjAjj areas of triangles shown in figures 7 and 11 

Ag,A0,A^ transformation matrices for rotation about Y-, X-, and Z-axis, 
respectively (see eqs. (2) to (4)) 



a,b,c 

h 

Li> L jj L k 

h 


M 

N 


s 

w 

w c 

X,Y,Z 

Xo>Y 0 ,Z 0 


x,y,z 

x,y,z 


x 0 Y 0 > z o 


x p>y p ,Zp 


ft 


0 

0 

* 


length of sides of triangles shown in figure 11 
radius of circle inscribed in triangle 
triangular area coordinates defined in equation (8) 

distance between a particular grid point and centroid of ith surrounding 
element 

number of grid points contained by a finite element 

number of elements which connect at a particular grid point (eq. (7)) 

perimeter of a triangle 

control variable represented by contour plots 
constant value of a desired contour line 
coordinate system fixed in model 
coordinate system containing viewing planes 
coordinates of a point in X,Y,Z coordinate system 
coordinates of element centroid 

coordinates of a point in Xq/Y^Zq coordinate system 
coordinates of generic point p 

weighting function to interpolate data values from ith element centroid to a 
particular grid point 

Euler angle rotation about Y-axis, performed second 
Euler angle rotation about X-axis, performed third 
Euler angle rotation about Z-axis, performed first 
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GENERAL PROGRAM FEATURES 


Programs to generate graphical displays of large amounts of output data are an 
important aspect of any analysis system. Many of the programs which have been devel- 
oped for this purpose are reviewed in references 1 and 2. The two programs which are 
described in detail in this paper are distinguished by provision for generality and ease of 
use, as well as numerous display options, computational speed, and efficient use of core 
storage. These programs were developed initially for the graphical presentation of 
cumbersome output data resulting from structural analyses by finite -element methods. 
However, the programs have been extended for display of a general class of data. A 
general discussion of the programs is given in reference 3 and their incorporation into a 
computer-aided design system is discussed in reference 4. 

In structural analysis, finite-element or finite-difference methods are capable of 
providing displacements, stresses, and vibration or buckling mode shapes at many dis- 
crete points on a structural model. Manual reduction of the data is a time-consuming 
task which sometimes necessitates reducing only that data associated with particular 
regions of the overall model considered to be critical. Errors in the analyst's judgment 
can prevent the detection of other critical regions. Graphical display of analytical infor- 
mation by oblique orthographic projections and contour plots over large portions of the 
model, on the other hand, allow the analyst to assimilate accurately and evaluate the data 
with a consequent reduction in time and manpower cost. For computer-aided design 
programs as described in references 4 and 5, verification of data between steps in the 
design process is often needed. Visual display of oblique orthographic projections or 
contour plots by means of a cathode-ray tube (CRT) allows timely assessment of design 
information needed to make intelligent changes in the design variables during program 
execution. 

Oblique orthographic projections are useful in checking input data describing a 
numerical model. The topology of a finite- element model for structural analysis is 
described in a user-prepared input deck which includes a set of cards containing grid 
point identification numbers and corresponding spatial coordinates (x,y,z) and a set of 
cards containing grid point connections for rod, beam, and plate elements. The program 
described herein includes options for plots of numerical models annotated with grid point 
or element numbers. One option of the program allows boundaries of an isolated portion 
of the model to be specified by cutting planes, and detailed inspection can be made of that 
selected region. Also, exploded views can be generated which separate the elements in 
a model to provide clarity in detecting the absence or presence of model elements. 
Structural deformations calculated from a particular analysis can be superimposed on 
the grid point coordinates of the undeformed structure. These displacements can also be 
represented as vectors extending from the grid points. The program accommodates 
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lineal elements with two grid points and triangular or quadrilateral plate elements with 
three or four grid points. 

Another method of displaying output data is in the form of a contour plot. A com- 
puter program is described herein that generates contour plots of data relative to a 
planar mathematical reference surface which can have irregular boundaries or internal 
cutouts. The program input data require a set of grid points and connecting triangular 
and/or quadrilateral elements defining a planar region as well as the data to be repre- 
sented by contour lines which can be specified at either the grid points or centroids of 
the connecting elements. The contour lines are labeled to provide quantitative informa- 
tion in a graphical form. For the most general case, the grid points and finite elements 
of a model are located by three-dimensional coordinates. As developed for the plotting 
routine herein, the reference surface is merely a planar projection of the three- 
dimensional surface. This arrangement is adequate for nearly flat surfaces; however, 
the presentation of contour plots on the two-dimensional projection of a three-dimensional 
surface having a large curvature may be difficult to interpret. In this case, a two- 
dimensional planar surface must be developed from the three-dimensional surface by the 
user outside the program. 

The input decks for the two programs are similar and thus allow for the programs 
to be used as a set. Both programs contain options for selecting various plotting equip- 
ment including CalComp, Varian, and CRT displays. Finally, this paper includes com- 
plete listings of the programs which contain comment cards to aid the user. 

OBLIQUE ORTHOGRAPHIC PROJECTIONS 
Specification of View 

A three-dimensional analytical model consists of a user-prepared set of grid points 
with given spatial coordinates (x,y,z) and a set of elements (e.g., rod, beam, triangular or 
quadrilateral elements) connected at the grid points. An example of an oblique ortho- 
graphic projection of a finite-element structural model of an airplane is shown in fig- 
ure 1. This model was used in preliminary design studies focused on the wing structure. 
The fuselage is represented by a structural model having a simplified rectangular geo- 
metric cross section to reduce model complexity, but the model does provide proper 
overall stiffness and mass distributions to represent symmetric behavior of the entire 
airplane. Oblique orthotropic projections allow a model to be viewed in any selected 
orientation. Euler angle transformations are used to specify orientation of the model to 
be projected. This transformation resolves the coordinate system of the model to a 
principal plane (i.e., viewing plane) on which the display is to be plotted. The model 
coordinate system is coincident with the coordinate system containing the viewing plane 
when all the rotation angles are zero. 
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Various approaches can be taken in formulating the Euler angle transformations. 
Herein, the viewing planes are fixed in space and the model rotated about its body axes. 
In this approach the principal plane is the plane normal to the user's line of vision. The 
model can then be rotated in space until the Euler angles giving the desired view are 
determined. The rotations (cp,6,p) of the body about the model axes (X,Y,Z) for this 
approach are shown in figure 2. The order of the Euler angle rotations is taken as p, 
6, and p in the program and transform coordinates of a point on the model (x,y,z) to 
viewing plane coordinates (x 0 ,y 0 ,z 0 ) as 







where 



COS p 

sin p 
0 


-sin p 0 

cos p 0 

0 1 


( 1 ) 


( 2 ) 



cos 9 
0 

-sin 6 


0 sin 6 

1 0 

0 cos 6 


(3) 



1 0 0 

0 cos p -sin cf > 

0 sin p cos p 


(4) 


These transformations are equivalent to those used in reference 6. The option to specify 
which of the principal planes is to be used as the viewing plane determines which two rows 
of equation (1) are used. The model can be projected on either the X q -Y o , X q -Z o , or 
Y q -Z o planes and is shown projected on the X Q -Z o plane in figure 2. 
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In addition to plots of the undeformed structure, oblique orthographic projections 
can be used to display the deformed structure by adding given displacements to the coor- 
dinates of the grid points before transformation and subsequent plotting of elements. An 
option is also available to plot the undeformed model and represent displacements by 
vectors extending from the grid points in the undeformed structure to the location of the 
grid points in the deformed structure. 

Exploded Plots 

Oblique orthographic projections may be used to check the topology of an analytical 
model. Grid points with erroneous coordinates are easily detected. Often, however, 
absence or presence of elements in the model cannot be determined from a conventional 
oblique orthographic projection. For example, a line element coincident with an edge of 
a triangular or quadrilateral element could not be detected. To show clearly each ele- 
ment, an algorithm for generating exploded oblique orthographic projections was devel- 
oped. Such a plot is shown in figure 3 for a selected portion of the model shown in 
figure 1. 

The algorithm for generating exploded oblique orthographic projections follows: 

For each element in the model - 

(a) Calculate the centroid of the element by using 


M 

M 

M 


' t y i 

1 Z i 

i=l 

y = — - — , and 

z = i=1 

M ’ 

J M 

M 


where M is the number of grid points contained by a finite element 

(b) Multiply coordinates of the element by a reduction factor to give the location of 

the reduced element centroid and vertices 

(c) Translate the reduced element until its centroid matches the location of the 

centroid of the element before reduction. 

Although this procedure provides exploded plots, the resulting sides of reduced tri- 
angular elements generally are not equidistant from the sides of the original element. 

To correct this behavior, the incenter (center of a circle inscribed in a triangle) is used 
for the reference point of the triangles during translation. The equations for the incenter 
of a general triangle are given in appendix A. 
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Sectioning or Cutting Planes 


It is often desirable when checking input for a model to be able to isolate a portion 
of the model for detailed examination. The program for generating oblique orthographic 
projections described herein contains the option to specify sectioning or cutting planes to 
isolate such a portion. For example, in the case of a finite-element model of an aircraft 
structure, a particular fuselage frame could be generated by placing cutting planes imme- 
diately in front of and behind the frame. Annotation of grid point and element numbers 
show clearly on such a plot but might be obscured on a plot of the entire model. This 
sectioning option was used to select the structural members of a fuselage frame-wing 
spar cross section shown in figure 4 from the complete model shown in figure 1. 

CONTOUR PLOTS 
General Description 

A contour plot is a graphic representation of a scalar variable as function of two 
independent variables, w = f(x,y). Here w is a quantity such as deformation, pressure, 
or stress which is distributed over some specified area of the X-Y plane. The set of 
values of w are herein referred to as control variables and have been measured or cal- 
culated at given locations called control points on a reference or datum surface. The 
locus of lines of constant values of w (contour lines) over the reference surface can be 
determined from the given values of w at the discrete control points. This section 
describes a method for generating contour plots which is particularly adaptable for dis- 
playing the results of a structural analysis by finite-element techniques. The data to be 
represented by contour plots can be referenced to either the grid points or to the element 
centroids. These alternatives are particularly convenient since finite -element analyses 
generally evaluate displacements at grid points and stresses at the element centroids. A 
contour plot of a selected vibration mode shape for the model presented in figure 1 is 
shown in figure 5(a). Oblique orthotropic projections of the same mode shape are also 
included in figures 5(b) and 5(c). 

Definition of Contour Surface 

The data to be represented by contour lines can be visualized as a set of points 
located perpendicular to and at a distance equal to the magnitude of the data value from 
each control point. A surface which passes through this set of data points will be 
referred to as a contour surface as shown in figures 5(b) and 5(c). Contour lines (lines 
connecting constant values of w) are projections onto the datum surface of intersections 
of regularly spaced planes, parallel to the datum surface, with the contour surface. The 
regular spacing of the planes, hence contour interval, occurs in units of 10 in figure 5(a). 
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The primary difference between various contour plotting programs is the method 
of representing the contour surface. Options for this representation include least- 
squares fitting of the data with polynomial surfaces and surface splines (ref. 7). Such 
representations result in smooth contour lines, however, considerable computational 
time is often required for the surface fitting, and difficulties are sometimes encountered 
with irregular boundaries or cutouts. Herein, the contour surface is idealized as a 
series of triangular planes with their vertices connected at points where the control vari- 
ables are located. The control variables may be located at either the grid points or ele- 
ment centroids of the model. A computationally efficient algorithm which uses linear 
interpolation is used to determine the location of contour lines on each of the triangular 
planes. 

Prior to generating contour lines, a triangular mesh must be generated over the 
datum surface and the magnitude of the control variable at each of the vertices in the 
mesh determined. The triangular and/or quadrilateral elements in the model are used 
to triangulate the datum surface. When control variables are given at the grid points, the 
quadrilateral elements are divided into two triangles. Such a triangular mesh is shown 
in figure 6(a). 

Control variables located at the centroid of the elements, on the other hand, require 
a different method of triangulating the datum surface since the control points (element 
centroids) do not coincide with the grid points. For this case, a triangular mesh with 
control points at both the centroid of the two-dimensional elements and the grid points of 
the structure is generated. The centroid of each element can be readily calculated from 
the data describing its corner points. Since the control variables are given only at the 
centroidal control points, a suitable method is required to interpolate this data to each of 
the grid control points. The following steps describe the method of interpolation used in 
this algorithm and is illustrated in figure 6(b): 

(a) Determine the elements that are connected to a particular grid point G (five ele- 

ments are connected to the grid point G in fig. 6(b)) 

(b) For each element connected to the grid point, calculate the centroid which is the 

location of the given control values 

(c) Calculate a weighted value of the contour value at the grid point Wq based on 

the contour values at the surrounding element centroids as follows: 

N 

W G = 2 W i^i (6) 

i=l 

where w^ is the given data value at the centroid of the ith element, ^ is a 

weighting function for that value, and N is the total number of elements con- 
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nected to grid point G (N = 5 in fig. 6(b)). The weighting function selected 
for use herein is 

1_ 

*- 4 - < 7 > 



where ^ is the distance between the grid point G and the centroid of the 
ith element surrounding the grid point. (See fig. 6(b).) 


Constructing Contour Lines 

The construction of contour lines which are projections onto the datum surface of 
the intersections of planes parallel to the datum surface and the contour surface is 
described in this section. For convenience of discussion, the datum plane will be 
assumed to lie in the X-Y plane although it could lie in another plane defined by two 
independent variables. The contour lines are generated by considering one triangular 
element of the contour surface mesh at a time. Area coordinates of a triangular plane, 
as described in reference 8, are used to give simple expressions for the location of a 
point on such a plane. 

The area coordinates Lj, Lj, and of a point p on the triangle shown in fig- 
ure 7 are given by 



A 



(8) 



j 


and 


^ + Lj + L k - 1 


(9) 


The coordinates of a generic point p on the triangular plane (one facet of the contour 
surface) with vertices (x^y^w^, ( x j»y and ( x k ,y k’ w k) can be ex P ressed as 
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x p = ^ + Lj Xj + L k x k 

(10) 

y p = L^ + L jYj + L k y k 

(11) 

w p = L^ + L jWj + L k w k 

(12) 


where x and y are measured in the reference surface and w is the magnitude of the 
control value measured perpendicular to the reference surface. 


For a contour line, w is equal to w c where w c is the constant value of a 
desired contour line. Procedures contained in the computer program require that the 
data be scaled so that the contour lines have positive or negative integer values. The 
number of contour lines, if any, and their associated magnitudes that are contained in a 
particular triangular plane are determined from the minimum and maximum values at 
the vertices and the prescribed contour interval. For those triangles which contain con- 
tour lines, the intersection of a particular contour line w c with an edge of the triangular 
plane can be determined through the use of equation (12) and the definition of area coordi- 
nates given in equations (8). For example, along the edge i-j, L k is equal to zero and 
Lj equals 1-1^ from equation (9). Solving equation (12) for along the edge i-j 
yields 


Li = 


w c- w j 


W; 


W, 


(13) 


Similarly along the edge j-k, 


Lj — 0 — 1 - Lj 


Thus, 


U = 


w - w, 
c k 


w 


j ' w k 


Along the edge k-i, 


L j=° 


Li - 1 - Lk 


which yields 


L k = 


w c - Wi 


w k 


w. 


(14) 


(15) 
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Two of the intersections given by equations (13) to (15) of the contour line and tri- 
angle edges will be between the vertices of the triangle and the third will be outside the 
triangle. The point outside the triangle is indicated by values of area coordinates which 
are less than zero or greater than unity and this point is disregarded. The x-y coordi- 
nates at the intersection of the contour line and an edge are solved from equations (10) 
and (11); for example, the Cartesian coordinates at an intersection along the edge i-j 
are given as 


x = X t + (Xj - x i )Lj 

(16) 

y = Yi + <yj - y^Lj 

(17) 


The end points (x-y coordinates at intersections) of each contour line within the triangular 
facet are then connected by straight lines. This procedure is repeated for each triangular 
facet over the contour surface. Segments of the contour lines on adjacent elements meet 
at the common element sides and result in all contour lines being connected properly 
over the entire surface. 

The plot is completed by drawing the boundary of the surface and labeling the mag- 
nitude of the contour lines. Surfaces with irregular boundaries or cutouts are allowed. 

The boundary lines are determined by summing the corner angles of each triangular facet 
connected to a common grid point. If the total angle is less than 360° (point not completely 
surrounded by elements), the point is obviously on a boundary. All elements in the datum 
plane are then tested and adjacent vertices which are boundary points are connected to 
form the complete boundary. Each contour line which intersects the boundary is labeled 
with its corresponding magnitude. An option is also available to label the contour plot at 
specified locations within the boundary of the contour plot. This option is necessary to 
allow labeling of contour lines which do not intersect the boundary of the datum surface. 

DESCRIPTION OF COMPUTER PROGRAMS 

Computer Hardware and System Requirements 

The plotting programs described herein are written in FORTRAN IV for the CDC 
6000 series computers. The programs are operational on the CDC SCOPE 3.2 system 
(NASA Langley Research Center version) and make use of subroutines from its graphic 
output system. This system includes CalComp, Varian, and CRT capabilities. These 
programs have been utilized on all the above equipment. These programs may not be 
operational directly on other computer systems without some modification. Therefore, 
a description of the graphics subroutines that are used is given on comment cards in the 
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subroutine called DOCMNT. This information is adequate for conversion to another 
graphics system. 


Storage Allocation 

Dynamic storage allocation is used in both programs for efficient accommodation of 
models of varying complexity. All large arrays used by the programs are stacked in 
blank COMMON designated in the MAIN program as ZZZ. The starting location of each 
array is calculated in the MAIN programs and the total blank COMMON requirement is 
printed out. The amount of blank COMMON is determined in general by the number of 
grid points, number of elements, and number of displacements considered at the grid 
points. The total required field length is the address of the first location in blank 
COMMON plus the blank COMMON length requirement. 

Since the field length is problem dependent, the user must estimate the required 
field length for any new problem. An approximate formula (in decimal) for the field 
length required for the oblique-orthographic-projection program is 

FL 10 = 18 000 + (4 + NUDISP + NVDISP + NWDISP)NNDEST 

where NNDEST is equal to the number of grid points and NUDISP, NVDISP, and NWDISP 
are the indicators for displacements in the x-, y-, and z-direction, respectively (must be 0 
or 1). The field length for this program is independent of the number of elements which 
are stored on an auxiliary storage unit. A corresponding formula for the contour-plot 
program is 


FL 10 - 17 000 + 5(NNDEST) + 6(NE LEST) 

where NNDEST is equal to the number of grid points and NELEST is the number of ele- 
ments. For operation of either of the programs on the CRT display at Langley Research 
Center an additional storage allocation of 3500 decimal words is required. 

Programs and Subroutines 

The two programs for oblique orthographic projections and contour plots are sep- 
arate but intended to be used as companion programs. The input decks are designed to 
be similar, and many of the variables controlling options in the programs have identical 
names and purposes. 

Complete listings of the programs are given in appendixes B and C as well as 
information on flow logic for the programs. Comment cards are incorporated in the 
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computer programs to give the purpose of each subroutine and to indicate the operations 
that are performed. 


COMPUTER PROGRAM USE 
General Setup of Input Deck 

In general the input data decks for both the oblique-orthographic-projection pro- 
gram and the contour -plotting program consist of six separate groups of data as shown 
schematically in figure 8. These groups are as follows: 

(1) A single card containing any desired title information 

(2) NAMELIST OPTION containing values to allocate storage in blank 

COMMON and control values specifying various program options 

(3) a geometry deck containing the specification of grid points and con- 

necting elements in the model 

(4) an optional single card used to identify the deck of data to be 

plotted 

(5) the deck of data to be plotted 

(6) NAMELIST PICT containing values to specify the type of plot 

desired and what information is to be included on the plots. 

Selected groups of this basic input deck can be repeated to make different plots of the 
same data or to input additional groups of data to be plotted during a single run. Details 
of the data contained in the input decks for both programs are described in the two sub- 
sequent sections of this report. 

Options are provided to input geometry data and data to be plotted in several dif- 
ferent forms by selection of appropriate input subroutines. Generality is provided in 
subroutines GEOM1 and DATA1 by user specification of the format by which the data is 
to be read directly in the input deck. A specialized subroutine GEOM2 provides for 
direct input of bulk data decks for the NASTRAN (ref. 9) program. The plotting programs 
herein handle only one- or two-dimensional elements (no solid elements) containing up to 
four grid points. In addition, only a single coordinate system can be used and therefore 
models defined using alternate coordinate systems must be transformed to the basic 
coordinate system of the model before plotting. Grid point numbers and element numbers 
need not be sequential since they are reordered internally in the programs. Data to be 
plotted can be read from magnetic tape by using subroutine DATA5. In addition, dummy 
subroutines GEOM9 and DATA9 are included for user-prepared routines which read 
geometry (grid points and connecting elements) information and data to be plotted. 
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Input for Oblique Orthographic Projections 

The input data deck is illustrated in the sequence shown schematically in figure 8. 
The data must be constructed in the order shown in figure 8 and is described in detail 
in this section. For clarity, no zeros appear in variable names described in this section. 

Title card .- This single card contains any desired alphanumeric information in 
columns 1 to 80. The title will appear at the beginning of the plots. 

NAMELIST OPTION .- This NAMELIST contains values to allocate storage in blank 
COMMON and control values specifying various program options. 


FORTRAN 

name 

Default 

value 


Description 

NNDEST 

200 

Estimated number of grid points, must be equal to or 
greater than the actual number of grid points 

NUDISP 

0 

0 

no displacement data in x-direction 



1 

data including displacements in x-direction 

NVDISP 

0 

0 

no displacement data in y-direction 



1 

data including displacements in y-direction 

NWDISP 

0 

0 

no displacement data in z -direction 



1 

data including displacements in z-direction 

KGEOM 

1 

Specifies the subroutine and corresponding method of 


input for model geometry 

1 subroutine GEOM1 to read in grid points and elements 
from cards with user- specified/format 

2 subroutine GEOM2 to read in NASTRAN bulk data deck 
with data in column widths of 8 

9 subroutine GEOM9, a user-supplied subroutine 

KDATA 1 Specifies the subroutine and corresponding method of 

input for displacement data 

1 subroutine DATA1 to read in displacement data from 
cards with user-specified format 

5 subroutine DATA5 to read in displacement data from 
TAPE 20 

9 subroutine DATA9, a user-supplied subroutine 

NVALUS 0 Used if KDATA = 5; specifies the number of grid points 

at which displacement data is to be read from TAPE20 
(must be less than or equal to NNDEST) 

IRESEQ 1 Grid point numbers are stored in the program from 1 to 

the total number of grid points 


14 



FORTRAN 

name 

Default 

value 

Description 

IRESEQ 

1 

0 no internal resequencing of grid points necessary; 
they are already ordered in ascending order starting 
with 1 

1 resequence grid points in the same order as they are 
input 

KPLOT 

1 

Specifies the type of output device to be used 

1 CalComp 

2 CalComp with plotting speed reduced for Leroy pens 

3 VARIAN 

4 cathode -ray -tube (CRT) console (set up for CDC 250 
scopes at Langley Research Center) 

XSPACE 

10.0 

Space between plots in x-direction, in inches 

PSIZE 

25.0 

Paper size in y-direction, in inches (used in scaling of 
plots to insure this dimension is not exceeded) 

IDCASE 

0 

0 no identification card preceding decks of displacement 
values 

1 identification card preceding decks of displacement 
values 


Geometry input data deck .- This portion of the input deck contains the geometry 
definition of grid points and connecting elements. The deck has one of the following 
forms, depending on the value of KGEOM specified in NAMELIST OPTION. 

KGEOM = 1 

(a) A single card containing the word FORMAT in columns 1 to 6 and a variable 
format corresponding to the format of the grid point cards with the left paren- 
thesis starting in column 11 and up to column 80 may.be used. 

(b) Deck of grid point cards. Each card contains 4 values: grid point number 
(integer), x-coordinate (real), y-coordinate (real), and z-coordinate (real). 

The format is specified in (a) above. 

(c) A single card containing the word ENDGRID in columns 1 to 7. 

(d) A single card containing the word FORMAT in columns 1 to 6 and a variable for- 
mat corresponding to the format of the element cards with left parenthesis start- 
ing in column 11 and up to column 80 may be used. 
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(e) Deck of element cards. Each card contains 5 integer fields which are the ele- 
ment number and grid point numbers at the vertices of the elements. For tri- 
angular elements the last integer field must be blank or zero. For rod or beam 
elements the last two integer fields must be blank or zero. The format is 
specified in (d) on the preceding page. 

(f) A single card containing the word ENDGEOM in columns 1 to 7. 

KGEOM = 2 

(a) A single card containing the word LINEEL in columns 1 to 6 and up to 9 
NASTRAN lineal element connection names, which are adjusted to the left in 
field widths of 8, starting in column 9 (cols. 9 to 16, 17 to 24, . . ., 73 to 80). 
This card can be omitted if lineal elements are not used for the plot. 

(b) A single card containing the word TRIAEL in columns 1 to 6 and up to 9 
NASTRAN triangular element connection names, which are adjusted to the left 
in field widths of 8, starting in column 9 (cols. 9 to 16, 17 to 24, . . ., 73 to 80). 
This card can be omitted if triangular elements are not used for the plot. 

(c) A single card containing the word QUADEL in columns 1 to 6 and up to 9 
NASTRAN quadrilateral connection names, which are adjusted to the left in field 
widths of 8, starting in column 9 (cols. 9 to 16, 17 to 24, . . ., 73 to 80). This 
card can be omitted if quadrilateral elements are not used for the plot. 

(d) A NASTRAN bulk data deck. Only the GRID cards and the element connection 
cards with names matching those given on the LINEEL, TRIAEL, and QUADEL 
cards will be used for the plot. All other cards in the NASTRAN bulk data deck 
will be ignored. 

(e) A single card containing the word ENDGEOM in columns 1 to 7. 

KGEOM = 9 

Calls subroutine GEOM9 which is prepared by the user to read geometry data. 

Case identification card .- This single card is omitted if IDCASE = 0 is specified 
in NAMELIST OPTION. If present, this card contains any desired alphanumeric informa- 
tion in columns 1 to 80. The identification will appear on all plots of the case. 

Deck of data to be plotted . - This deck contains displacement sets at grid points for 
the oblique-orthographic-projection program. A displacement set for each grid point is 
defined to contain from 2 to 4 values (i.e., a grid point number and displacements corre- 
sponding to NUDISP, NVDISP, or NWDISP equal to 1). The deck has one of the following 
forms, depending on the value of KDATA specified in NAMELIST OPTION. 
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KDATA = 1 


(a) A single card containing the word FORMAT in columns 1 to 6 and a variable 
format for the data cards with the left parenthesis starting in column 11 and up 
to column 80 may be used. If displacements are included for more than one 
grid point per card, the number of grid points per card must be entered as an 
integer in column 8. 

(b) Deck of displacement sets. There can be multiple displacement sets per card 
or the set can extend to more than one card (often the case with NASTRAN 
punched output) which can be handled with a format for reading multiple cards. 

(c) Blank card or cards to end the data deck. The number of blank cards must cor- 
respond to the number of cards read at one time by the specified variable format. 

KDATA = 5 

Reads NVALUS (from NAMELIST OPTION) displacement sets from TAPE20. Each 

displacement set must have been written on TAPE 20 as an unformatted record. 

KDATA = 9 

Calls subroutine DATA9 which is prepared by the user to read displacement data. 

NAMELIST PICT .- This NAMELIST contains values to specify the type of plot 
desired and the information that is to be included on the plots. 


FORTRAN 

name 

Default 

value 

Description 

KHORZ 

1 

Integer designating the horizontal axis of the 
viewing plane where 1 = X Q , 2 = Y 0 , and 
3 = Z 0 

KVERT 

2 

Integer designating the vertical axis of the view- 
ing plane where 1 = X Q , 2 = Y 0 , and 3 = Z Q 

PHI 

0.0 

Angular rotation of model about its X-axis in 
degrees (must be performed third) 

THETA 

0.0 

Angular rotation of model about its Y-axis in 
degrees (must be performed second) 

PSI 

0.0 

Angular rotation of model about its Z-axis in 
degrees (must be performed first) 

NEWFR 

1 

1 frame change before plotting (a frame change 
resets the x-origin past previous plot by XSPACE 
given in NAMELIST OPTION and resets the 
y-origin at 0.0) 
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FORTRAN 

name 

Default 

value 

Description 

NEWFR 

1 

0 no frame change before plotting 

ISC ALE 

1 

1 automatic computation of proper origin 
location and scaling of plot 

2 user-specified origin and scaling 

PLOTSZ 

10.0 

Maximum dimension desired on completed 
plot, in inches (used for scaling if 
ISCALE = 1) 

XORGN 

0.0 

x-location of plot origin (used if ISCALE = 2) 

YORGN 

0.0 

y-location of plot origin (used if ISCALE = 2) 

PSCALE 

1.0 

Model size reduction factor (i.e., PSCALE is 
equal to actual model size divided by desired 
plot size (used if ISCALE = 2)) 

NOTAT 

0 

0 no numbering on plots 

1 numbering of grid points 

2 numbering of elements 

XLHT 

0.15 

Height of integers specified by NOTAT, in 
inches (must be 10.07) 

KDISP 

0 

0 plot of undeformed structure 

1 plot of deformed structure 

2 exploded plot 

3 displacements represented by vectors 

IDMAG 

2 

1 direct magnification of displacement data 
by DMAGS 

2 scaling of displacement data to a maxi- 
mum value of DMAGS 

DMAGS 

1.0 

Magnification of displacements (if 
KDISP = 1 or 3) 

Reduction factor of elements (if KDISP = 2) 

KSYMXY 

0 

1 symmetry about X-Y plane 

KSYMXZ 

0 

1 symmetry about X-Z plane 

KSYMYZ 

0 

1 symmetry about Y-Z plane 


Symmetries are performed consecutively (i.e., a plate quadrant with KSYMXZ and 
KSYMYZ equal to one would yield a complete plate). 
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FORTRAN 

name 


Default 

value 


Description 


XXMAX, YYMAX, ZZMAX 1.0 E+20 
XX MIN, YYMIN, ZZMIN -1.0 E+20 

NDMAX 9999999999 


Locate cutting planes parallel to principal 
planes (X-Y, X-Z, Y-Z) to limit plot 

Maximum grid point identification number to 
be included in plot 


NDMIN 


0 Minimum grid point identification number to 

be included in plot 


NELMAX 

NELMIN 

KODE 


9999999999 Maximum element identification number to 
be included in plot 


0 Minimum element identification number to be 

included in plot 

0 Specifies control option after plot is complete 

0 last plot, exit from program 

1 read another NAMELIST PICT 

2 read a new set of displacement data, 
including a case identification card if present 

3 read a complete new set of input data, 
including a title card. 


This section describes a complete basic set of input data if KODE = 0 in NAME- 
LIST PICT. For KODE = 1, 2, or 3, additional sections of the basic deck must be 
repeated. The deck must end with NAMELIST PICT having a value of KODE = 0. 


Input for Contour Plots 

The input data deck is illustrated in the sequence shown schematically in figure 8. 
The data must be constructed in the order shown in figure 8 and is described in detail in 
this section. For clarity no zeros appear in the variable names described in this section. 

Title card .- This single card contains any desired alphanumeric information in 
columns 1 to 80. The title will appear at the beginning of the plots. 


NAMELIST OPTION .- This NAMELIST contains values to allocate storage in blank 
COMMON and control values specifying various program options. 


FORTRAN 

name 


Default 

value 


Description 


NNDEST 200 Estimated number of grid points, must be 

equal to or greater than the actual number of 
grid points 
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FORTRAN 

name 


NELEST 


KGEOM 


KDATA 


NVALUS 


IRESEQ 


Default 

value 


200 


1 


1 


0 


1 


Description 


Estimated number of elements, must be 
equal to or greater than the actual number of 
elements 

Specifies the subroutine and corresponding 
method of input for model geometry 

1 subroutine GEOM1 to read in grid points 
and elements from cards with user-specified 
format 

2 subroutine GEOM2 to read in NASTRAN 
bulk data deck with data in column widths of 8 
9 subroutine GEOM9, a user -supplied 
subroutine 

Specifies the subroutine and corresponding 
method of input for control variable data to 
be represented by contour lines 

1 subroutine DATA1 to read in data to be 
plotted from cards with user-specified 
format 

5 subroutine DATA5 to read in data to be 
plotted from TAPE20 
9 subroutine DATA9, a user-supplied 
subroutine 

Used if KDATA = 5; specifies the number of 
given control variable points at which data to 
be plotted is to be read from TAPE20 (can be 
either the number of grid points or the num- 
ber of elements) 

Grid point numbers are stored in the program 
from 1 to the total number of grid points 

0 no internal resequencing of grid points 
necessary; they are already ordered in 
ascending order starting with 1 

1 to resequence grid points in the same 
order as they are input 
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FORTRAN 

name 

Default 

value 

Description 

KPLOT 

1 

Specifies the type of output device to be used 

1 CalComp 

2 CalComp with plotting speed reduced for 
Leroy pens 

3 V ARLAN 

4 cathode -ray -tube (CRT) console (set up 
for CDC 250 scopes at Langley Research 
Center) 

INFOR 

1 

1 data to be plotted is specified at grid 
points 

2 data to be plotted is specified at element 
centroids 

XSPACE 

10.0 

Space between plots in x-direction, in inches 

KSIGN 

1 

-1 change signs of y-coordinates 
1 do not change signs of y-coordinates 

IDCASE 

0 

0 no identification card preceding decks of 
data to be plotted 

1 identification card preceding decks of 
data to be plotted 


Geometry input data deck .- This portion of the input deck contains specifications of 
grid points and elements describing the reference surface. The deck has one of the fol- 
lowing forms, depending on the value of KGEOM specified in NAMELIST OPTION. 

KGEOM = 1 

(a) A single card containing the word FORMAT in columns 1 to 6 and a variable 
format corresponding to the format of the grid point cards with the left paren- 
thesis starting in column 11 and up to column 80 may be used. 

(b) Deck of grid cards. Each card contains 3 values: grid point number (integer), 
x-coordinate (real), and y-coordinate (real). The format is specified by (a) 
above. 

(c) A single card containing the word ENDGRID in columns 1 to 7. 

(d) A single card containing the word FORMAT in columns 1 to 6 and a variable for- 
mat corresponding to the format of the element cards with left parenthesis start- 
ing in column 11 and up to column 80 may be used. 
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(e) Deck of element cards. Each card contains 5 integer fields which are the ele- 
ment number and grid points at the vertices of the element. For triangular 
elements the last integer field must be blank or zero. The format is specified 
in (d) from the preceding page. 

(f) A single card containing the word ENDGEOM in columns 1 to 7. 

KGEOM = 2 

(a) A single card containing the word TRIAEL in columns 1 to 6 and up to 9 
NASTRAN triangular element connection names, which are adjusted to the left 
in field widths of 8, starting in column 9 (cols. 9 to 16, 17 to 24, . . ., 73 to 80). 
This card can be omitted if triangular elements are not used for the plot. 

(b) A single card containing the word QUADEL in columns 1 to 6 and up to 9 
NASTRAN quadrilateral connection names, which are adjusted to the left in 
field widths of 8, starting in column 9 (cols. 9 to 16, 17 to 24, . . ., 73 to 80). 
This card can be omitted if quadrilateral elements are not used for the plot. 

(c) A NASTRAN bulk data deck. Only the GRID cards and the element connection 
cards with names matching those given on the TRIAEL and QUADEL cards will 
be used for the plot. All other cards in the NASTRAN bulk data deck will be 
ignored. 

(d) A single card containing the word ENDGEOM in columns 1 to 7. 

KGEOM = 9 

Calls a subroutine GEOM9 which is prepared by the user to read geometry data. 

Case identification card .- This single card is omitted if IDCASE = 0 is specified 
in NAMELIST OPTION. If present, this card contains any desired alphanumeric informa- 
tion in columns 1 to 80. The identification will appear on all plots of the case. 

Deck of data to be plotted .- This deck contains sets of a control variable value cor- 
responding to each grid point or element (depending on INFOR in NAMELIST OPTION) to 
be represented by contour lines. The deck has one of the following forms, depending on 
the value of KDATA specified in NAMELIST OPTION. 

KDATA = 1 

(a) A single card containing the word FORMAT in columns 1 to 6 and the variable 
format for the data cards with the left parenthesis starting in column 11 and up 
to column 80 may be used. If control variables are included for more than one 
grid point or element per card, the number of grid points or elements per card 
must be entered as an integer in column 8. 
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(b) Deck of data to be plotted. There can be multiple data value sets per card or 
the set can extend to more than one card (often the case with NASTRAN punched 
output) which can be handled with a format for reading multiple cards. 

(c) Blank card or cards to end the data deck. The number of blank cards must cor- 
respond to the number of cards read at one time by the specified variable format. 

KDATA - 5 

Reads NVALUS (from NAMELIST OPTION) sets of control variable values from 
TAPE20. Each set of control variables must have been written on TAPE20 as an 
unformatted record. 

KDATA = 9 

Calls subroutine DATA9 which is prepared by the user to read data to be plotted. 

NAMELIST PICT .- This NAMELIST contains values to specify the type of plot 
desired and what information is to be included on the plots. 

Description 

Specifies the type of plot to be generated 

1 layout of elements in reference surface 
without grid point or element labels 

2 element layout with grid point labels 

3 element layout with element labels 

4 contour plots without symbols at grid 
points 

5 contour plots with symbols at grid points 

x-location of origin of first plot, in inches 

y-location of origin of first plot, in inches 

Model size reduction factor (i . e . , PSCALE is 
equal to actual model size divided by desired 
plot size) 

Method of scaling control variable data to be 
plotted; contour lines have only integer val- 
ues annotated on the plot and, thus, the data 
must be scaled such that these integers will 
contain the desired number of significant 
digits; the definitions of WMAGS and ICNTRS 
depend on the value of ISCALE 
1 user specification of scale factors 


FORTRAN 

name 

NPLOT 


Default 

value 


XORGN 

YORGN 

PSCALE 


0.0 

0.0 

1.0 


ISCALE 
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FORTRAN 

name 

ISCALE 


WMAGS 


ICNTRS 


XLHT 
NX LAB 


XLAB 


Default 

value 


3 


100.0 


10 


0.15 

0 


all zeros 


Description 

2 program calculation of scale factors to 
give the user-specified number of significant 
digits in annotation of the maximum absolute 
contour line 

3 program calculation of scale factors to 
give WMAGS as the maximum value of data 

If ISCALE = 1, magnification of control vari- 
ables for annotation of contour lines on plot 
If ISCALE = 2, number of significant digits 
in annotation of maximum absolute contour 
line (WMAGS = 1.0, 2.0, 3.0, etc.) 

If ISCALE = 3, maximum value of scaled 
data, WMAGS =2.0 (maximum contour 
line is integer truncation of WMAGS) 

If ISCALE = 1, user-specified contour 
interval (difference in integer values of 
adjacent contour lines) 

If ISCALE = 2 or 3, approximate number of 
different contour line values; the contour 
interval is calculated by the program 

Height of integers to be annotated on plots, 
in inches (must be =0.07) 

Total number of lines parallel to Y-axis 
along which contour lines are labeled, must 
be =10 (all contour lines are labeled where 
they intersect with these selected lines); 
these labels are in addition to those auto- 
matically provided at boundaries of the con- 
tour surface 

Array of distances in x-direction from the 
origin to lines parallel to the Y-axis along 
which contour lines are labeled; there must 
be NXLAB of these values and they must be 
in units of the original (unsealed) model 
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FORTRAN 

name 

Default 

value 

Description 

NY LAB 

0 

Same as NXLAB for label locations parallel 
to X-axis 

YLAB 

all zeros 

Same as XLAB for label locations parallel 
to X-axis 

KODE 

0 

Specifies control option after plot is complete 


0 last plot, exit from program 

1 read another NAMELIST PICT 

2 read a new set of control variable values 
to be plotted including a case identification 
card if present 

3 read a complete new set of input data, 
including a title card 

The above NAMELIST comprises a complete basic set of input data if KODE = 0 
in NAMELIST PICT. For KODE = 1, 2, or 3, additional sections of the basic deck must 
be repeated. The deck must end with NAMELIST PICT having a value of KODE = 0. 

EXAMPLES OF DATA PREPARATION 

An example problem is presented to illustrate preparation of input data for both the 
oblique-orthographic-projection program and contour -plotting program as well as printed 
and plotted output. The problem presented herein is designed to exercise and illustrate 
use of various options of the programs. The example selected is the graphical display of 
the normal displacements of a square flat plate under uniform lateral pressure. A finite - 
element structural model of the flat plate with 36 grid points and 25 quadrilateral plate 
elements is used to reduce the amount of input data required to illustrate the use of the 
plotting programs. The input data for the programs are presented in appendixes D and E. 
Differences between these input decks occur primarily in the quantities contained in the 
NAMELIST OPTION and NAMELIST PICT. Output listings are similar for both programs, 
therefore, only the output for the contour program is included herein. 

Oblique Orthographic Projections 

The input deck for oblique orthographic projections of the flat plate example is given 
in appendix D. The geometry portion of the input deck includes the entire NASTRAN case 
control and bulk data deck used for analysis of the plate. Option KGEOM = 2 is used to 
read this data; the TRIAEL card is omitted since no triangles are used, and the CQUAD1 
element type is the only one specified on the QUADEL card. Therefore, in this case, all 
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cards in the NASTRAN deck are ignored except GRID cards and CQUAD1 cards. The 
displacement data are punched output from a NASTRAN analysis. Two cards are punched 
for each grid point containing displacements in the x-, y-, and z- directions and rotations 
about the X-, Y-, and Z-axes. The specified format reads two cards at a time but retains 
only the displacements in the x-, y-, and z-directions which are given on the first card 
(x- and y-displacements are zero at all grid points for this example). Two blank cards 
are needed to indicate the termination of the displacement data for this particular format. 

Two NAMELISTS PICT are used to generate oblique orthographic projections. The 
first NAMELIST contains data for generating the deformed shape shown in figure 9(a), 
and contains KODE = 1 which causes the second NAMELIST to be read. In the second 
NAMELIST, the value of KDISP is changed from 1 (as specified in the first NAMELIST) 
to 3 which results in a plot with the displacements represented by vectors as shown in 
figure 9(b). The magnitude of the displacements are arbitrarily exaggerated for visual 
clarity. The second NAMELIST PICT terminates the run with KODE = 0. 

Contour Plot 

The input deck for a contour plot of the flat plate example is given in appendix E. 
Default values are used for most of the NAMELIST parameters. Use of the option, 
KGEOM = 1, to input grid point and element data with a user-specified format is illus- 
trated. GRID cards and CQUAD1 cards from a NASTRAN deck are used, however, any 
specified card formats could have been used. In the case of the GRID cards, only the x- 
and y-coordinates are read. Thus, if nonzero z-coordinates were used to locate the grid 
points, they would be neglected and only the projection would be used. The data to be 
represented by contour lines are the displacements normal to the plate and are given at 
four grid points per card. In NAMELIST PICT, the option to specify labels on contour 
lines on the interior of the plotting surface (see values of NXLAB, XLAB, NYLAB, and 
YLAB) is used since the contour lines do not intersect the boundary. 

The contour plot which was generated by using this input is shown in figure 10. The 
displacements given in the data are symmetric about the two orthogonal principal planes 
through the center of the plate. However, the contour lines generated by the program are 
slightly unsymmetric as shown in figure 10. This slight nonsymmetry is introduced by 
the process of breaking the quadrilaterals into triangles to define the contour surface. 

That is, the triangularization of the surface, which is done automatically, is not quite sym- 
metric since the resulting right triangles face the same direction in all quadrants of the 
plate. This effect would have diminished with a more refined model. Default scaling of 
the data was used to assign a value of 100 to the point of maximum displacement at the 
center of the plate. A listing of the contour-program output for this example is presented 
in appendix F. Titles and headings are provided to indicate the meaning of the printed 
output. 
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CONCLUDING REMARKS 


User and programer documentation is presented for two programs for automatic 
plotting of digital data. One of the programs generates oblique orthographic projections 
of three-dimensional numerical models and the other program generates contour plots of 
data distributed in an arbitrary planar region. The user documentation gives a general 
description of the computational algorithms, user instructions, and complete listings of 
the programs. Several plots are included to illustrate various program options and a 
single example is described in detail to facilitate learning the use of the programs. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., November 18, 1974. 
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APPENDIX A 


INCENTER OF A GENERAL TRIANGLE 

The incenter of a triangle is the center of an inscribed circle as shown in figure 11. 
If the incenters of two similar triangles are coincident, the perpendicular distance between 
corresponding sides of the triangles will be equal. This geometric relationship is desired 
when locating the reference point for triangular elements in an exploded plot. 

The coordinates of the incenter of a general triangle can be derived by using area 
coordinates. Areas of triangles Aj, Aj, and A^ shown in figure 11 are 



By use of equations (8) and (Al), the area coordinates are determined as 



where s = a + b + c. 

The area coordinates can be transformed to Cartesian coordinates by using equa- 
tions (10) to (12) and equation (A2) 


(A3) 
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APPENDIX A - Concluded 


where 
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APPENDIX B 


LISTING OF COMPUTER PROGRAM FOR OBLIQUE 
ORTHOGRAPHIC PROJECTIONS 

An overall flow chart for this program is given in figure 12. The MAIN program is 
used to allocate blank COMMON storage and to call other subroutines necessary to read 
input data and to generate the desired plots. The purpose of each subroutine is described 
in comment cards in the listing. A subroutine called DOCMNT consisting entirely of 
comment cards is included in the program. Subroutine DOCMNT contains (1) a directory 
of selected variables used in the program, (2) user-input instructions, and (3) a descrip- 
tion of plotting subroutines which are required from the Langley Graphic Output System 
library. 


30 



o o o o o n o 


APPENDIX B - Continued 


PROGRAM MAI N! I NPUT=201 , GUTPUT=201 , T AP E5= 1 NPUT, T AP E6=0UT PUT , 
LTAPE9, TAP El 0, TAP £20 = 20 1J 


*** THIS IS THt MAIM PROGRAM WHICH CALLS OTHER SUBROUTINES 
INTEGER NUMPT , XPT »YPT , Z.PT , UPT »VPT , WPT 

COMMON/ CONTRL/ KGEGM.KUATA, KP LUT ,KS YMXY , KSY MXZ , KSYMY Z , NOT AT , XLHT, 
1KMURZ.K VERT, PHI , THETA.PS I , NEWF R , I. SC AL E , PLOT SZ , XORGN , YORGN , 
2PSCALE, KOI SP,DMAG,KOOE 

COMMON/ L I MI TS/ XXMAX , Y YMAX , ZZMAX , XXM I N , YYMI N , ZZMI N , NO MAX, NOMI N, 
1NELMAX, NELMIN 

COMMON/ CORGN/ XOABS , YOABS, XPMAX, XSPACE ,PSIZE 
COMMUN/GLOUP/ ILCCP 
CUMMON/ABLK/ A(3,3) 

COMMON/ SAVEV/ DMAGS , IDMAG 

COMMON/ KOUNT/ NNGOE , NNDEST , NUD 1 SP ,N VO ISP ,NWDI SP 
COMMUN/ VALUES/ NVALLS 
COMMON/ CASE 10/ ICCASE 
COMMON ILL 11) 

DIMENSION DS AV ! 3) 

DIMENSION ABCOJB) 

NAMELIST/ PI Cl/ KHOR Z»KVERT»PH£, THETA.PS I, NEWFR, I SCALE, 

1PL0TSZ , XORGN, YORGN, P SC ALE, NOT AT, KOI SP , IDMAG ,0 MAGS ,KODE, 
2KSYMXY,KSYMXZ ,KSYMYZ, XXMAX, Y YMAX, ZZMAX,XXMIN, 

3YYMIN, ZZM IN, NOMAX, NOMI N , NELMAX ,NE LM IN , XLHT 

*** Tu ZERO NODE AND ELEMENT SUMMATION COUNTERS 


I LOOP = 0 
NNUOE - 0 
XOABS = 0.0 
YOABS = 0.0 
XPMAX =0.0 
300 CONTINUE 

XSTRT = 0.0 
YSTRT = 0.0 
WRI TE ( 6 ,8 ) 

6 FORMAT! 1H1 ) 

*** TO READ TITLE CARD FOR RUN 

RE AO ( 3 , 10 1 ABC 0 

10 FORMAT! BA 10 1 

1 F ! EuF , 3 ) 2222,3333 

2222 CALL PSTUP 
3333 CONT I NO E 

WRITEi6.ll) ABCO 

11 FORMAT! ///,20X,8A10,///» 

CALL IN IT AL 

IF(rPL0T.Nc.A-.0R.KDATA.N£.3) GO TO 111 
REWIND 20 

CALL ME SAGE ! 1 , 24HEN TER I NG GRTHOGC PROGRAM, 2A) 
Call MESAGE ( 1 , 19HSELECT DESIRED MODE, 19) 

CALL MESAGE ( 1, IBHBY PRESSING FN KEY, IB) 

CALL NEAUNKEY) 

IMGuE = NKEY-I 
IF ( I MODE. ED. 0) GC TO 111 
DO 37 J — 1 , I MODE 
KOUNT = 0 
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I F ( NUD ISP.NE.O ) KOUNT - KOUNT+ 1 
IF(NVOISP.NE.O) KOUNT = KOUNT+1 
IF(NWOISP.NE.O) KGUNT = KOUNT+1 
DO 58 I=1,NVALUS 

HEAD 120) IDUM, (DSAV(K) ,K=l, KOUNT) 

58 CONTINUE 
37 CONTINUE 
lit CONTINUE 

HEIGHT = 0.15 

XSTRT = XSTRT+2.0*HEIGHT 

YSTRT = 1.0 

CALL NO TAT El XSTRT, YSTRT .HEIGHT ,ABiCD,9 0.0, 80) 

CALL CALPLT (12. 0,0. 0,-3) 

C 

c *** TO SET POINTERS FOR BLANK COMMON STORAGE ZZZ 
C *** (WITH INTEGER NAMES OF ARRAYS USED IN CALLED SUBROUTINES) 

C 

NUMPT = l 

XPT = NUMPT +NNDEST 

YPT = XPT +NNDEST 

ZPT = YPT+NNDEST 

UPT = ZPT+NNDEST 

I F( NUDI SP.EQ.O) VPT = UPT*-1 

l F ( NUD ISP.NE.Oi VPT = UPT+NNDEST 

IF(NVUlSP.Ey.O) iuPT = VPT+l 

IF(NVDiSP.NE.O) WPT = VPT+NNDEST 

IF(NWDISP.EQ.O) NEND = WPT+1-1 
IFINWOISP.NE.O) NEND = WPT+NNDEST-1 
WRITE(6,15) NEND 

15 FORMAT ( /// ,20X,*8LANK COMMON STORAGE ZZZ REQUIRES AT LEAST *,I6, 

1* LOCATIONS FOR THIS CASE*///) 

IF (KGEOM. EQ. 1) CALL GEOMl 

1( ZZZ (NUMPT) , ZZZ ( XPT ) , ZZ Z l YP T ) , ZZZ ( ZPT ) , LLL ( UP T ) , Z ZZ 1 VPT ) , ZZ Z ( WP T ) ) 

I F ( KGEOM . EQ . 2 ) CALL GEUM2 

l (ZZZ (NUMPT) ,ZZZ( XPT) , ZZZ( YPT) , ZZZ (ZPT) ,Z ZZ ( UPT ) , Z ZZ ( VPT ) , ZZZ ( WP T ) ) 
I F ( KGEOM. EQ .9 ) CALL GEUM9 

1( ZZZ (NUMPT) ,ZZZ(XPT) , ZZZ (YPT) , ZZZ (ZPT) , ZZZ (UPT) , ZZZ (VPT) , ZZZ (WPT) ) 
CALL PN TOUT ( 1 , 

1 ZZZ (NUMPT ), ZZZ (XPT) , ZZZ ( YP T ) , ZZZ.liZPT ) , ZZZ (UPT) ,ZZZ( VPT) ,ZZZ ( WPT ) ) 
600 CONTINUE 

1F( IDCASE.EQ.O) GO TU 650 
KtAUl 5, 10) ABCD 
WRITE(6,ll) ABCD 
650 CONTINUE 
CALL ZEROO 

l (ZZZ (NUMPT) , ZZZ (XPT) , LLLIYPT I , LLH LPT i , ZZZ (UPT) ,ZZZ( VPT) ,ZZZ(WPT)) 
IF(nDATA.EQ.I) CALL DAT A1 

1 (ZZZ (NUMPT) , ZZZ (XPT ) ,ZZZ( YPT) , ZZZ (ZPT ),ZZZ( UPT) ,ZZZ( VPT) ,ZZZ( WPT) ) 
I F ( KuAT A . EQ . 5 ) CALL DATA5 

l (ZZZ (NUMPT) , ZZZ (XPT) , ZZZ (YPT) , LLH aPT. ) , ZZZt UPT ) , ZZZ ( VPT ) ,ZZZ ( WPT ) ) 
IFIKUATA.EQ.9) CALL UATA9 

1 (ZZZ (NUMPT) , ZZZ l XPT ) ,ZZZ( YPT ) , ZZZ ( ZPT ) , ZZZ( UP T) ,ZZZ I VPT) , ZZZ I WP T ) ) 
CALL PnTOUT (2, 

1ZZZ ( NUMPT) , ZZZ(XPT ) ,ZZZ(YPT ), ZZZ (ZPT) ,ZZZ(UPT) , ZZZ ( VPT ) ,ZZZ l WPT ) ) 
700 CONTINUE 

IFtKPLOT.EQ.A.ANO.lLGOP.NE.O) GO TO 6000 
READ(5,P1CT) 

WRI TE (o ,P ICT) 
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oOQu CONTINUE 

IFIKPLQT.EQ.4) CALL CCRTL 
CALL US CALL 

l (ZZZI NU MPT) , Z ZZ I XPT ) » Z ZZ I Y PT ) » ZZZIZPT ) , ZZZIUPT ) , ZZZ( VPT ) ,ZZZl WPT ) ) 
CALL BUUND 

1 (ZZZI NUMPT) , ZZZ IXPT) ,ZZZ( YPT ) , ZZZ I Z PT ) , ZZZI UP T ) , ZZZ I VPT ) , ZZZ I WP T ) ) 
CALL RQTAT 
CALL PcUT 

l (ZZZINUMP1) , ZZZ I XPT), ZZZ I YPT), ZZZIZPT), ZZZIUPT), ZZZ I VPT), ZZZI WPT) > 
IFlNOTAT.fcU.l) CALL NDLET 

i (ZZZI NUMPT) , ZZZ IXPT) ,ZZZ( YPT) , ZZZ IZPT I , ZZZI UPTi , ZZZI VPT) , ZZZ (WPT) ) 

If GKPLGT.EQ.4) CALL CAL PL T I 12 . 0 ,0 .0 , -3 ) 

IFUPLUT.EQ.4) CALL CCKT2 
(LOOP = ILDUP+1 
If (KUOE.EU.O) GC TO 800 
GU TO 1700,600,500) KUOE 
800 CONTINUE 
CALL PSTOP 
END Uf MAIN 
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SUBROUTINE DGCMNT 
C 

C *** THIS SUBROUTINE CONTAINS PROGRAM DOCUMENTATION 


C 

C 

c 

C DESCRIPTION OF INPUT DATA CAROS 

C 
C 
C 

C TITLE CARO - CONTAINS ANY DESIRED ALPHANUMERIC INFORMATION IN COLS. 1-30. 

C NAMELIST OPTIUN - CONTAINS VALUES TO ALLOCATE STORAGE IN BLANK COMMON ZZZ, 
C AND CONTROL VALUES NEEDED BY THE PROGRAM. 

C THE FOLLOWING VALUES ARE INCLUDED 

C NNDEST = ESTIMATE NUMBER OF GRID POINTS TO BE USED. VALUE MUST BE 

C GREATER THAN OR EQUAL TO THE ACTUAL NUMBER OF GRID POINTS. 

C ** DEFAULT = BUG ** 

C NUDISP = 0 FOR NO DISPLACEMENT DATA IN X-DIRECTION. 

C = 1 FOR DATA INCLUDING DISPLACEMENTS IN X-DIRECTION. 

C DEFAULT - 0 ** 

C NVOI3P = 0 FuR NO DISPLACEMENT DATA IN Y-DIRECTION. 

C =1 FOR DATA INCLUDING DISPLACEMENTS IN Y-DIRECTION. 

C *<= DEFAULT = U ** 

C NWDISP = 0 FOR NO DISPLACEMENT DATA IN Z-DIRECTIUN. 

C =1 FOR DATA INCLUDING 01 SPLACEMENTS IN Z-DIRECTION. 

C ** DEFAULT = 0 ** 

C 

C KGtOM SPECIFIES SUBROUTINc AND CORRESPONDING METHOD OF INPUT FOR 

C MODEL GEOMETRY. 

C KGEUM = 1 FOR GRID POINTS AND ELEMENTS READ FROM CARDS WITH USER 

C SPECIFIED FORMAT. 

C = 2 FOR NASTRAN DECK WITH CARD IDENTIFIERS LEFT ADJUSTED 

C AND DATA IN COLUMN WJDTHS OF 8. 

C = V FOR USER SUPPLIED SUBROUTINE - GEOM9 . 

C ** OtFAULT = I ** 

C KOATA SPECIFIES SUBROUTINE AND CORRESPONDING METHOD OF INPUT FOR 

C DISPLACEMENT DATA. 

C KDATA = 1 FOR SUBROUTINE DATA! TO READ IN DISPLACEMENT DATA 

C FROM CARDS WITH USER. SPECIFIED FORMAT. 

C =5 FOR SUBROUTINE DA TAB TO READ IN DISPLACEMENT DATA 

C FKCM TAPE2G. 

C =9 FOR SUBROUTINE DATA9, A USER SUPPLIED SUBROUTINE. 

C ** DEFAULT = l ** 

C NVALUS - USED IF KDAT A = 5 TO SPECIFY THE NUMBER OF SETS OF 

C DISPLACEMENTS TO BE READ FROM TA PE 20. 

C DEFAULT ■= 0 ** 

C IRESEQ = 0 FOR NO RESEQUENCING UF GRID POINT NUMBERS. 

C = 1 TO RESEQUENCE GRID POINT NUMBERS IN SAME ORDER 

C AS THEY ARE INPUT. 

C ** DEFAULT = I ** 

C KPLul SPECIFIES THE TYPE OF OUTPUT DEVICE TO BE USED. 

C KPEUT = 1 FOR CALCOMP. 

C =2 FUR CALCOMP WITH PLOTTING SPEED REDUCED TO USE LEROY PENS. 
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APPENDIX B - Continued 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


= 3 FOR VAR IAN. 

= 4 FOR CRT (OSE CDC250 SCOPES AT IRC) 
** DEFAULT - 1 ** 

XSPACt = SPACE BETWEEN PLOTS IN X-UIRECTIGN, 
** DEFAULT = 10.0 ** 

PS1ZE = PAPER SIZE IN Y-DIRECMON IN INCHES, 
PLOTS TO INSURE THIS DIMENSION IS 
** DEFAULT = 25.0 ** 

IDCASE = 0 FOR NO TITLE CARD PRECEDING 

DECKS OF DISPLACEMENT VALUES. 

= 1 FOR TITLE CARL PRECEDING 

DECKS OF DISPLACEMENT VALUES. 

** DEFAULT = 0 ** 


IN INCHES. 

USED IN SCALING OF 
NOT EXCEEDED. 


MODEL GEOMETRY IS NOW INPUT IN ONE OF THE FOLLOWING FORMS, 

DEPENDING ON THE VALUE OF KGEOM SPECIFIED IN NAMELIST OPTION. 


USE IF KGEOM = 1 


U 

c 

0 

c 

c 

c 

c 


(A) 


( B ) 


ICI 

( 0 ) 


(E) 


(FI 


A SINGLE CARD CUNT AIRING THE WORD FORMAT IN COLUMNS 1-6 AND 
A VARIABLE FORMAT CORRESPONDING TO THE FORMAT OF THE GRID 
POINT CARDS WITH LEFT PARENTHESIS STARTING IN COLUMN 11 
AND UP TO 80 COLUMNS MAY BE USED. 

OtCK OF GRID POINT CARDS. EACH CARD CONTAINS 4 VALUES, GRID 

POINT NUMBER (INTEGER), X-COORDINATE (REAL), Y-COORDINATE 
(REAL) AND Z-CUOROINATE (REAL). THE FORMAT IS SPECIFIED 
IN (A) ABOVE. 

A SINGLE CARD CONTAINING THE WORD ENDGRID IN COLUMNS 1-7. 

A SINGLE CARD CONTAINING THE WORD FORMAT IN COLUMNS 1-6 AND 
A VARIABLE FORMAT CORRESPONDING TO THE FORMAT OF THE 
ELEMENT CARDS WITH LEFT PARENTHESIS STARTING IN COLUMN 11 
AND UP TO 30 COLUMNS MAY 81 E USED. 

DECK OF ELEMENT CAROS. EACH CARO CONTAINS 5 INTEGER FIELDS 
WHICH ARE THE ELEMENT NUMBER, AND GRID POINT NUMBERS AT THE 
VERTICES OF THE ELEMENTS. FOR TRIANGULAR ELEMENTS THE 
LAST INTEGER FIELD MUST BE BLANK OR ZERO. FOR ROD OR BEAM 
ELEMENTS THE LAST TWO INTEGER FIELDS MUST BE BLANK OR ZERO. 
THE FURMAT IS SPECIFIED IN (D) ABOVE. 

A SINGLE CARU CONTAINING THE WORD ENDGEOM 


IN COLUMNS 1-7. 


C USE IF KGEOM = 2 


C 

c 

c 

c 

c 

c 

c 

c 

o 

c 


(A) A SINGLE CARD CONTAINING THE WORD LINEEL IN COLUMNS 1-6 AND 
UP TO NINE NAST.kAN LINEAL ELEMENT CONNECTION NAMES, 

WHICH ARE LEFT-ADJUSTED IN FIELD WIDTHS OF 3, STARTING IN 
COLUMN 9 (CULS. 9-10, L7-24, ..., 73-30). THIS CARD CAN 
BE OMITTED IF LINEAL ELEMENTS ARE NOT USED FOR THE PLOT. 

Id) A SINGLE CARD CONTAINING THE WORD TRIAEL IN COLUMNS 1-6 AND 
UP TO NINE NASTRAN TRIANGULAR ELEMENT CONNECTION NAMES, 
WHICH ARE LEFT-ADJUSTED IN FIELD WIDTHS OF 8, STARTING IN 
COLUMN 9 (COLS. 9-16, L7-24, ...» 73-80). THIS CARD CAN 

BE OMITTED IF TRIANGULAR ELEMENTS ARE NOT USED FOR THE PLOT 
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APPENDIX B - Continued 


c 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 


(CJ A SINGLE CARO CONTAINING THE WORD QUAD EL IN COLUMNS 1-6 AND 
UP TO NINE NASTRAN QUADRILATERAL CONNECTION NAMES, WHICH 
ARE LEFT-ADJUSTED IN FIELD, WIDTHS OF 8, STARTING IN COL. 9 
(COLS. 9-16, 17-24, ...» 73-80). THIS CARD CAN BE OMITTED 
IF QUADRILATERAL ELEMENTS ARE NOT USED FOR THE PLOT. 

(OJ A NASTRAN BULK DATA DECK. ONLY THE GRID CARDS AND THE 

ELEMENT CONNECT ION CARDS WITH NAMES MATCHING THOSE GIVEN ON 
THE LINEEL, TRIAEL, AND QUADEL CARDS WILL BE USED FOR THE 
PLOT. ALL OTHER CARDS IN THE NASTRAN BULK DATA DECK WILL 
BE IGNORED. 

(E) A SINGLE CARD CONTAINING THE WORD ENDGEOM IN COLUMNS 1-7. 


C USE IF KGEUM - 9 


CALL SUBROUTINE GE0M9 WHICH I SI PREPARED BY THE USER TO READ 
GEUMETRY DATA. 


CASE IDENTIFICATION CARD. 


THIS CARD IS OMITTED IF I0CASE=0 IS SPECIFIED IN SOPTION. 

IF PRESENT, THIS CARD CONTAINS ANY DESIRED ALPHANUMERIC 
INFORMATION IN COLS. 1-80. WILL APPEAR BEFORE EACH DATA PLOT. 


DATA TU BE PLOTTED IS NOW INPUT IN ONE OF THE FOLLOWING FORMS, 

DEPENDING ON THE VALUE OF KDATA SPECIFIED IN NAMELIST OPTION. 


USE IF KDATA = 1 

(A) A SINGLE CARD CONTAINING THE WORD FORMAT IN COLUMNS 1-6 AND 
A VARIABLE FORMAT FOR THE DATA CAROS WITH LEFT PARENTHESIS 
STARTING IN COLUMN 11 AND UP TO 80 COLUMNS MAY BE USED. IF 
DISPLACEMENTS ARE INCLUDED FOR MORE THAN ONE GRID POINT PER 
CARD, THE NUMBER OF GRID POINTS PER CARD MUST BE ENTERED 
AS AN INTEGER IN COLUMN 8. 

( B i DECK OF DISPLACEMENT SETS. THERE CAN BE MULTIPLE DISPLACE- 
MENT SETS PER CARD OR THE SET CAN EXTEND TO MORE THAN ONE 
CARD (OFTEN THE CASE WITH NASTRAN PUNCHED OUTPUT) WHICH CAN 
BE HANDLED WITH A FORMAT FOR READING MULTIPLE CARDS. 

A DISPLACEMENT SET FOR EACH GRID POINT IS DEFINED TO 
CONTAIN FROM 2 TO 4 VALUES, A GRID POINT NUMBER AND 
DISPLACEMENTS CORRESPONDING TO NUDISP, NVDISP, OR NWDISP 
EQUAL TU 1. 

(C) BLANK CARD OR CARDS TO END DATA DECK. THE NUMBER OF BLANK 
CAROS MUSI CORRESPOND TO THE NUMBER' OF CARDS READ AT ONE . 
TIME BY THE SPECIFIED VARIABLE FORMAT. 


USE IF KDATA = 5 

READS NVALUS (FROM NAMELIST OPTION) DISPLACEMENT SETS FROM 
TAPE20. EACH DISPLACEMENT SET MUST HAVE BEEN WRITTEN ON TAPE20 
AS AN UNFORMATTED RECORD. 
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APPENDIX B - Continued 


C USfc It KOATA = 9 

CALL SUBROUTINE DAT A9 WHICH IS PREPARED 8Y THE USER TO READ 
DISPLACEMENT DATA. 


NAMELIST PICT - CONTAINS VALUES NEEDED TO GENERATE PLOTS. 
THE FOLLOWING VALUES ARE INCLUDED 


C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


KHURZ = INTEGER DESIGNATING HORIZONTAL AXIS OF VIEWING PLANE, 
WHERE 1 = X» 2= Y , 3=Z. 

** DEFAULT = L *# 

KVERT = INTEGER DESIGNATING VERTICAL AXIS OF VIEWING PLANE, 
WHERE 1 = X, 2= Y , 3= Z . 

** DEFAULT = 2 ** 

PHI = ANGULAR ROTATION OF MODEL ABOUT ITS X-AXIS, IN DEGREES 
(MUST BE TAKEN THIRD). 

** DEFAULT = 0.0 ** 

THETA = ANGULAR ROTATION OF MODEL ABOUT ITS Y-AXIS, IN DEGREES 
(MUST BE TAKEN SECOND). 

** DEFAULT = 0.0 ** 

PSI = ANGULAR ROTATION OF MODEL AdOUT ITS Z-AXIS, IN DEGREES 
(MUST BE TAKEN FIRST). 

** DEFAULT = 0.0 ** 

NEWER = 1 FOR FRAME CHANGE BEFORE PLOT IS MADE. 

(A FRAME CHANGE RESETS THE X-ORGIN PAST PREVIOUS PLOT 
BY X SPACE ANU THE Y-GRGIN AT 0.0). 

N Ev» PR .NE. i FOR NO FRAME CHANGE BEFORE PLOTTING. 

** DEFAULT = 1 ** 

ISCALE = 1 FOR INTERNAL ORGIN LOCATION AND SCALING. 

= 2 FOR USER SPECIFIED ORGIN AND SCALING. 

** DEFAULT = 1 ** 

PlOTSZ = MAXIMUM DIMENSION DESIRED ON COMPLETED PLOT. 

(USED FOR SCALING IF ISCALE = I) 

** DEFAULT = 10. 0 ** 

XORGN = X— LUC AT I CN OF PLOT ORI.GIN (USED IF ISCALE = 2). 

** DEFAULT = 0.0 ** 

YURGN = Y-LOCATICN OF PLCT ORIGIN (USED IF ISCALE = 2). 

*<= DEFAULT = 0.0 ** 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


PSCALE = MODEL SIZE REDUCTION FACTOR, PSCALE = ACTUAL MODEL 
SIZE/DESIRED PLOT SIZE (USED IF ISCALE = 2). 

** DEFAULT = 1.0 *■* 

NUTAT = 0 FOR NO NUMBERING ON PLOTS. 

= I FOR NUMBERING OF GRID POINTS. 

- 2 FOR NUMBERING OF ELEMENTS. 

** DEFAULT = 0 ** 

XL HT = HEIGHT OF INTEGER S SPEC I F I E D BY NOTAT, IN INCHES. 

** DEFAULT = 0.15 ** 

KDISP = 0 FOR UNDEFORMED PLOT. 

= 1 FOR DEFORMEO PLOT. 

= 2 FOR EXPLODED PLOT. 

- 3 FOR DISPLACEMENTS REPRESENTED BY VECTORS. 

** DEFAULT - 0 ** 

IDMAG = 1 FOR DIRECT SCALING OF DATA BY DMAGS. 

= 2 FOR SCALING OF DATA TO A MAX. VALUE OF DMAGS. 

** DEFAULT = 2 ** 
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C 

c 

C 

C 

c 

r* 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


DMAG S = MAGNIFICATION OF DISPLACEMENTS (IF KDISP=I). 

= K EUUCT ION FACTOR OF ELEMENTS (IF KDISP=2). 

** OEFAULT = 1.0 

KSYMXY = I FOR SYMMETRY ABOUT X-Y PLANE. 

** DEFAULT = 0 ** 

KSY/-1XZ = 1 FOR SYMMETRY ABOUT X-Z PLANE. 

** DcFAULT = 0 ** 

RSYMYZ = 1 FOR SYMMETRY ABOUT Y-Z PLANE. 

** DEFAULT = 0 ** 

XXMAXf YYMAX , ZZMAX,XXMIN» YYMIN»ZZMI N LOCATE CUTTING PLANES 
PARALLEL TO PRINCIPAL I X-Y, X-Z, Y-Z) PLANES 
TO LIMIT PLOT. 

** DEFAULT XXMAX=YYMAX=ZZMAX= 1.0E+20 ** 

** DEFAULT XXMIN=YYMIN=ZZMIN=-1.0E + 20 ** 

NUMAX = MAXIMUM GRIU PT. TO BE INCLUDED IN PLOT. 

** DEFAULT = 9999999999 ** 

NO MIN = MINIMUM GRID PT. TO 0E INCLUDED IN PLOT. 

** DEFAULT = 0 

NELMAX = MAXIMUM ELEMENT NUMBER TO BE INCLUDED IN PLOT. 

** DEFAULT = 9999999999 ,** 

NELMIN = MINIMUM ELEMENT NUMBER TO BE I NCLUDtD IN PLOT. 

** DEFAULT = 0 ** 

RODE SPECIFIES CONTROL OPTION AFTER PLOT IS COMPLETE. 

RODE - 0, LAST PLUT, EXIT FROM PROGRAM. 

= I, READ ANOTHER NAMELIST PICT. 

= 2, READ A NEW SET OF DI.SPLACEMENT DATA, INCLUDING A 


C ' CASE IDENTIFICATION CARD IF PRESENT. 

C = 3, READ A COMPLETE NEW SET OF INPUT DATA, 

C INCLUDING A TITLE CARD. 

C ** DEFAULT = 0 ** 


C THE AouVE COMPRISES A COMPLETE BASIC SET OF INPUT DATA IF 
C ROUE = U IN iPICT. FOR ROOt = I, 2, OR 3, ADDITIONAL SECTIONS OF 

C THE BASIC DECK MUST BE REPEATED. THE DECR MUST END WITH 

C NAMELIST JPICT HAVING RODE = 0. 

C 

C 

c 


c 

c 

C DESCRIPTION OF GRAPHICS SUBROUTINES 

C 

C 


v. 

n 

U 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


GRAPHICS 
subRour ine 

CALCUMP 

LEROY 

PSEUUu 

CDC250 


USED BY FOLLOWING 
PLCTTING DEVICE 

CALCCMP 
CALCCMP 
VARI AN 
CRT 


CALLED BY FOLLOWING 
PROGRAM SUBROUTINES 

INITAL 

LNITAL 

INITAL 

INITAL 
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C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

L 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

C 

c 

c 

c 


CALPLT 

CALCOMP .VAR I AN .CRT 

MAIN, CCRT3,PST0P, PLOT ,GARROW 

NOTATE 

CALCOMP, VARIAN, CRT 

MAIN 

number 

CALCOMP, VARI AN, CRT 

PLOT , NDLE T 

NFRAME 

VARIAN 

PLOT 

NEXT 

CRT 

MAIN.CCRT I.CCRT 2.CCRT 3 

MES Au E 

CRT 

MAIN.CCRT 1,CCRT2,CCRT3 

PAR AMS 

CRT 

CCRT I 

KFORMAT 

CRT 

CCRT1 


SUBROUTINE CALCOMP 

PURPOSE THIS IS THE NORMAL MQOt PROCESSOR. THE NECESSARY 

PARAMETERS AND LINKAGE ARE SET UP TO OUTPUT A TAPE 
FOR THE CALCUMP 780/763 0. 0 1 0/0. 005- I NCH STEP PLGTTERi 

USE CALL CALC CMP 

COMMENTS THIS CALL MUST BE GIVEN BEFORE THE FIRST CALL TO A 

PLOTTING ROUTINE. 

SUBROUTINE LEROY 

PURPOSE THE PARAMETERS NECESSARY TO ACCOMODATE PLOTTING WITH THE 

LIQUID INK PEN ARE SET UP BY CALL LEROY. 

USE CALL LEROY 

CUMMENTS THIS CALL SHOULD O.nLY BE USED WITH THE CALCOMP PROCESSOR. 

IN ADDITION TO REDUCING THE SPEED OF THE PLOTTER FOR 
ALL PLOTTING MOVEMENTS. THE NUMBER OF PLOT VECTORS IN ANY 
ANNOTATION IS CONSIDERABLY INCREASED. 

SUBROUTINE PSEUDO 

PURPOSE INITIALIZES PLOT VECTOR FILE FOR VAR I AN PLOTTER 

USE CALL PSEUDO 

SUBROUTINE CDC250 

PURPOSE INITIALIZES CATHODE RAY TUBE CONSOLE. 

USE CALL CDC260 
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c 

c 

c 

c 

c 

C PURPOSE 

C 

C 

c 

c use 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

C COMMENTS 

C 
C 
C 

c 

c 

c 

c 

C PURPOSE 

C 
c 
c 

C USE ■ 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


SUBROUTINE CALPLT 


TO MOVE THE PLOTTER. PEN TO A NEW LOCATION WITH PEN UP 
Of DOWN AND TO SIGNAL THE END OF A JOB SEGMENT BY 
INCREMENTING THE BLOCK ADDRESS NUMBER. 

CALL CALPLT(X,Y, IPEN) 

WHERE 

X* Y ARE THE FLOATING POINT VALUES FOR PEN MOVEMENT. 

IPEN = 2 PEN DOWN 
= 3 PEN UP 

NEGATIVE IPEN WILL ASSIGN X=0, Y=0 
AS THE LOCATION OF THE PEN AFTER MOVING THE 
X , Y (CREATE A NEW REFERENCE POINT) AND 
INCREASE THE BLOCK NUMBER BY ONE. 

ALL X AND Y COORDINATES MUST BE EXPRESSED AS FLOATING 
POINT INCHES (ACTUAL PAGE DIMENSIONS) IN DEFLECTION FROM 
THE ORIGIN. 


SUBROUTINE NOTATE 


TO DRAW ALPHANUMERIC INFORMATION FOR ANNOTATION AND LABEL INI 

AND PROVIDE SPECIAL! CENTERED SYMBOLS FOR ANNOTATION OF 

DATA POINTS. 

CALL NOT ATE (X,Y, HEIGHT , BCD, THETA, N) 

WHERE 

X i Y ARE THE FLOATING POINT PAGE COORDINATES OF THE 

FIRST CHARACTER. FOR ALPHANUMERIC CHARACTERS 
THE COORDINATES OF THE LOWER LEFT-HAND CORNER 
OF THE CHARACTERS ARE SPECIFIED. 

HEIGHT SPECIFIES CHARACTER SIZE AND SPACING IN FLOATING 
POINT INCHES FOR A FULL-SIZE CHARACTER. THE 
WIDTH OF A CHARACTER WILL BE (4/7)*HEIGHT AND THE 
SPACE BETWEEN CHARACTERS IS (2/7 )*HE IGHJ . 

BCD IS THE STRING OF ALPHANUMERIC CHARACTERS TO BE 

DRAWN. 

THETA IS THE ANGLE IN FLOATING POINT DEGREES AT WHICH 
THE INFORMATION IS TO BE DRAWN. 

N IS THE NUMBER OF CHARACTERS, INCLUDING BLANKS, IN 

THE LABEL. 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


PURPOSE 

USE 


COMMENTS 

PURPOSE 

USE 

PURPOSE 

USE 


APPENDIX B - Continued 

SUBROUTINE NUMBER 


TO CONVERT A F LUATI.NG NUMBER TO BCD (EXPRESSED IN 
F FORMAT), AND DRAW- THE RESULTING APL HANUMER IC CHARACTERS. 

CALL NUMBER! X, V, SIZE. FPN, THETA, NI 

WHERE 

X , Y ARE THE COORDINATES IN FLOATING POINT INCHES 

OF THE LEFT LOWER CORNER OF THE FIRST DIGIT 
OF OUTPUT. 

SIZE IS THE HEI.GHT OF THE PLOTTED NUMBER IN FLOATING 

POINT INCHES. 

FPN IS THE FLOATING POINT NUMBER TO BE DRAWN. 

THETA IS THE ANGLE IN FLOATING POINT DEGREES AT WHICH 

THE NUMBER IS TO BE DRAWN. 

N IS THE NUM8ER OF DECIMAL DIGITS TO THE RIGHT 

OF THE DECIMAL POINT FOR OUTPUT. 

N = -1 AND N = 0 BOTH SPECIFY NO DECIMAL PLACES, 
HOWEVER, -1 SUPPRESSES THE DECIMAL POINT. 

THE NUMBER IS RESTRICTED TO A MAXIMUM OF 12 DIGITS. 

THE ROUTINE TRUNCATES THE FLOATING POINT NUMBER AT THE 
REQUIRED DECIMAL PLACE. 


SUBROUTINE NFRAME 


USED BY VARIAN PLOTTER TO ADVANCE PLOTTING FRAME. 
CALL NFRAME 


SUBROUTINE NEXT 


PROVIDES A BREAK PO<INT OR HALT DURING APPLICATION PROGRAM 
EXECUTION. OPERATOR MUST PRESS FUNCTION KEY TO RESUME, 
AND NUMBER OF KEY I.S RETURNEO IN CALLING PARAMETER. 

CALL NEXT ( NKE Y I 

WHERE 

NKEY IS NUMBER OF FUNCTION KEY PRESSED. 
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APPENDIX B - Continued 


c subroutine mesage 

c 

c 

C PURPOSE PRCVIOES THE CAPABILITY TO DISPLAY A MESSAGE ON ThE CRT250. 

C 

C USE CALL MESAGE I l .BCD, N) 

C 

C WHERE 

C 

C I INDICATES INTENSITY OF CHARACTER DISPLAY. 

C 

bCD IS ADDRESS OF ARRAY CONTAINING THE MESSAGE 
IN HOLLERITH FORM. 

N IS THE NUMBER OF CHARACTERS IN THE MESSAGE (LESS THAN 50} 


C 

C 

C 

c 

D 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

G 

c 

c 

c 

c 

c 

c 

c 

c 

c 


SUBROUTINE P ARAMS 


PURPOSE 

USE 


C OMME NTS 


PURPOSE 

USE 


USED TO GENERATE A TABLE OF SYMBOLIC NAMES THAT CAN BE 
ACCESSED USING THE ALPHANUMERIC KEYBOARD ON THE CRT250. 

CALL PARAMS (BCD, VAR ) 

WHERE 


BCD 

VAR 


IS 

IS 


THE HOLLERITH REPRESEN TAT I 
PROGRAM VARIABLE REFERRED 


ON 

OF 

SYMBOL 

IC 

NAME 

TO 

BE 

SYMBOL 

IC 

NAME 


UP TO 3 PAIRS OF VARIABLES MAY BE SPECIFIED IN A SINGLE 
CALL TO PARAMS. TABLE HAS CAPACITY FOR A2 PAIRS. 


SUBROUTINE KFQRMAT 


ALLOWS PROGRAMMER TO CHANGE FORMAT FOR KEYBOARD INPUT. 
CALL KFORMAT<NHdCD} 

WHERE 

N IS THE NUMBER OF CHARACTERS IN BCD. 

H IS REQUIRED. 

BCD IS THE REQUIRED FORMAT ( I A , FA . 2 , A 10, ETC . ) 


RETURN 

END UF OOCMNT 
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APPENDIX B — Continued 


SUBROUTINE CCRT1 
C 

C *** EUR CHANGING VALUES INPUT BY $PICT USING CRT. 

CUMMON/CONTRL/ KGEOM ,KDAT A , KPLOT , KS YM XY , K SYMXZ , KSYM YZ , NOT AT , XLHT , 
IKHORZ » KVERT ,PHl,THETA,PSI,NEwFR, I SC AL E , PLOT SZ , XORGN , YORGN. 
2PSCALE* KOISP,OMAG« RODE 

C UMMUN / L I MI TSZ XXMAX , Y Y M AX , ZZMAX , XXM I N , Y YMI N » ZZM I N, NOMAX , NDM I N, 
INELMAX, NELM IN 
CGMMUN/ SAVLV/ OMAGS , IOMAG 

*** FOR INTEGER CONTROL VALUES 

CALL KFORMATI3HI 1C) 

CALL PARAMS 

Call params i slkhorz ,khorz, slkvert , kvert) 

CALL PARA NS 1 5LNEWFR » NEWFR , 6 L I SCALE , I SCALE) 

CALL PARAMS ( 5L NCTAT , NOT AT »5LKDISP»KDI SP , 5L i DMAG 1 1 DMAG) 

CALL PARAMS! bLKSYMXY.KSYMXY ,6LKSYMXZ, KSY MXZ , 6 LKSY MY Z , KS YMYZ ) 

CALL PARAMS ( 5 L NDM AX , NDM AX , 5 L NDM IN , NDM IN) 

CALL PARAMS (6LNEL MAX, NELM AX ,6LNEL MIN, NELMIN) 

CALl MESAGE I 1 ,32hT0 CHANGE INTEGER CONTROL VALUES, 32) 

CALL MESAGE ( 1 , 13FVAR IABLES ARE, 13) 

CALL MESAGE! I, 27HKHGRZ, KVERT, NEWER, 1SCALE,27) 

CALL MESAGE! 1, I9HN0TAT , KOISP, IDMAG,19) 

CALL MESAGE! 1.22HKSYMXY, KSYMXZ, KSYMY Z ,22 ) 

CALL MESAGE! 1, 28HNCMAX, NUMIN, NELMAX, NELMIN, 28) 

CALL MESAGE t 1 , I 7HANY KEY CONTI NUES, 17) 

CALL NEXT (NKEY ) 

*** FUR FLOATING POINT CONTROL VALUES' 

CALL KF QRMA T ( SHF 1 0. 3 ) 

CALL PARAMS 

CALL PARAMS! 3L YAw.PSI , 4LR0LL , PH I , 5L P LTCH , THET A) 

CALL PARAMS 1 bL PLOT S Z , PL GT SZ , 6LP SC AL E , PSC ALE ) 

C AL L PARAMS ! 5L XORGN , XURGN , 5L YORGN;, YORGN) 

CALL PARAMS 15LDMAGS , OMAGS ) 

CALL PARAMS (A-L XL FT, XLHT ) 

CALL PA RAMS ( 5LXXMAX , XXMAX, 5LXXMIN, XXM IN) 

CALL PARAMS ! 5L Y YMAX , YYMAX, 5L YYM IN , Y YM IN) 

CALL PARAMS t 5L ZZM AX, ZZM AX, 5LZZ MIN, ZZM IN) 

CALL MESAGE ( 1 ,39HT0 CHANGE FLOATING POINT CONTROL VALUES, 39) 

CALL MESAGEll, 1 3FVAR IABLES ARE, 13) 

CALL MESAGE l l , 32HYAW, ROLL, PITCH, PLOTSZ, PSCALE, 32) 

CALL MESAGE (1 , 25HXORGN, YORGN, OMAGS, XLHT, 25) 

CALL MESAGE! 1, 19HXXMAX, YYMAX, ZZM A X , 1 9 ) 

CALL MESAGE! 1,19HXXMIN, YYMIN, ZZMIN,19) 

CALL MESAGE ( 1 , 17HANY KEY CUNT I NUE S , 1 7 ) 

CALL NExTtNKEY) 

RETURN 

ENO OF CCR T 1 
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APPENDIX B - Continued 


SUBROUTINE CCRT2 

*** FOR SELECTING CONTROL OPTION. KODE, AT END OF JOB USING CRT. 


COMMON/ CONTRL/ KG EOM, KD AT A , KPLOT , KS YM XY, K SYMXZ , KSYMYZ , NOT AT , XLHT , 
1KHOkZ,KVERT,PHI , THETA. PSI .NE wFR. LSC AL E , PLOT SZ , XORGN , YORGN, 


2P SC AL E . KO ISP.DMAG, KODE 
10 CONTINUE 

CALL ME SAGE ( 1 . 1H ,11 
CALL MESAGEI 1.32HFN KEY 12 
CALL M £ SAGE ( 1 , 42HFN KEY 24 
CALL MESAGE (1.37 FFN KEY 36 
CALL MESAGEI 1. 18HFN KEY 48 
CALL NEXTINKEY 1 
KOOE = 10 

IFlNKEY.EQ.48) KCOE = 0 
IFINKEY.EQ.12J KOOE = 1 
lFINKEY.Ey.24) KOOE = 2 
IFlNKEY.EQ.36) KCDE = 3 
IFIK0DE.EQ.10) GO TO 10 
RETURN 


TO ALTER EXISTING PLOT, 32) 

TO READ NEW SET OF D I SPL ACEMENT S , 42 ) 
TO READ A COMPLETE NEW CASE, 37) 

ENDS JOB, 18) 


END OF CCR T 2 


SUBROUTINE CCR T 3 

*** RtMINOER TO PUT EOF ON PLOT FILE WHEN USING CRT. 

CALL MESAGEI 1 , 1H ,1) 

CALL MESAGEI 1,45F.LAST REMINDER TO PUT EOF ON PLOT FILE, IF ANY, 45) 
CALL MESAGEI 1 , 1H ,1) 

CALL ME SAGE 1 1 ,'tOF.CO IT AT NEXT PLOT FILE COMPLETE MESSAGE, 40) 

CALL MESAGEI 1 , 1H ,1) 

CALL MESAGEU,17HANY KEY CONTINUES, 17 ) 

CALL NEXTINKEY ) 

CALL CALPLTI 12. 0,0. 0,-3) 

RETURN 

ENU OF CCRT3 


SUBROUTINE PSTOP 

*** TO TERMINATE JOB. 

C 

COMMON/ CONTRL/ KGEGM ,K0 A T A , KPLOT , KS YMXY , KSYMXZ , KSYMY Z , NOT AT , XLHT , 
1KHGRZ,KVERT , PHI , THE T A, P S I , NEWFR , LSC AL E , PLOT SZ , XORGN , YORGN , 
2PSoALE,KDISP,0MAG,K0DE 
CALL CALPLTIO. 0,0. 0,999) 

1FI KPLOT. EQ. 4) CALL CCRT3 
STOP 

END OF PSIUP 
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APPENDIX B — Continued 


SUBROUTINE INITAL 

*** TO SET UP VALUES FOR CONTROL PARAMETERS 

COMMON/ CONTRL/ KGECM , KO AT A , KPLOT , KS YM XY , K SY MXZ , K SYM YZ ,NGT AT , XLH T , 
1KHORZ , N VERT ,PHI , THE T A , P S I , NE WFR , I SC AL E , PLOT SZ , XURGN , YORGN , 
2PSCALE,KDlSP,UMAG,KOOE 

COMMON/ LIM I TS/ XXMAX , YYMAX , ZZMAX , XXM IN , YYMI N , ZZMl N, NOMAX, NUMIN, 
1NELMAX.NELMIN 

COMMUN/ CORGN/ XC A6S , YO ABS , XPMAX , X SP AC E , P S I.ZE 
COMMON/ SAVE V/ OM AG S , IOMAG 

COMMON/ KOUNT/ NNOCE , NNUEST , NUD I SP , N VO I SP , NWD I SP 
COMMON/SttiNCE/ IRESEg 
COMMON/VALUES/ NVALUS 
COMMON/CASE 10/ IOCASE 

NAMELIST/ uPTION/ NNOES T ,NUO IS P , NVD I SP , NWD IS P , 

1KGE0M, KDAT A, NVALUS, I RE SEQ , KPLOT , XSP AC E , P S I Z E , IOCASE 

*** DESCRIPTION OF VALUES IN iOPTION GIVEN - IN SUBROUTINE DOCMNT 


*** TU SET DEFAULT VALUES FOR iCPTION 

NNOEST = 200 
NUOISP = 0 
NVDISP = 0 
NwOlSP = 0 
KGt QM — I 
KDATA = 1 
NVALUS = 0 
IRESEJ = l 
KPLOT = I 
X SP ACE = 10,0 
PSIZE = 25-0 
IOCASE = 0 


*** TO SET DEFAULT VALUES FOR iPICT 


KHORZ = 1 
KVERT = 2 
PHI = 0.0 
THETA = 0.0 

PSl = 0.0 
NEWER = 1 
I SCALE = l 
PLOTSZ = 10. 0 
XURGN = 0.0 
YORGN = 0.0 
PSCALE = 1.0 
NUTAT = 0 
XLH T = 0.15 
KDISP = 0 
l DM AG = 2 
DM AGS = 1.0 
KSYMXY = 0 
KSYMXZ = 0 
KSYMYZ = 0 
XXMAX = 1.0E2O 
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100 
2 00 


YYMAX = 1.0E20 
ZZMAX = 1.0E20 
X XM I N - -1.0E20 
YYMIN = — L . 0t2 0 
ZEMIN ^ -1.0E20 
NOMAX ■= 9099999999 
NDMIN = 0 


NcLMAX = 9999999999 
NELMIN = 0 
NODE = 0 
REA0I5, OPTION) 
IF(cOF,5) 100,200 


CONTINUE 
CALL PS TOP 
CONTINUE 
IFIKPLQT.LE 
IP I KPLOT. EQ 
I F ( KPLO T. EQ 
IF ( KPLOT. EQ 
w K I TE I b 1 0P1 


.2) 

CALL 

CALCOMP 

.2) 

CALL 

LEROY 

.3) 

CALL 

PSEUDO 

.4) 

CALL 

CDC2S0 

ION) 




RETURN 

END OF INITAL 


Continued 
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APPENDIX B - Continued 


SUBROUTINE GEOM i ( NUMPT »XPT»YPT»ZPT»UPT*VPT,WPT) 

*** TU READ GRID POINT INFORMATION AND ELEMENTS 
*** FROM CARDS WITH USER SPECIFIED FORMAT. 

COMMON/CONTRL/ KGEOM , KDAT A , KPLOT , KS YMXY , KSYMXZ, K S YMYZ , NOT AT , XLHT , 
lKHURZjRVERT *PHI ,THETA,PSI,NEWFR,i SC AL E , PLOT SZ , XORGN , YORGN , 

2 P SC ALE , KDISP»DMAG t KODE 

COMMON/ KUUNT/ NNODE » NNDEST *NUDI SP , N VD I SP , NWD I SP 
COMMON/SEQNCE/ I RES EQ 

DIMENSION NUMPI ( 1 ) »XPT( 1) , YPT (1) » ZP T ( L) ,UPT ( l) ,VPT( 1) ,WPT(1 ) 
DIMENSION ABCD ( 8 ) 

DIMENSION FORMTI7) 

DIMENSION NUDE ( A ) 

REw IND 10 

TEST1 = 10HENDGRID 
TEST2 = lOHtNDGECM 
TESTA = 10HFORMAT 

*** TO READ GRID INFORMATION 

100 CONTINUE 

RtAOI 5 , 10 I ABCD 
10 FURMAT ! 8 A 1 0 ) 

I F t AdCu(l).EQ.TESTl) GO TO 1000 
IFIABCD(l).Ew.TESTA) 50,60 
50 CONTINUE 

DECUDEI80,55,ABCD} FORMT 
55 FORMAT l 10X, /AlO) 

GO TO 100 
60 CONTINUE 

NNODE = NNUDE+l 

DECODE! BO, FORMT , ABCD) NUMPT (NNODE ), XP T! NNODE ) , YPT ( NNODE ) 

1 ,ZPT(NNODE) 

GO TO 100 
1000 CONTINUE 
C 

C *** TO READ ELEMENT CONNECTION INFORMATION 
C 

200 CONTINUE 

kEAD( 5, 10 i ABCD 
IF(AbCD(lJ.E0.TEST2) 275,285 

275 CONTINUE 

END FILE 10 
GO TO 2000 
285 CONTINUE 

IFIABCDI 1 ). EQ. TESTA) 250,260 
250 CONTINUE 

DECODE! 80, 55, ABCD) FORMT 
GO 10 200 
260 CONTINUE 
IDUM = 1 

uECODc! 80, FORMT , ABCD) NUMEL , NODE l , NOD E2 , NODE 3 , NODEA 
IF ( ikESEQ.NE. 1 ) GO TU 700 
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APPENDIX B — Continued 


**# TO RENUMBER ELEMENT NOOES 

NODE(i) = NOOEl 
NODE I 2 ) = N0DE2 
NODE ( 3) = NODE 3 
N0DEI41 = N0DE4 
DO 500 1=1,4 

I F ( NODE ( 1 1 . EQ . 0 } GO TO 550 
DO 510 J= 1 1 NNOD E 
IF (NODE (1 I.EQ.NUMPTl J) ) 511,510 
511 CONTINUE 

NODE! I) = J 
GO TO 500 
510 CONTINUE 
500 CONTINUE 
550 CONTINUE 

NOD El = NODE ( 1 ) 

N0Ut2 = NODE l 2 ) 

NODE3 = NODE 131 
N0DE4 = NODE ( 4 I 
700 CONTINUE 

CALL RECOUT ( 10 , 1 , 0 , NUMEL , NODE l , NODE 2, NODE 3, N0DE4J 
GO TO 200 
2000 CONTINUE 
RETURN 
END OF GEOM1 
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APPENDIX B - Continued 


SUBRUU FINE GEUM2 INUMPT , XPT , YPT , ZP T , UP T , VP T , W P T ) 

*** FOR INPUT OF NASTRAN DECK TO DESCRIBE GEOMETRY. 

COMMON/CONTRL/ KGECM ,KDAT A , KPLOT, KS YMXY , K. SYMXZ, KSYMY Z , NOT AT , XL HT , 
IKHDRZ,KVERT ,PHI , T HET A , P S I , NEWER , I.SC AL E , PLOT SZ , XORGN , YORGN , 
2PSCALE*KDISP»DMAG, NODE 
COMMON/ SEQNCE/ IRESEQ 

COMMUN/KUUNT/ NNOUE »NNUEST »NUOI SP , N VO I SP , NWD TSP 
D I MENS I UN NUMPT! 1) ,XPT ( 1) ,YPT 1 1) , ZPT! 1) ,UPT (1) ,VPTi 1) ,WPT( L ) 
DIMENSION ABCDI8) 

DIMENSION NODE (4) 

DIMENSION L INEEL 19) , TR l AEL ! 9) , QUADE L ! 9) 

10 FORMAT ( 8A 1 0 ) 

REWIND 9 
TESTu = I OH 
TESTl = 10HGR1D 
TEST2 = 10HLINEEL 
TEST3 = 10HTRIAEL 
TEST4 = lOHOUAUEL 
TESTS = 1 OHENOGECM 
DO 15 1=1,9 
LINEELII) = 10H 
TRI AELI I) = 10 H 
UUADELl I) = I OH 
15 CONTINUE 
800 CONTINUE 

RE AD ( 5 , 1 0 ) A8CU 
DECODE! 80, 50, A8CD) W0RD1 
50 FORMAT! A8) 

IFIWORDI. Eg. TESTO) GO TO 800 
IFIWORDI. EQ. TESTl) GO TO 100 
If (W0RD1.EU.TEST2) 60,61 

60 DEC0DE!80,lll,ABCG) IL INEEL ! I ) ,1 = 1,9) 

111 FORMAT t 6X.9A8) 

GU TU 800 

61 CONTINUE 

IFIWORDI. ED. TEST3) 62,63 

62 UEC0DEt80,lll,ABCD) ! TR IAEL (I) , 1 = 1 ,9 ) 

GO TO 800 

63 CONTINUE 

I F ( WORD 1. Eg. TESTA) 64,65 

64 DECGOElSO.lil, ABCD) (QU ADEL ! I ) , 1= 1, 9) 

GU TO 800 

65 CONTINUE 

DO 70 1=1,9 

IFIWORDI. EQ.LINEELl I)) GO TO 200 
IFlWURDl. ED. TRIAEL! D) GO TO 300 
IFCWURDl.Ew.gUADEH I ) ) GO TO 400 
70 CONTINUE 

I F ( WORD 1. EG. TESTS) GO TO 2000 
GO TO 800 

*** TO READ GRID CAROS ^ 

100 CONTINUE 

NNODE = NnODE+1 

UEC0DE<80, 101, ABCD) NUMPT ! NNODE ) ,XP T < NNODE ) , YPT ( NNODE ) 

1 , ZPT ( NNODE) 
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101 F0RMATi8X,Ao,ax,3F6.0) 

CALL I R 1 TE I NUMPT ( NNGOE I ) 

GU Tu 300 

*** TO RE AO CARDS CONTAINING ELEMENTS WITH 2 GRID POINTS 

200 CONTINUE 

TtSTlOu = 1 OHP LOT EL 
lEIWOKOl.EG.TESTlOOJ GO TO 250 
DECODE I 80 ,201 ,A8CCI NEL ,NDl,ND2 

201 FORMAT l 8X ,A8, 8X.2A81 
LALL IRITEINELI 
CALL IRITEIND1J 
CALL 1RITEIN02) 

NU3 = 0 
NU4 = 0 

CALL RE COOT l 9 , l , 0 , N EL , NO 1 , ND2 , ND 3 , ND4 I 
GO TO dOQ 

*** TO RcAO PLOTEL CAROS 

250 CONTINUE 

UECOUtl 80.251, A8CUI NEL,ND1,N02 

251 FORMAT! 8X,3A8) 

CALL IRITEINELI 
CALL I R I T E ( NO l I 

CALL I R I T E( N02 I 
NG3 = 0 
N04 = 0 

CALL RECOuT ( 9» 1 , 0 , NEL , ND1 , ND2 , NO 3 , ND41 
GU TO a 00 

*** TO RcAO CAROS CONTAINING ELEMENTS WITH 3 GRID POINTS 

300 CUNTINUE 

DECODE! 80,301, A 8CD I NEL ,ND1 , ND2 ,ND3 

301 FORMAT !8X,A8,6X, 3A8) 

CALL IRITEINELI 
CALL IRITEIND1I 
CALL IR l TE! N02 I 
CALL I R I T E ! N03 I 

N04 = 0 

CALL RECOUT !9, L ,0 , NEL, ND1 , N02 , ND3:, ND4 1 
GO TO 300 

*** TO READ CARDS CONTAINING ELEMENTS WITH 4 GRID POINTS 

400 CONTINUE 

OECOOEl 80,401 , ABCO) NEL ,NDl ,ND2 , ND3 , N04 

401 FORMAT l 8X ,A8, 8X,4A8) 

CALL IRITEINELI 
CALL IRITEIND1J 
CALL IR ITE! ND2 1 
CALL I R I T E ! NO 3 I 
CALL IR ITE! NU4 1 

CALL RECOUT! 9, 1 , 0 , NEL , N01 , N02 , ND-3 , ND4 I 
GO TO 800 
2000 CONTINUE 

END FILE 9 
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*** TO RENUMBER ELEMENT NODES ANO WRITE ON TAPE 10 

REWIND 9 
REWIND 10 
600 CONTINUE 

CALL RECINl 9, 1 , 5 , NUMEL , NODE l , N0DE.2 , NUDE3 , NODE 4) 

I F ( EOF , 9 ) 450,451 
451 CONTINUE 

IFIIRESEQ.NE.1) GO TO 700 

NODE III = NUDE 1 
N0DEI2) = N0DE2 
NODE ( 3) = NODE 3 
NODE ( 4 ) - N0DE4 
DO 500 1-1,4 

1 F ( NUDE (II.EQ.O) GO TO 550 
OU 510 J= 1 , NNOD E 
IFINOOEI U.EQ.NUMPTt J)) 511,510 
511 CONTINUE 

NODE (II = J 
GO TO 500 
510 CONTINUE 
5 00 CONT I NU E 
550 CONTINUE 

NOD El - NODE ( 1 J 
NODE 2 = NODE-121 
NOD E 3 = NODE ( 3 ) 

N0DE4 = NOD E ( 4 J 
700 CONTINUE 

CALL RECUUT ( 10 , 1 , 0 , NUMEL , NODE 1 , NODE 2, NODE 3, NODE 4 J 
GO TU 600 
450 CONTINUE 

END FILE 10 

RET URN 

END OF GEUM2 
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SUBROUTINE IRITE(NUM) 

*** TU RIGHT ADJUST INTEGERS IN A FIELD WIDTH OF EIGHT 

DIMENSION N(8J 
LANK = 1H 

DECOuE! 3, 1,NUM) NSAVE 

1 FORMAT! IB) 

OtCODE! 3 , 2, NUM) (N( I ) *1 = 1 » 3 i 

2 FORMAT (BAD 
DU 10 1=1, B 
II = 9-1 

IFINIII J.NE.LANK) GO TO 20 
10 CONTINUE 

20 NUM = N SAVE/ I 10**1 S— II ) ) 

RETURN 
END OF IRITE 


SUBROUT INE GE0M9 ( NUMPT , XP T , YPT , ZP T , UP T , VP T , WP T) 

*** USER SUPPLIED GEOMETRY INPUT SUBROUTINE. 

COMMON/ CONTRL/ KGEOM , KO AT A , KPLQT, KS YM XY , K SYMXZ , K S YMYZ , NOT AT , XLH T , 
1KH0RZ tKVERT ,PHI , THETA, PS I , NEWER , I SCAL E , PLOT SZ , XORGN , YORGN , 
2PSCALE, KOISP,OMAG,KOOE 

COMMON/KUUNT / NNOOE , NNDEST ,NUOISP»NVDISP»NWDI SP 

DIMENSI ON NUMPT! 1 ) , XPT! 1) ,YPT (ll.ZPTI 1 ) , UPT < 1 J , VP T ! 1) ,WPT! 1 ) 

*** INSERT ROUTINE HERE 

RETURN 

END OF GE0M9 
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SUUROUT INE BOUND ( NUMPT , XPT, YPT, ZPT ,UPT, YPT, WPT) 

C r 

c *#* TO DETERMINE MAXIMUM DIMENSIONAL LIMITS OH BODY FOR USE 
C IN SCALING PLOTS 

C 

COMMON/ CUN TRL/ KGEOM ,KDAT A , KPLUT , KS YM XY , KSY MXZ , KSYMYZ , NOT AT , XLHT , 
IKHORZ.RVERT ,PHl , THET A , P S I , NEwFR, I.SC AL E , PLOT SZ , XORGN , YORGN , 
2PSCALE*NDISP,0MAG«K0DE 

COMMON/ LIMITS/ XX MAX , Y Y MAX , ZZM AX , XX MI N, YYM I N , ZZM 1 N, NDMAX, NDM I N , 
1NELMAX, NELM IN 

COMMON/ XYZLI M/ XY ZMAX ( 3 ) , XY ZM I Nt 3 ) 

COMMON/ KOUNT/ NNGDE ,NNUEST , NUOI SP ,NVD ISP.NWDl.SP 

DIMENSION NUMPT I 1} ,XPT (I) ,YPT( 1) rZPTl l ) ,UPT < 1 } , VP T ( 1 ) , WPT ( I ) 

DIMENSION NODE ( 4 ) 

DO 5 1=1,3 
X Y ZM I N ( L) = +1 .UE20 
XYZMAXl I I = -1.0E20 
3 CONTINUE 
R c W I N D 10 
100 CONTINUE 

CALL RE C INI 10, 1 , 5 , NUMEL , NUDE 1 , NODiEZ , NODE 3 , N0DE4 ) 

I FI EOF , 10) 1000,200 
200 CONTINUE 

IFINUMEL. LT.NELMIN. OR. NUMEL. GT.NELMAX) GO TO 100 

NODE ID = NODE 1 

N0DEI2J = NODc 2 

NODE I 3 J = NUDE 3 

N0UEI4) = NUDE 4 

DU 10 1=1,4 

NO = NUDE! I j 

IFINODEI U.EQ.Oi GO TO 15 

IFINUMPTINO) .LT.NOMIN. OR. NUMPT ( NO ). GT . NDMAX) GO TO 100 
10 CONTINUE 
is continue 
DO- 20 1=1,4 

I F I NODE 1 1 ) , EQ . 0 ) GO TO 25 
ND = NODE I I ) 

IFIXPTINDKGT.XXMAX) GO TO 20 
IF(XPTINO).LT.XXMIN) GO TO 20 
IF I YPTI ND) .GT.YYMAX) GO TO 20 . 

IF IYPT IND) .LT.YYMINJ GO TO 20 
IF! ZPT l NO) .GT.ZZMAX) GO TO 20 
IFIZPT IND) .LT.ZZMINi GO TO 20 
IF! XPTI ND } .GT. XYZMAXl 11 ) XYZMAXtl) = XPTIND) 

IF IXPTINDl.LT. XYZMIN(l) ) XYZMINU) = XPTIND) 

IF! YPTINO) .GT.XYZMAXI2) ) XYZMAXI2) = YPTIND) 

IF I YPTI ND) . LT. XYZMINI2 ) ) XYZMIN12) = YPTIND) 

IFIZPT (ND).GT. XYZMAX I 3 ) ) XYZMAXI3) = ZPTIND) 

IFIZPT I ND ) . L T . XYZMI N I 3 ) ) XYZMINI3) = ZPTIND) 

20 CONTINUE 
25 CONTINUE 
GO TO 100 
1000 CONTINUE 

DO 300 1=1,3 

IF! I .ED. 1 . AND. KSYMYZ.NE . 1) GO TO 300 
IF! I.EQ.2.AND.KSYMXZ.NE.1) GU TO 300 
IF! I.EQ. 3 .AND. KSYMXY.NE.l) GU TO 300 
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XYZdlG = Ab S l XY ZMAX ( I ) ) 

IFIAdSlXYZMlNl II ) .GT.XYZBIG) XYZ6 IG = AB S I X YZM INI IJ ) 
XYXMAX(I) = XYZeIG 
XYZrtlNII) = -XYZBIG 
300 CONTINUE 
KbTUkN 

cNO UF BOUND 


SUB ROUT INE ZEROO ( NUMPT , XPT , YP T , ZP T , UPT , VP T , WPT I 

C 

C *** INITIALIZES ALL DISPLACEMENTS TO ZERO. 

C 

COMMGN/KOUNT/ NNODE * NNDES T , NUD i SP , NVD I S P , AIM D I SP 
DIMENSION NUMPT I 1 J , XP T (IJ , YP T I l J , ZPTI U ,UPT (II , VPT( LI ,WPTI 1 J 
IF! NUOISP.EQ.O) GO TO 200 
DU 150 1=1* NUD ISP 
UP1UI = 0.0 
150 CONTINUE 
200 CUNTINUE 

IF ( NVOISP.EQ.O) GO TO 300 
DO 250 1= l f NVD I SP 
VPT (1) =0.0 
250 CONTINUE 
300 CONTINUE 

IF1 NWDISP.EQ.O) GO TO 400 
DO 350 1 = 1 * NWU ISP 
WPTI U =0.0 
350 CONTINUE 
400 CONTINUE 
RE I URN 

END OF ZEROD 
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SUBROUTINE DATA1!NUMPT,XPI, YPT,ZPT,UPT,VPT,WPT) 

*** Tii READ DISPLACEMENT DATA FROM CARDS 
*** WITH USER SPECIFIED FORMAT. 

COMMON/ CONTRL/ KGEOM f KO AT A , KPLOT , KSYMXV , KSY MXZ , KSYMYZ , NOT AT , XL HT , 
1KH0RZ,KVERT ,PH1 , THETA,PSI ,NEWFR, I.SCAL E , PLOT SZ , XORGN , YURGN , 

ZP SCALE, KQISP,DMAG,KQOE 

COMMON/KOUNT/ NNCDE , NNDEST , NUD ISP , N VD I SP , NWD I SP 

cummun/ Sconce/ ireseu 

DIMENSION NUMPT I 1 ) ,XPT I 1J ,YPT III ,;ZPT ( 1 ) ,UPT < 1 ) , VP T (1) , WPT U ) 

DIMENSION ABCDI8) 

DIMENSION FORMTI7) 

DIMENSION ISAV! 10) ,DSAV(10,3) 

10 FORMAT l 8A10) 

*** TO READ DISPLACEMENT INFORMATION 

TEST = oHFORMAT 
READ ( 5 * 10) ABC D 
DECODE!80,45,ABC0) WORD , K VALU 
AS FORMAT ( A6 , IX, I 1 ) 

IFtKVALU.Eg.O) KVALU = 1 
I F < WORD .Ey.TEST) 300,200 
2 00 WRITE(6,2U) 

20 FORMAT! 1H1,/// ,20X,* SORRY, FORMAT FOR DATA NOT GIVEN*) 

S I UP 

300 CONTINUE 

DECOO E ! BO , 3 U , ABCO ) FGRMT J 

50 FORMAT l IOX, 7A10 ) 
lOu CONTINUE 
KOUNT = 0 

IrlNUuISP.NE.O) KOUNT = KQUNTa 1 
IF t NVOISP.NE.O ) KOUNT = KOUNTaI 
IFINwuiSP.NE.O) KOUNT - KOUNTaI 

RE ADI 5, FGRMT) ( I S AV ! K) , I OSAV ( K, I ) , 1=1 , KOUNT ) , K= 1 , KVALU) 

DO AGO K=i, KVALU 
IDUM = ISAV(K) 

IF! IOUM.EQ.O) GO TO 1000 
IF! IRESEQ.NE.l ) GO TO 700 

*** FOR KE SEQUENCED GRIO POINTS 

00 500 J=1»NNG0E 
I F l NUMPT ( J) .EU.IUUM) 501,500 
501 CONTINUE 
kJUNT = 1 

I F ! NUD ISP.NE.O) 610,620 
610 UPT(J) = OSAV! K , KOUNT) 

KOUNT = KOUNTaI 
620 CONTINUE 

IF! NVOISP.NE.O) 630,640 
o3U VPItJ) = OSAV! K ,KCUNT) 

KOUNT = KOUNTaI 
o40 CONTINUE 

IF(NWUISP.NE.O) 650,660 
o50 WP T l J ) = OSAV! K , KOUNT) 

660 CJhT 1 NUE 
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GO TO 550 
500 CONTINUE 
550 CONTINUE 
oU TO 400 
700 CONTINUE 

*** FOR NO RESEQUENCE OF GRID POINTS 
KOUNT = I 

IF(NUDISP.NE.O) 1610,1620 
1610 UPTIIDUM) = DS AV ( K » KOUNT J 
KOUNT = KOUNT *■ 1 
1620 CONTINUE 

IFINVDISP.NE.O) 1630,1640 
1 o30 VPT(IDUM) = DS AV ( K , KOUNT ) 

KOUNT = KOUNT + 1 
1640 CONTINUE 

IFINrtDISP.NE.O I 1650,1560 
1650 WPT(IDUN) = OSAVIK, KOUNT) 

1660 CONTINUE 
400 CONTINUE 
GO TO 100 
1000 CONTINUE 
RETURN 

END UF DA I A 1 
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SUB ROUT liMt OAT A5INUMPT, XPT, YPT, ZP T,UPT, VPT, WPT) 

C 

*** TO RLaD DISPLACEMENT DATA FROM TAPt20. 

COMMON/ CONTRL/ KGEUM ,KO AT A, KPLOT ,*.S YMXY , KS YMXZ , KSYMYZ , NOT AT ,XLHT, 
1KHURZ»kVEkT,PHI , THE T A, P SI ,N twFR, I SC ALE , PLOT SZ , XORGN , YORGN , 

2PSCALE, KOISP.DMAG.KGDE 

CUMMGN/ KOUNT / NNOOE , NNOEST , NUO I SP . N VO I SP » NWDI.SP 
COMMON/ SEONCE/ IKESEU 
COMMON/ VALUES/ NVALUS 

DIMENSION NO MPT ( 1 ) , XPT I 1 1 , YPT ( 1 ) ,ZPT( 1 1 , UP T ( 1 ) , VPT ( L ) , WPT ( I ) 

DIMENSION OSAVI3) 

DO LO 1 = 1, NVALUS 
KOUNT = 0 

IF! NODI SP.NE.O ) KUUNT = KOUNT-M 
IFINVDl SP.NE.O) KOUNT = KOUNT+1 
IFINwDlSP.NE.O) KUUNT = KOUNT*-! 

REAU120) IDUM, (OSAVIK) ,K=L, KOUNT) 

KUUNT = l 

I FINUUi SP.NE.O) 610,620 
610 UUUM = OS AV ( KUUNT ) 

KUUNT = KOUNT+1 
o20 CONTINUE 

IFINV01SP.NE.GJ 630,640 
630 VOUM = OS A V I KOUN T J 
KUUNT = KUUNT+1 
690 CONTINUE 

IF! NrtDI SP.NE.O ) 630,660 
630 iwOUM = USA V (KOUNT ) 

6b0 CONTINUE 

I F I IRESEU.NE.1I GO TU 550 
OU 600 J=l, NNOOE 
IF (NUMPTI J) .Ew. IOUMI 501,500 
501 CONTINUE 
10UM = J 
GO TO 550 
500 CONTINUE 
550 CONTINUE 

IF(NODISP.NE.O) UPT(IDUM) = UDUM 
IFINVDI SP.NE.O) VPTIIOUM) = VDUM 
IFlNWDISP.NE.O) WPT ( IDUM ) = WDUM 
10 CONTINUE 

RE TURN 
END OF OAT A 5 

SUBROUT IN E DAT A9 INUMPT , XPT, YPT.ZPT , UPT , VPT, WPT) 

*** USER SUPPLIED DISPLACEMENT INPUT SUBROUTINE. 

COMMON/ CUNT RL/ KGEOM ,KD AT A , KPL OT ,,KS YMX Y , KSVMXZ , K SYMYZ ,NOT AT , XLHT, 
lKHGkZ»KVEKT,PHi,THErA*PSI,N£wFR»ISCALEfPLOTSZ,XORGN»YORGN, 

2PSCALE,K0lSP,0MAG,KUDt 

CUMMQN/RGUNT/ NNOuE , NNOEST , NUD l SP ,NVO 1 SP ,NWO I SP 
DIMENSION NOMPT(l),XPTIl),YPT( l ) , ZPT ( 1 ) , UPT (I) , VPTI 1 ) , WPT II ) 

*** INSERT ROUTINE HERE 

RETURN 

END OF DAT A 9 
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SU6R0UI I ME PNTOU T ! I OUT , NUMPT , X PT ,Y PT , ZPT , UPT , VP T, WP T ) 

*** FUR PRINTED OUTPUT UF INFORMATION IN BLANK COMMON - ZZZ 

C UMMJN / KOUN T/ NNODE , NND E ST , NUD I SP » M VD I S P , NWOI SP 

DIMENSION NUMPT I II ,XPTI 1) ,YPT(I) , ZPTI I) , UP T I U » VPT< 1) ,WPT( 1 ) 

GU TO 11000,20001 IQUT 
1000 CONTINUE 

*** FOR OUTPUT OF GEOMETRY INFORMATION 
wR1TE(o,16) 

lo FORMAT! ///, 5X, *GRID POINT INFORMATION*,///) 

WRITE! o,l 71 

17 FORMAT ( 5X,*RESEQUENCED*,4X,*USER INPUT*/ 

15X,*0RID POINT* ,5X,*GRID POINT*/ 

2 5 X , *NUM BE R* , 9X , *N UMBER* f l3X,*X*,14X,*Y*,l4X,*Z*//> 

00 30 1=1, NNODE 

WRI IE <6 , IB) I , NUMPT < I I , XP T ( 1 1 , YP T ( i I , ZP T ( I I 

18 FORMAT ( 2X,I 10, 5X , 1 10, 3X ,3E1 5 .4 ) 

30 CONTINUE 

wRITE<6,19) 

1 9 FORMAT!///, 5X,*ELEMENT INFORMA T I OiN - WITH RESEQUENCED GRID POINTS 
1*///) 

WRI TE (o ,21) 

21 FORMAT! 5X,*REScQUfcNCED*,4X,*USER INPUT* , 19X ,*GR I D POINTS*/ 
13X,*ELEMENT*,8X»*ELEMENT*/ 

28X,*NUM8ER*,9X,*NUM8ER*,13X,*1*,9.X, *2* , 9X ,* 3* , 9X, *4* / / ) 

REWIND 10 

1 = 0 

35 CONTINUE 

1 = 1 + 1 

CALL RtC I N! 10 , 1 * S , NUMEL » NODE 1 , N0DE2 * N0DE3 » NOD E4 I 
1 F t EOF , 10 I 36,37 

36 RET URN 

37 CONTINUE 

WRI Tc! t, ,22) I, NUMEL , NODE 1 , NODE 2 , NUD E 3 , NOU E4 

22 FORMAT 12X, 110, 5X»I10,4X, 4110 
GO TU 3 5 

2000 CONTINUE 

*** FUR OUTPUT UF DISPLACEMENT DATA 
*R1 TE(b,210) 

210 FORMAT! /// , 5X , *D l SPLACEMENT S TO BE PLOTTED*,///) 

W R1 TE ( 6 , 1 7) 

DU 230 1=1, NNODE 
U = 0.0 

IF! NUD I SP.NE.O) L = UP T ( I ) 

V,= 0.0 

lPtNVUl SP.NE.O) V = VPT(I) 

W = Q 
w = 0.0 

IF(NWUISP.NE.O) W = WPT(I) 

WRITE(u,ld) I , NUMPT! I ) ,U, V, W 
230 CONTINUE 
RET URN 

end of pntout 
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SUBRUUr INE DSCALE(NUMPT,XPT , YP T , Z-PT , UPT , VP T , WPT ) 

*** THIS SUBROUTINE DETERMINES THE SCALE FACTOR FOR DISPLACEMENTS 

CuMMON/CONTRL/ KG ECM , KD AT A , KP LGT , KS YMXY , K SY MX Z , K SYM YZ , NOT AT , XL HT , 
1 KHORZ,kVERT,PHI,THETA,PSI ,NEwFR, I.SC ALE , PLOT SZ ,XORGN,YORGN, 
UPSCALE, KD ISP, DMAG, NODE 
COMMON/ SAVEV/ DMAGS » IDMAG 

COMMON/ KOUNT/ NNCDE , NNDEST , NUO I SP ,N VD I SP , NWD I SP 
DIMENSION NUMPTI1) ,XPT( 1) , YPTC li ,ZPT( 1 1 » UPT ( 1 ) » VP T( 1 ) ,WPTC 1 ) 
IF(KDI3P.Ew.0.QR.KDlSP.EU.2) GO TO 10 
GO TO (10,201 IDMAG 
10 CONTINUE 

DMAG ■= DMAGS 
GO TO 30 
20 CONTINUE 
UMAX = 0.0 
DO 100 1=1, NNCDE 
I F ( NO DISP.EU.OJ GO TO 500 

IF ( ABSl UPTC I ) I .GT.DMAX) DMAX = ABiS (UPT ( I ) ) 

500 CONTINUE 

IFl NVUISP.EQ.O) GO TO 501 

1F( ABS(VPT( I)) .GT.DMAXJ DMAX = ABiS(VPKI)) 

501 CONTINUE 

1FINUDISP.EQ.U) GO TO 502 

IF( ABS( WPT( I)) .GT.DMAX) DMAX = ABS(WPTd)) 

502 CONTINUE 
100 CONTINUE 

DMAG = UMAGS/DMAX 
30 CONTINUE 
RETURN 

END OF DSCALE 
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SUBROUTINE PLU T (NUMPT »XPT » Y PT , ZP I ,UPT , VP T , WPT J 
*** FOR GENERATING PLOTS. 

COMMON/ CONTRL/ KGEOM ,KOATA , KPLOT, KS YM XY, K SYMXZ , KSYMYZ , NOT AT , XLHT , 
1KHORZ ,KVERT ,PHI , THETA, PS I , NEWF R, I.SC AL E , PLOT S Z , XORGN , YORGN , 
UPSCALE, KOI SP, DMA G.KOOE 

COMMON/ LI MI TS/ XXMAX , Y YMAX , ZZMAX , XXMI N, Y YM I N , ZZM1 N, NDMAX, NDMIN, 

IN EL MAX, NELM IN 

„ COMMON/XYZLIM/ X YZMAXI 3 ) , XY ZM I N t 3) 

COMMON/ CORGN/ XG AbS , YOABS , XPMAX, XSP AC E , P S I ZE , 

COMMON/GLOOP/ I LOOP 
COMMON/ ABLK/ A (3, 3) 

CGMMON/KOUNT/ NNODE .NNDEST , NUDI SP ,N VD l SP , NWDi SP 
COMMON/PDELS/ OELX,DELY 

DIMENSION NUMPT (1) , XPT l II ,YPT(i),ZPT(l» .UPTIIi ,VPTt I) ,WPT(I) 
DIMENSION NODE (4),X(4),Y(4) ,ZI4),XDISP(4J ,YDISP (4) , ZDISPI 4) , 
1XR0TI4J ,YRUT(4) 

DIMENSION TEST (3) 

DIMENSION ABCDI3) 

1 FORMAT I 8A1G) 

2 FORMAT! IX , 8 A 10 I 
C 

C *** TO MAKE ALL GRID POINT NUMBERS NEGATIVE- 
C 

Du 50 I = 1 , NNODE 
NUMPT! I ) = -NUMPT! I J 
50 CONTINUE 

PI = 3.1415926 
XMQVE = 0.0 

IFINEWFR.EQ.IJ XMOVE= XPMAX+XSPACE 

YMOVE = -YUAdS 

CALL CALPuT ( X MOV E » YMOV E » -3 ) 

XUA8S = XOABS+XMUVE 
YUABS = YUABS+YMCVE 
GO TO ! 701, 701 , 703, 701 I KPLOT 
701 CONTINUE 
GU TO 710 
703 CONTINUE 

IF! NEnFR.EQ.l) CALL NFRAME 
710 CONTINUE 
DEL X = 0.0 

dely - 0.0 

IFIISCALE.LQ.D CALL XYSCAL 
CALL CALPLT l XORGN , YORGN ,— 3 ) 

XOABS = XOA BS+ XORGN 
YOABS = YUA8S+Y0RGN 
XSH1FT = 0.0 
YSH1FT = 0.0 
Z SHIFT = 0.0 
XPMAX = -1.GE20 

0 

C *** LOOPS TO ACCOUNT FOR SYMMETRY 
C 

ZSIGN = +1.0 
DO 300 11=1,2 

IF! II. EQ. 2. AND. KSYMXY.NE.il GO TO 3 00 
IF! I I .EQ.2.AND.KSYMXY.EG. 1) ZSIGN: = -1.0 
YSIGN = +1.0 
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DO 5 10 JJ=1,2 

I F t JJ-EQ.2.AN0.KSYMXZ.NE. 1) GO TO 510 
IFIJJ.EQ.2.AND.KSYMXZ.EQ.L) YSIGN = -L.0 
XSIGN = +1.0 
DO 520 KK=1 , 2 

IFIKK.EQ.2.AND.KSYMYZ.NE. 1) GO TO* 520 
IFIKK.EQ.2.AN0.KSYMYZ.EQ.1 ) XSIGN = -1.0 
C 

*** TO DETERMINE PROJECTED COORDINATES OF ELEMENTS 

REWIND LO 
TOO CONTINUE 

CALL RE C INI 10, 1 , 5 , NUMEL , NODE 1 » N00E2 , N0DE3 ,N00E4 ) 
IFIE0F.10) 1000,200 
200 CONTINUE 

IFINUMEL.lT.NELMIN.OR.NUMEL.GT.NELMAXJ GO TO 100 
NODE ( I ) = NODE 1 
NODE I 2 ) = NODE 2 
NODE! 3 ) = NODE 3 
NODE 14 j = NUDE4 
DO 10 1=1,4 
NO = NODE ( I ) 

IFINUOEI I I.EO.Oi GO TO U 

*** TO MAKE GRID POINT NUMBERS CONNECTED BY ELEMENTS POSITIVE 
NUMPT I NO j = I A8S I NUMPT I NOJ J 

IFI NUMPT ( NO j .LT.fsCMlN. OR. NUMPT (NO).GT .NDMAX ) GO TO 100 
NENO = I 

10 CONTINUE 

11 CONTINUE 
I = KHURZ 
J = KVERT 
DO 20 N=1 , NEND 
ND = NODE IN) 


IFI XPI l ND) 

.GT 

. XXMAX ) 

GO 

TO 

100 

IFIXPTINO) 

• LT 

. XXM IN ) 

GO 

TO 

100 

IF I YPT I ND) 

• GT 

. YYMAX) 

GO 

TO 

100 

IF I YPT I ND) 

.LT 

. YYMIN) 

GO 

TO 

100 

IF ( ZPT(ND) 

• GT 

. ZZMAX ) 

GO 

TO 

100 

1FIZPTIND) 

.LT 

. ZZM IN ) 

GO 

TO 

100 

XDISPIN) = 

0. 

0 




YDISPlN) = 

0. 

0 




ZDISPIN) = 

0. 

0 





IFIKDISP.EQ.l. AND.NUDISP.NE.O) XDISPIN) = UPT.(ND) 
IFIKOISP.EQ.I.AND.NVDISP.NE.OJ YDISPINi = VPTIND) 
IFIKDISP.EQ.L. AND. NWDiSP.NE.OJ ZDISPIN) = WPTIND) 
XIN) = XSIGN*! XPTINDJ+XOISPIN) *0M AG+X SHI F T ) / PSC AL E 
YIN) = Y S I GN* I YPT INDI+YDISPINJ *DMAG+Y SHIFT ) /PSC ALE 
ZIN) = ZSIGN*IZPTINU)+ZDISPIN)*DMAG+ZSHIFT)/PSCALE 
20 CONTINUE 

I FI KUlSP.tQ.2) CALL XP LCD I NENO ,X , Y, Z J 
XCENT =0.0 
YCENT = 0.0 
DO 25 N= 1 , NEND 

XROTIN) = All, 1)*XIN)+A(1,2)*YIN)+AII ,3)*Z(N) 
YROTIN) = AIJ, 1 ) *XINJ+AIJ,2)*Y (NJ+AIJ , 3 ) *Z ( N ) 

XCENT = XCENT+XRCiT I N) 

YCENT = YCENT+YRCTIN) 
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XROT(N) = XROT ( N i +OELX 
YRQT(N) = Y ROT ( N J +0E LY 
l F ( XROT (N ) . GT . XPMAX 1 XPMAX = XROT IN) 

25 CONTINUE 

XCENT = XCENT/FLUAT(NEND)-(6.0/7.0)*XLHT 
YCENT = YCENT/FLGAT (NENOl-XLHT/2.0 
XCENT = XCENT+OElX 
YCENT = YCENUUELY 
AL = NO MEL 

I F ( NOT A T . Ew • 2 1 CALL NUMBER { XC ENT, YC ENT , XLHT , AL , 0. 0, -1 1 

*** TO PLOT ELEMENTS 

CALL CALPLT IXROT ( 1 1 »YROT( 1) ,3) 

GO 30 N=2,NENU 

CALL CALPlT (XROT (N) ,YRUT ( N ) ,2) 

30 CONTINUE 

CALL CALPLT UROT(NEND) ,YRGT ( NEND) * 3 ) 

IFINEN0.GT.2) 35,36 

35 CALL CALPLT (XROT ( 1) ,YRGT (1) ,2) 

CALL CALPLT (XROT ( l) , YROT ( 11,3) 

3 6 CONTINUE 

go ru too 

1000 CONTINUE 

IFIK01SP.EQ.31 600,650 
600 CONTINUE 

*** TO PLOT VECTORS AT GRID POINTS 

UO 601 NQ = l , NNUOE 
IFiNUMPT(NU) .LE.O) GO TO 601 

I FI NUMPT ( NO) »LT . NCMIN.OR. NUMPT I NDi) . GT . NDMAX 1 GO TO 601 

IFlXPTINDl.GT.XYZMAXm > GO TO 601 

IFtXPKNOJ .LT.XYZMINJ li) GO TO 601 

L F l YPT l NO ) . GT . XY ZMAX (21) GO TO 601 

IFIYPT l NO 1 . L T . X Y ZM l N ( 2 ) ) GO TO 601 

IFIZPUND).GT.XYZMAX(3) ) GO TO 601 

IF(ZPT(ND).LT.XYZMIN<31) GO TO 601 

X(l) = XSIGNU XPUNOUXSHIFT) /PSCALE 

Y(l) = YS IGN*(YPT (NOl+YSHIFT 1 / PSC AL E 

Zll) = Z S I GN* ( Z P T ( NU 1 + Z SHI F T 1 / PSC AL E 

XDISPti 1 - 0.0 

Y01SPI 11 = 0.0 

ZOI SP ( 1 1 = 0.0 

I F ( NUu t SP . NE . 0 ) XOISPll) = UP UNO:) 

IF(NVDISP.NE.G) YD1SP11) = VPTiND) 

1FINWOISP.NE.01 ZDISP(l) = WPTCND.) 

X ( 2 1 = XS IGN* (XPT (NOl+XOI SP ( 1 1 *OMAG+X SHI F T 1 /PSC AL E 
Y ( 2 1 = YSIGNU YPUNOI+YUI SP ( 1) *DMAG + XSHIFT) /PSCALE 
Z 1 2 1 = ZSIGN^IZPT (NCUZDISPI 1) *0MAG+ZSH1FT ) / PSC AL E 
l = KHJRZ 
J * KVERT 
00 o05 N= 1,2 

XKUT(M) = A([,1J*X(N) + AU'Z)*Y(N)+A(I,3)*Z(N) 
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605 


601 
o50 
520 
510 
5 00 


YROTIN) = Al J, 1 l*X(Ni+A(J,2 1*YIN)+A(J ,3)*ZI N) 

XKOTIN) = XROTtN)+DELX 
YROTIN) = Y ROT I N ) +DELY 
CONT I NOE 
XARW = 0.06 
YARW = XARW / 3.0 

CALL GARRuwIXRUTll) ,YROTU) .XR0TI2) , YROT 1 2 1 , 1 ,XARW, YARW 

CONTINUE 

CONTINUE 

CONTINUE 

CONTI NUE 

CONTINUE 

RE TURN 

ENO OF PLOT 


SUBROUTINE RuT A f 


SETS UP COEFFICIENTS OF ROTATION MATRIX 

COMMON/ CONT RL/ KG ECM ,KQAT A , KPLOT ,KS YMXY , KSYMXZ, KSYMY l , NOT A T ,XLHT, 
1KH0RZ,KVERT,PHI,THETA,PSI , NEWER, I.SC ALE, PLOT SZ ,XORGN, YORGN, 
2PSCALE, KOI SP,OMAG,KOO£ 

COMMuN/ABLn/ A ( 3,3) 

PI = 3. 1^1592o53b 
S 1 N PH I = SIN(PHI*PI/180.0) 

CUSPHl = CUS(PHI*Pi/ldO.O) 

SINTHc = SINI THE T A#P I / 180.01 
COS THE = CGSlThETA*Pi/180.0) 

SINPSI - SINI PSI *PI /180.0I 
CUSPSl = C0S(PSI*PI/160.0) 

A I 1 , 1 ) = COSTHE’i'COSPSI 

All, 21 = CUSPS I*SINT.HE*SINPHI-S1NPS1*C0SPHI 
All, 3) = SINTHE*C0SPHl*C0SPSH-SINPHI#SINPSl 
AI2.11 ~ S 1NPS I *CCSTHE 

AI2.21 = S1NTHE*SINPHI*SINPSI + C0SIPHI*C0SPSI 

A I 2 , 3 1 = SINTHE*COSPHI*SINPSI-SINPhI*COSPSI 

A l 3 , 1 1 ^ -SINTHE 

A ( 3 , 2 1 = COST HE*S INPHI 

A I 3 , 3 1 = COST ric^CCSPHI 

RETURN 

end of rotat 
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SUBROUTINE XYSCAL 

*** TO DETERMINE SCALE FACTOR FOR MODEL GEOMETRY. 

COMMON/CONTRL/ KGEOM ,KD AT A , KPLOT ,XS YM XY, KSYMXZ , KSYMYZ , NOTAT , XLHT , 
1 KHORZ.K VERT , PHI , THE T A , P S I , NEWFR, I.SC AL E , PLOT SZ , XORGN , YORGN , 
2PSCALE,KUISP,OMAG,KODE 
COMMON/XYZLIM/ XYZMAXI 3 ) , XYZMINI 3 ) 

COMMON/ CORGN/ XOABS , YQABS , XPMAX , XSRAC E ,PSI ZE 
COMMON/ ABLK/ AI3.3) 

COMMON/ PUEL S/ UELX.DELY 
I = KHURZ 
J = K.VERT 
UMAX = 0.0 
DO 5 N= 1 » 3 

VDOM = ABS( XYZMAXINI-XYZMININ)) 

IFl VDUM.GT.DMAX) UMAX = VDUM 
5 CONTINUE 

PSCALE = UMAX/PLOTSZ 
DO 10 L=l,2 
DO 10 M=1 , 2 

00 10 N=I,2 

X = XYZMINI li 

IF( L.EQ.2) X = XYZMAXI I i 

Y = XYZMINI2) 

IFIM. EU.2) Y = XYZMAXI 2) 

Z = XYZMINI 3 ) 

IFIN. EU.2) Z = XYZMAXI 3 ) 

XROT = A( I , l J*X+A( I , 2)*Y+A( 1 , 3)*Z 
YROT = AIU, 1J*X+A(J,2J*Y*-A(J»3J*Z 
IFIL*M*N.EQ.U 20,30 
20 CONTINUE 

XRMIN - XROT 
XRMAX = XROT 
YRMIN = YROT 
YRMAX = YROT 
30 CONTINUE 

IF! XROT.GT.XRMAXi XRMAX = XROT 
IF IXROT .LT. XRMIN ) XRMIN = XROT 
I F I YROT. GT. YRMAX) YRMAX = YROT 
IFIYRUT .LT. YRMIN) YRMIN = YROT 
10 CONTINUE 

YR = ABSI YRMAX— YRMIN) 

1 F ( YR/PSCALE.OT .PSIZE) PSCALE = Y:R/PS1ZE 
XRMAX = XRMAX/ PSCALE 

YRMAX = YRMAX/ PSCALE 
XRMIN - XRMIN/PSCALE 
YRMIN = YRMIN/ PSCALE 
DEL X = -XRMIN 
DELY = -YRMIN 
XORGN = 0.0 
YORGN = 0.0 
RETURN 

END OF XYSCAL 
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SUBROUTINE XPLOD(NEND,X,Y,Z) 

c 

C *** FUR GENERATING EXPLODED PLOTS. 

C 

COMMON/ CONTRL/ KGEOM.KDATA, KPLOT , ,K S YM XY, K SYMXZ , K SYMYZ , NOT AT , XLH T , 
1KHORZ ,KVERT ,PHl , THETA, P SI , NEWFR , I, SC AL E , PLOT SZ , XORGN , Y ORGN , 
2PSCALE,KDISP,OMAG,KQDE 
DIMENSION X(4l ,Y(4) ,Z(4» 

C 

C *** TO CALCULATE THE INCENTER OF TRIANGLES 
C 

I FIN END. EQ. 31 10,20 
10 CONTINUE 

A = SURTC (X(2)-X(3H**2 + IY(2»-Y(3))**2<-(Z(2)-Z13» )**2) 

B = SQRT({.X(l)-XI3))**2*(YCn-Yt3i)**2+(Zm-Z(3))**2) 

C = SuRTl(X(ll-X(2))**2HYm-Y(2l)**2HZlll-il2ll**2) 

AC 1 = A/(A*-B + C) 

AC2 = B/(A+B<-C) 

AC 3 = C / I A+ B + C ) 

XOC ■= AC1*XUH-AC2*X(2)+AC3*X< 3) 

YOC = ACI*YU)+AC2<‘Y12)+AC3*Y13) 

ZUC = ACl*Z(l)+AC2*ZI2)+AC3*Z(3) 

GO TO 190 
20 CONTINUE 
C 

C *** TO CALCULATE THE CENTROID OF RODS, BARS, AND QUADS 
C 

XOC - 0.0 
YOC =0.0 
ZOC = 0.0 
DO 100 I=1,NEND 

xuc = xoc+xm 

YOC = YOC+YI I ) 

ZOC = ZOCtZlI) 

100 CONTINUE 

XOC = XOC /FLOAT INEND) 

YOC = YOC/FLOAT (NtNDi 
ZUC = ZGC/FLOAT I NEND1 
190 CONTINUE 
C 

C *** TO REDUCE THE SIZE OF THE ELEMENT 
C 

DO 200 1=1, NENU 
X ( I ) = XI I) *DMAG 

Y ( I J = YII) *DMAG 
Z(I» = Z( I ) *DMAG 
200 CONTINUE 
C 

C *** TO CALCULATE THE CENTROID OF THE REDUCED ELEMENT 
C 

XRC = XOC*DMAG 
YRC = YUC*DMAG 
ZRC = ZOC*DMAG 
C 

C *** SHIFT CORNERS OF OkIGINAL AND REDUCED TO MAKE CENTROIDS MATCH 
C 

DO 400 I=1,NEND 

xm = xi i j ♦ i xoc-xrc ) 
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400 


❖ * £ 


iOOO 

5000 


rm = Y( I 1 + ( YOC— YRC ) 

n n = zi i > *■ ( zgc-zrc ) 

CONTINUE 

RETURN 

END OR XPLOD 


SUB ROUT IN E G ARROW IXl.Yl , X2 , Y2 , NC , XH t A D , YHE AD ) 

TO DRAW ARROWS FROM Xl*Yi TO X2, Y:2. 

DEN = SORT I l X2-X L)**2+1Y2-Y 1 ) * *2) 
IFIDEN.EQ.O.O) GO TO 5000 
C = (X1-X21/DEN 
S = (Y1-Y21/DEN 
CALL CALPLT(Xi,Yl,3i 
CAlL CALPLT! X2 ,Y2, 2) 

IFtNC.LT.ii GO TO 1000 
XA = X2 + I C*XHEAD-S*YHEAO) 

YA = Y2*< S=S=XFIEAU+C*YHEAD) 

CAlL CALPLT(XA,YA,2) 

1FINC.LT. 2) GO TO 1000 
X 8 = X2 + 1 C«XHEAU-S*t-YHEAD) ) 

YB = Y2M S*XHEA0 + C* I-YFIEAD) 1 
CALL CaLPLT (XB,Y8, 21 
IFINC.LT. 3) GO TC 1000 
CALL CALPLT(X2,Y2,2» 

IFINC.LT. 4) GO TO 1000 
XC = A2+(-S*YHEA0i 
YC = Y2+-I +C*YHEAD) 

CALL CALPLT (XCiYC»2J 
IFINC.LT. 5) Gu TO 1000 
XD = X2+I - S*(-THEAO) J 
YD = Y2 + I + C# I — Yh E AD ) i 
CALL CALPLT (XD, YD, 2 ) 

CONTINUE 

CALL CALPLT (X2,Y2, 31 

CONTINUE 

RETURN 

END OF G ARROW 
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SUB ROUT IN E NOL E T ( NUMPT , XP T , YP T , ZP T, UP T, VPT, WRT) 

*** FOR ANNOTATING GRID PUINT NUMBERS ON PLOTS. 

COMMON/ CUNTRL/ KGEOM ,KD AT A« KPLOT , KS YMXY , KSYMXZ , K SYM Y Z , NOT AT ,XLHT, 
iKriORZ > X VERT » PH l » THETA, P SI , N EWFR , I.SC AL E , PLOT SZ , XGRGN , YORGN , 
2PSCALE»ROISP»QMAG* KOOE 

COMMON/LIMITS/ XXMAX , Y YMAX , ZZMAX, XXMI N , Y YMI N , ZZ MI N , NDMAX, NOMI N, 
INELMAX, NELMIN 

COMMON/XYZLIM/ XYZMAX(3),XYZMIN< 3) 

COMMON/ ABLk/ A ( 3 , 3 ) 

COMMON/ KOUNT/ NNOOE , NNDEST , NUU l SP , N VO I SP , NWD l SP 
COMMON/Pl/ElS/ DELX.DELY 

01 MEN SI ON NUMPT (1) ,XPT I 1) ,Y PT ( 1 ) , ZPT (I) ,UPT 11) , VPT< 1 ) ,WPT U) 

II = KHORZ 

JJ = KVERI 

XSHIFT = 0.0 

YSHIFT = 0.0 

ZSHIFT = 0.0 

DO 300 1=1, NNOOE 

IF < NUMPT UI.LE.OJ GO TO 500 

IFINUMPTI 11 .LI .NOMIN.OR.NUMPTI I I.GT.NDMAX) GO TO 500 

IFIXPTI II .GT.XYZMAXllJ i GO TO 500 

IFIXPTI ll.LT.XYZMINt 1) I GO TO 500 

IFIYPTI I1.GT.XYZMAX12I) GO TO 500 

IFIYPT ( I) .LT.XYZMINI2) I GO TO 500 

IFIZPTI Il.GT.XYZMAX(i) I GO TO 500 

IF (ZPTl II .LT.XYZMINI 31 1 GO TO 500 

X = ( XP f( I) +XSHIFT) /PSCALE 

Y = (YPTC II +YSHIFD/PSCALE 

Z = (ZPT( 1) +ZSHIFTI/PSCALE 

X RUT = A< II , LJ*X + A( II,2I*Y + A<II,3)*Z 

YKOT = A(JJ,1J#X+A(UJ,2)*Y+A(JJ,3I*Z 

XL = XROT + XLHT / 2 .0 

YL = YRGT+XLHT/2.0 

XL = XL+OELX 

YL = YL+OELY 

AL = NUMPT ( l I 

CALL NUMBER (XL ,YL,XLHT ,AL ,0. 0,-1) 

500 CONTINUE 
RETURN 

END OF NQLET 



APPENDIX C 


LISTING OF COMPUTER PROGRAM FOR CONTOUR PLOTS 

An overall flow chart for this program is given in figure 13. The MAIN program is 
used to allocate blank COMMON storage and to call other subroutines necessary to read 
input data and to generate the desired plots. The purpose of each subroutine is described 
in comment cards in the listing. A subroutine called DOCMNT consisting entirely of 
comment cards is included in the program. Subroutine DOCMNT contains (1) a directory 
of selected variables used in the program, (2) user-input instructions, and (3) a descrip- 
tion of plotting subroutines which are required from the Langley Graphic Output System 
library. 
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PROGRAM MAI M INF CT=20l ,0UTPUT=201 ,TAPE5= INPUT ,TAPE6=0UT PUT , 

IT AP E20= 201) 

C 

C *** THIS IS MAIN PROGRAM WHICH CALLS OTHER SUBROUTINES 
C 

c 

c *** USER UUCUMENTAT ION IS GIVEN IN SUBROUTINE DGCMNT 
C 

INTEGER NUMPT, XFT , Y P T , WP T ,R ACPI , NUMEL , NODE 1 , NOUE2 , NODE3 , NOD E4 , 
1CENTR 

COMMON/ KOUNT/ NNODE , NEL MT , NNDEST , NELE SI 

COMMON/ CUNT RL/ INFORjKGECM, KDATAt K.S I GN , KPLOT . XURGN , YORGN, 

IP SC ALE, I SCALE, wMAG, I CN T K , XLhT , NXL AB , X L AB I 10 1 , NYL AB, YL AB 1 10 ) 
COMMON/ CRUDE/ KCCE.NPLOT 
CUMMON/ VALUES/ NVALUS 
COMMON/ CASE I D/ 1CCASE 
COMMON/ SAVE V/ WMAGS.ICNTRS 

CCMMuN/SlZEXY/ XM IN I M , XMAX l M , YhiN 1M , Y MAX I M, XSPACE , Y SPACE 
COMMON aZZII) 

DIMENSION ASCO ( B ) • BCOC 2 i 

NAM tL I S T/PICT/ N FLU T »XORGN, YORGN, PSCALE, I SC ALE , WMAG S , I CNTRS » XLHT , 
1NXLAB,XLAB, NYL AB , YL AB, KOUE 
500 CONTINUE 
I LOOP = 0 

*** TO ZERO NODE ANO ELEMENT COUNTERS AND SET TEST VALUES 

NNODE = 0 
NELMT = 0 
W R 1 TE(6»B) 

B FORMAT l.lHi) 

*** TO READ TITLE CARC FOR RUN 

Rc AO I 5 , 101 ABCO 

10 FORMAT! 8A101 
1 F ( EOF , 51 2222,3333 

2222 CALL PSTUP 
3333 CONTINUE 

WRI TE (o « 111 ABCD 

11 FORMA TC///, 20X,BA10,///) 

CALL INITAL 

IFIKPlOT.NE. A. OR. KUATA.NE. 51 GO TO 111 
REWIND 20 

CALL MESAGEC 1 ,2hHENTERI NG CONTOUR PROGRAM, 24) 

CALL MESAGEI t, 19FSELECT DESIRED M.UDE.19) 

CALL MESAGEI 1, 18HEY PRESSING FN KEY, 18) 

CALL NEXTINKEY) 
l MODE = NKEY-1 
I F I I MODE . EQ .0 ) GO TC 111 
DO 57 J= l » I MODE 
DO SB 1=1, NVALUS 
READ I 20 ) IOUM , WOCM 
58 CONTINUE 
57 CONTINUE 
ill CONTINUE 

HEIGHT = 0.15 

XSTRT = 2.0*HE IGFT + 2.0 
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YSTRT = 1.0 

CALL NO TaT£(XSTRT, YSTRT, HEIGHT , ABCD , 90. 0, BO ) 

CALL NFRAME 
C 

C *** TO SET POINTERS FCR BLANK. COMMON STORAGE ZZZ 
C *** (WITH INTEGER NAMES OF ARRAYS USED IN CALLED SUBROUTINES) 

C 

NUMPT = 1 

APT - rtUMP T + NNOE S I 
YPT = XPT + NNDEST 
WPT = YPT+NNDEST 
RAOPT = WPT^NNDEST 
NUMEL = RAUPT+NNCEST 
NUDE 1 = NUMEL+NELEST 
NUDE 2 = NOUEUNELEST 
NODES — N0UE2+NELEST 
NUDEA = NUDt3+NEL£ST 
CENTR — NUDE4*NEL£ST 
NEnD = CENTR+NELEST-1 
WRITE(o,15) NEND 

15 FORMAT! ///,20X,*8LANK COMMON STORAGE III REuUIRES AT LEAST *,I6, 

1* LOCATIONS Fur THIS CASE*///! 

*** CALL SUBROUTINES FCR INPUT OF CONTROL VALUES AND GEOMETRY 

IFIKwtUM. tO.l) CALL GEGMi! 

Lli.ll NUMPT ) .ZZZIXPT) .ZZZIYPT) tZZH WP T ) .ZZZIRAOPT), 

2ZZZI NUMEL), ZZZI NC DEI) .ZZZIN0DE2) .ZZZI NODE 3 ) till INODES) ,ZZZ( CENTR) ) 
I F ( KGeUM. EG . 2 ) CALL GECM2( 

1ZZZ l NUMPT ) , ZZZ ( XP T ) .ZZZ(YPF) .ZZZI WPT) .ZZZIRADPT ) , 

2 ZZZ I NUMEL ) » ZZZINODEl) »ZZZlNG0E2).ZZZ( NODE 3 i till ( NOD E4) , III I CENT R.) ) 
IF( KGEUM. tw. 9) CALL GECM9I 
Lilli NUMPT ), ZZZIXPT) .ZZZIYPT ),ZZZC WPT) .ZZZIRADPT), 

ILLLi NUMEL ).ZZ2(N0DEl),ZZZ(NGUE2)-,ZZZt N0UE3) , III ( NODE'. ) , III I CENTR ) ) 
CALL XY SC AL I 

Lilli NUMPT) ,ZZlI XPT) .ZZZIYPT ) , ZZZ! WPT) .ZZZIRADPT) , 

1 1 III NUMEL ),ZZZ(NCC£l),ZZZ(NCU£2)«ZZZiNGUE3) till I N0DE4) , ZZZI CENTR) ) 
CALL ENCLOSI 

1ZZZI NUMPT ) .ZZZIXPT ) , ZZZIYPT), ZZZ I. WPT) .ZZl(RADPT) , 

Zlll I NUM EL ) , ZZZ I NODE 1) .ZZZlNCDEZ) , III ( NUDES) .ZZZINUDEv) , ZZZ ( CENT R1 ) 
CALL PNTOUTIl, 

1ZZZ I NUMP T ) .ZZZIXPT) , ZZZI YPT ),ZZZI WPT) .ZZZIRADPT) , 

lllli NUMEL ),ZZZ(NCCEi).ZZZ( NODE 2) * ZZZ I N0DE3) tllli N0DE4) .ZZZ(CENTR) ) 
oOO CONTINUE 

IF t IDCASE.EU.O ) GC TO 650 
READ! 6,10) ABC D 
wKITtlo.il) ABCU 
HEIGHT = 0.15 
XSTRT ^ 2.0*hE IGHT + 5.0 
YSTRT ■= 1.0 

CALL NUT AT E l XSTRT .YSTRT .HEIGHT,. ABCD, 9 0.0, 80) 

CALL NFRAME 
650 CONTINUE 
Call zeruw 

1 ( ZZZI NUMPT) till l WPT) tllli NUMEL), III i CENTR)) 

1 F I RUAT A. EQ. 1 ) CALC CAT A1 
II ZZZ (NUMPT) .ZlZIWPT) ,Z ZZI NUMEL ), ZZZ IC ENTR ) ) 

I FIKCATA.E0.5) CALL DATA5 
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l(ZZZtNUMPT) ,ZZZUPT) ,ZZZ(NUREL)» ZZZ J CENTR ) ) 

If UDaTA.EO. 9) CALL QATA9 
KZZZ(NUHPT) ,ZZZ4hPTJ ,ZZ Zi NU REL J , ZZZ { CENT R )) 

CALL PN TOUT 12* 

iZZZlNUMPT ) , ZZZ (XPT i ,ZZZ(YPT),ZZZIWPT) ,ZZZ(RADPT), 

2ZZZI iNUMfcL ) , Z^ZtNODED, ZZZ( N00E2) , ZZZl NUOE3) ,ZZZ(N0DE4) , ZZZ ( CENTR) ) 
*** CALL SUBROUTINES FOR PLOTTING GEOMETRIC LAYOUTS 
700 CUNT I Nut 

IFInPLOT.EO. 4. ANC.ILOGP.NE.G) GO TO 6000 
RE AO l b * PLOT ) 

wRITEtofPICTJ 
6000 CONTINUE 

I F( RPL0T.c0.41 CALL CCRTi 

IF (NPLUT. EO.4.0R .NPLQT . EU. 5 ) 211, 212 

2 1 1 CONTINUE 

CALL SCALEWIZZZUPTJ , ZZZ (CENTR J i 
LABEL - dH WMAG = 

ENC0Gt( 20,3.33, BCC) LABEL, OMAG 
333 FORMAT ( Ad ,E12.4) 

HE i Grl f - 0.15 
XSTRT = X SP ACE 
YSTRT = 2.0*HL(GHI 

CALL NO TATE (XSTRT, YSTRT, HE 1GHT , BCD, 0. 0,20 ) 

212 CONTINUE' 

YSPACE * 0.0 

1 F ( Y SPA C£.LT.5.C*H EIGHT ) YSPACE = 5.0*HEIGHT 
XORGN = -XMINIM/PSCALE4-XSPACE 
YURGN ^ FLOAT ( KSIGM*(-YPIMMI / PSCALE+YSPACE 
call calplt (xorgn, yorgn, -3) 

1 F ( NPLU r . LE . 31 2C 1 , 202 

201 CONTINUE 

CALL LA YUUT ( NPLUT , 

IaZZ(NUMPT) ,ZZZ(XP f ) ,ZZZ( YPT),ZZZ(WPT) ,4ZZ(RA0PT) , 

2ZZZ t NUMEL ) , ZZZl NCDE IT, ZZZINU0E21 , ZZZl NUDE 3) ,ZZZINUDE4> , ZZZl CENTR)) 
GU TU 900 

202 CON I IN JE 


*** FUk contour PLOTS 

RPUINI = 1 
RING = -2 
1 31 ZL = 1 

I F ( NPluT . LO • 5 J CALL POINTS! R IND, 1 S IZE, ZZZ (XPT ) , ZZZ( YPT ) ) 

1 Fi lNFOR.Eo.ii CALL WNUOEI 
1 ZZZ l NUMPT 1 ,LlL (XPTl , ZZZ (YPT 1 .ZZZ.UPT 1 ,ZZZ(RADPT 1 , 

2ZZZ l NUM EL 1 , ZZZ l NCOE 1 ) ,ZZZ ( N0DE2J , ZZZ ( N0UE31 ,ZZZ(N0DE41 ,ZZZ(CENTR) J 
I F( 1NF0R.EU.21 CALL rtELRT ( 
iZZZi NOMPT i , ZZZ(XP T 1 ,ZZZ( YPT1 ,ZZZ( wPT) ,ZZZ(RAOPT) , 

2ZZZ INUMELl , ZlZ(NGDEI) , ZaZ ( NOO £2 1 , ZZZ(N00E31 , ZZZ (NUU E4 > , LIL ( C ENT R i 1 
CALL BGUN C ( 

IaZZ ( nUMPT 1 , ZZZ (XPT) , ZZZ (YPT ) ,ZZZlwP.T ) , ZZZIRADPT ) , 

ZZZzlNUMEL ) , ZZZ(NCDEll,iZZ(N00E2),ZZZ( NODE 3 J ,ZZZ(N0DE4) , ZZZ (CENTR) ) 
900 CONTINUE 
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Call ne.rame 

I Fl KPL0 T.cW.4J CALL CCRT2 
iloup = iloup+i 
lHNODc.tw.O) GO TO 150 
GO TO (700,600,500) KOD£ 
150 CONTINUE 

CALL PSTOP 
STOP 

END OP MAIN 
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SUBROUTINE DGCMNT 

0 

c *** this su tsRuuT i ne contains program documentation 
c 


c 

c 

C DESCRIPTION OF INPUT DATA CARDS 

C 
C 
C 

C TITLE CARD - CONTAINS ANY DESIRED ALPHANUMERIC INFORMATION IN COLS. 1-80. 


C NAMELIST OPTION - CCNTAINS VALUES TO ALLOCATE STORAGE IN BLANK COMMON LLL, 
C AND CONTROL VALUES NEEDED BY THE PROGRAM. 

C THt FOLLOWING VALUES ARE INCLUDED 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

o 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


NNDEST = ESTIMATE NUMBER OF GRID POINTS TO BE USED. VALUE MUST BE 
GREATER THAN UR EQUAL TO THE ACTUAL NUMBER OF GRID POINTS. 

** UcFAULT = 200 ** 

NELEST = ESTIMATED NUMBER OF ELEMENTS TO' BE USED. VALUE MUST BE 
GREATER THAN OR EQUAL TO THE ACTUAL NUMBER OF ELEMENTS. 

** DEFAULT = 20C ** 

KGEOM SPECIFIES SUBROUTINE AND CURRESPONDING METHOD OF INPUT FOR 
MuDEL GEOMETRY. 

KGEOM = 1 FOR GRID POINTS AND ELEMENTS READ FROM CARDS WITH USER 
SPECIFIED FORMAT. 

2 FOR NASTRAN OECK WITH CARD IDENTIFIERS LEFT ADJUSTED 
AND DATA IN COLUMN WIDTHS OF 8. 

= 9 FOR USER SUPPLIED SUBROUTINE - GE0M9 . 

** DEFAULT = 1 ** 

KDATA SPECIFIES SUBROUTINE AND CORRESPONDING METHOD OF INPUT FOR 
DATA T 0 BE REPRESENTED BY CONTOUR LINES. 

KDATA = 1 FOR SUBROUTINE DATA! TU READ IN DATA TO BE PLOTTED 
FROM CARDS WITH USER SPECIFIED FORMAT. 

■= 5 FOR SUBROUTINE OATAi TO READ IN DATA TO BE PLOTTED 
FROM TAPE20. 

= 9 FOR SUBROUTINE DATA*, A USER SUPPLIED SUBROUTINE. 

** DEFAULT = L ** 

NVALuS - USED IF KDATA - 5 10 SPECIFY THE NUMBER OF SETS OF 
DATA TO BE PLOTTED WHICH ARE READ FROM TAPE20. 

** DEFAULT = 0 ** 

IRESEQ = 0 FOR NG RESEQUENCING OF GRID' POINT NUMBERS. 

= 1 TO RESEQUENCE GRID POINT NUMBERS IN SAME ORDER 
AS THEY ARE INPUT. 

** DEFAULT = l ** 

KPLUT SPECIFIES THE TYPE OF OUTPUT DEVICE TU BE USED. 

KPLOT = 1 FOR CALCOMP. 

= 2 FOR CALCOMP WITH PLOTTING SPEED REDUCED TO USE LEROY PENS. 

= 3 FOR VAR1AN. 

•= 9 FOR CRT (USE CCC250 SCOPES AT LHC) 

** DEFAULT = l ** 

INFOR = I IF DATA TO BE PLOTTED IS SPECIFIED AT THE GRID POINTS. 

= 2 IF CAT A TU BE PLOTTED IS SPECIFIED AT THE ELEMENT CENTROIDS. 
* * DEFALLT = 1 ** 

XSPACE = SPACE BETWEEN PLOTS IN X-OIRECTION, IN INCHES. 

** DEFAULT = 10.0 ** 
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KSiGN = -1 ro CHANGE SIGN OF V ORDINATES. 

= 1 UG NCr change sign of y coordinates. 

** DEFAULT = 1 ** 

luCASt ■= U FOR NU TITLE CARD PRECEDING 
DECK CF DATA TO BE PLOTTED. 

= I FOR TITLE CARL PRECEDING 

DECK OF DATA TO BE PLOTTED. 

DEFAULT - 0 <■* 


MODEL GEJMcTR Y IS NOw INPUT IN ONE OF THE FOLLOWING FORMS. 

DEPENDING ON TEE VALUE OF KGE'GM SPECIFIED IN NAMELIST OPTION. 


C USE IF KGEOM — 1 


C (A i A SINGLE CARD CONTAINING THE WORO FORMAT IN COLUMNS 1-6 AND 

C A VARIABLE FORMAT CORRE SPOND I NG TO THE FORMAT OF THE GRID 

C POINT CARDS WITH LEFT PARENTHESIS STARTING IN COLUMN 11 

C AND UP TO SO COLUMNS MAY BE USED. 

C Ibi DECK OF GRID POINT CARDS. EACH CARD CONTAINS 3 VALUES, GRID 

C POINT NUMBER (INTEGER), X-CGORDINATE (REAL), AND 

C Y— COORDINATE (REAL). THE FORMAT IS SPECIFIED BY (A) ABOVE. 

C (Ci A SINGLE CARD CONTAINING THE WORD ENDGRID IN COLUMNS 1-7. 

C ID) A SINGLE CARD CONTAINING THE WORD FURMAT IN COLUMNS 1-6 AND 

C A VARIABLE FORMAT CORRESPONDING TO THE FORMAT OF THE 

C ELEMENT CARDS WITH LEFT PARENTHESIS STARTING IN COLUMN 11 

C AND UP TO BO COLUMNS MAY BE USED. 

C It) DECK OF ELEMENT CARDS. EACH CARD CONTAINS 5 INTEGER FIELDS 

C WHICH ARE THE ELEMENT NUMBER, AND GRID POINT NUMBERS AT THE 

o VERTICES CF ThE ELEMENTS. FOR TRIANGULAR ELEMENTS THE 

L LAST INTEGER FIELO MUST BE BLANK OR ZERO. 

C THE FORMAT IS SPECIFIED IN (D) ABOVE. 

C (F) A SINGLE CARD CONTAINING THE WORD ENDGEOM IN COLUMNS 1-7. 


C USE IF KGEOM = 2 

C (A) A SINGLE CARU CONTAINING JHE WORD TR I AEL IN COLUMNS 1-6 AND 

C UP TU NINE NASTRAN TRIANGULAR ELEMENT CONNECTION NAMES, 

C WHICH ARE LEFT- ADJUSTED IN FIELD WIDTHS OF B, STARTING IN 

C COLUMN 9 (COLS. 9-16, 17-24, ...» 73-80). THIS CARD CAN 

C BE OMITTED IF TRIANGULAR ELEMENTS ARE NOT USED FOR THE PLOT 

C Id) A SINGLE CARD CONTAINING THE WORO QUADEL IN COLUMNS 1-6 AND 

C UP TO NINE NASTRAN QUADRILATERAL CONNECTION NAMES, WHICH 

C ARE LEFT-ADJUSTED IN FIELD WIDTHS OF ’8, STARTING IN COL. 9 

C (COLS. 9-16, 17-24, ...» 73-BO). THIS CARD CAN BE OMITTED 

C IF QUADRILATERAL ELEMENTS ARE NOT USED FOR THE PLOT. 

C (C) A NaSTRAN BULK DATA CECK. ONLY THE GRID CAROS AND THE 

C ELEMENT CLNNECT ICN CARDS WITH NAMES MATCHING THOSE GIVEN ON 

G THE TRIALL ANO QUADEL CARDS WILL BE USED FOR THE PLOT. ALL 

0 CThER CARDS IN THE NASTRAN BULK DATA DECK WILL BE IGNORED. 

C (D) A SINGLE CARD CONTAINING THE wORD ENDGEOM IN COLUMNS 1-7. 
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C USE IE XGEOM = 9 

CALL SUBROUTINE GEUM9 WH ICH IS PREPARED BY THE USER TO READ 
GEOMETRY DATA. 


CASE IDt.MTIFICAT ION CARD. 

THIS CARD IS CMITTED IE IDCASE=Q IS SPECIFIED IN $QPTIQN. 

IF PRESENT, THIS CARD CONTAINS ANY DESIRED ALPHANUMERIC 
INFURMATLGN IN CCLS.1-80. RILL APPEAR BEFORE EACH DATA PLOT. 


DATA TO qE PLOTTED IS NOw INPUT IN ONE UF THE FOLLOWING FORMS, 

DEPENulNG ON THE VALUE OF KDAT A SPECIFIED IN NAMELIST OPTION. 


USE IF KDAT A = I 


(A) A SINGLE CARD CONTAINING THE WORD FORMAT IN COLUMNS 1-6 AND 
A VARUBLE FORMAT FUR THE DATA CARDS WITH LEFT PARENTHESIS 
STARTING IN COLUMN 11 AND UP TO 80 COLUMNS MAY BE USED. IF 
CONTRUL VARIABLES ARE INCLUDED FOR MORE THAN ONE GRIO 
POINT OR ELEMENT PER CARD, THE NUMBER OF GRID POINTS OR 
ELEMENTS PER CARO MUST BE ENTERED AS AN INTEGER IN COL. 8. 

IB) DECK OF DATA TO EE PLOTTED!. THERE CAN BE MULTIPLE DATA 

VALUE SETS PER CARD OR THE SET CAN EXTEND TO MORE THAN ONE 
CARO (UFTEN THE Case WITH NASTRAN PUNCHED OUTPUT) WHICH CAN 
BE HANDLED WITH A FORMAT FOR READING MULTIPLE CARDS. 

(Ci BLANK CARD UR CARDS TO END, DATA DECK. THE NUMBER OF BLANK 
CARDS MUST CORRESPOND TO THE NUMBER UF CARDS READ AT ONE 
TIME BY ThE SPECIFIED VARIABLE FORMAT. 


USE IF KDAT A = 5 

READS NV ALU S (FROM NAMELIST UPTIUN) SETS QF CONTROL VARIABLE 
VALUES FROM TAPE20. EACH SET OF CONTROL VARIABLES MUST HAVE 
BEEN WRITTEN C N TAPE2G AS AN UNFORMATTED RECORD. 


USE IF KDAT A = 9 

CALL SUBROUTINE CATA9 WHICH l Si PREPARED BY THE USER TO READ 
DATA TO BE PLOTTED. 


C NAMELIST PICT - CONTAINS VALUES NEEDED TU GENERATE PLOTS. 
C THE FOLLOWING VALUES ARE INCLUDED 
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NPLOT SPECIFIES TYPE OF PLOT JO BE GENERATED. 

NPCOT = 1 IF ELEMENT LAYCUT WITHOUT LABELS IS DESIRED. 

-2 IF ELEMENT LAYOUT WITH GRID PT. LABELS IS DESIRED. 

= 3 IF ELEMENT LAYOUT WJTH ELEMENT LABELS IS DESIRED. 

■=4 IF CONTOUR PLOTS ARE DESIRED WITHOUT SYMBOLS AT GRID POINTS. 
= 5 IF CCNTCUR PLOTS ARE DESIRED WITH SYMBOLS AT GRID POINTS. 

** DEFAULT = 4 ** 

XORGN = X-LGCATIGN OF ORIGIN OF FIRST PLOT, IN INCHES. 

♦A DEFAULT = 0.0 ** 

YORGN = Y-LOCATICN CF ORIGIN OF FIRST PLOT, IN INCHES. 

** DEFAULT = O.C ** 

PSCALE = MODEL SIZE REDUCTION FACTOR, PSCALE = ACTUAL MODEL 
SIZE/CESIREU PLOT SIZE. 

** DEFAULT = 1.0 ** 

ISCALE = METHOD OF SCALING CONTROL VARIABLE DATA TO BE PLOTTED. 

CUNTCUR LINES HAVE ONLY INTEGER VALUES ANNOTATED ON 
THE PLOT AND, THUS, THE DATA MUST BE SCALED SUCH THAT 
THESE INTEGERS WILL CONTAIN THE DESIRED NUMBER OF 
SIGNIFICANT DIGITS. THE DEFINITION OF WMAGS AND 
ICNTRS DEPENDS CN THE VALUE OF ISCALE. 

ISCALE = l, FOR USER SPECIFICATION OF SCALE FACTORS. 

= 2, FOR PROGRAM CALCULATION UF SCALE FACTORS TO GIVE 
THE USER SPECIFIED NUMBER OF SIGNIFICANT OIGITS IN 
ANNOTATION OF THE MAXIMUM ABSOLUTE CONTOUR LINE. 

= 3, FOR PROGRAM CALCULATION OF SCALE FACTORS TO GIVE 
WMAGS AS THE MAXIMUM VALUE OF DATA. 

** DEFAULT = 3 ** 

WMAGS - (FUR ISCALE = li , MAGNIFICATION OF CONTROL VARIABLES 
FOR ANNOTATION OF CONTOUR LINES ON PLOT. 

(FUR ISCALE = 2), NUMBER OF SIGNIFICANT DIGITS IN 
ANNOTATION OF MAXIMUM ABSOLUTE CONTOUR LINE. 

(WMAGS = i.O, 2.0, 3.0, ETC.) 

(FOR ISCALE = 3), MAXIMUM VALUE OF SCALED DATA, 

WMAGS MUSI BE LESS THAN OR EQUAL TO 2 (MAXIMUM CCNTQUR 
LINE IS INTEGER TRUNCATION UF WMAGS). 

** DEFAULT = 100. .0 ** 

ICNTRS - (FUR ISCALE =1), USER SPECIFIED CONTOUR INTERVAL 
(DIFFERENCE IN INTEGER VALUES OF ADJACENT CONTOUR 
LINES). 

(FOR ISCALE ■= 2 OR 3), APPROXIMATE NUMBER OF DIFFERENT 
CONTOUR LINE VALUES, THE CONTOUR INTERVAL IS 
CALCULATED BY THE PROGRAM. 

** DEFAULT = 1.0 ** 
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XLHT = HEIGHT OF INTEGERS TO BE ANNOTATED ON PLOTS, IN INCHES. 
** DEFAULT = O.Id ** 

NXLAB = TUTAL NUMBER OF LINES PARALLEL TO Y-AXIS ALONG WHICH 

CONTOUR LINtS ARE TO BE LABELED, MUST BE ,LE. 10, (ALL 
CONTOUR LINES ARE LABELED WHERE THEY INTERSECT. WITH 
THE Sc SELECTED LINES). THESE LABELS ARE IN ADDITION 
TO THOSE AUTOMATICALLY PROVIDED AT BOUNDARIES OF THE 
CUNTCUR SURFACE. 

** DEFAULT = 0 ** 

XLAd = ARRAY OF DISTANCES IN X-DIRECT1.0N FROM THE ORIGIN TO 

LINES PARALLEL TO THE .Y-AXIS ALUNG WHICH CONTOUR LINES 
ARE LABELED. THERE MUST BE NXLAB OF THESE VALUES AND 
THEY MUST BE IN UNITS UF THE URIGINAL (UNSCALED) MODEL. 
** DEFAULT = ALL ZEROS ** 

NYLAB - SAME AS NXLAB FOR LABEL LOCATIONS PARALLEL TO X-AXIS. 
** DEFAULT = 0 ** 
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YLAB = SAME AS XL AB FOR LABEL LOCATIONS PARALLEL TO X-AXIS. 

** DEFAULT = ALL ZEROS 

KOuE - SPECIFIES CONTROL OPTION AFTER PLOT IS COMPLETE. 

= 0, FOR LAST PLOT, EXIT FROM PROGRAM. 

= I, K E AC ANOTHER NAMELIST PICT. 

= 2, READ A Ntw SET OF CONTROL VARIABLE VALUES TO BE 

PLOTTED, INCLUDING A CASE IDENTIFICATION CARD IF PRESENT 
= 3, READ A COMPLETE NEW SET. OF INPUT DATA, INCLUDING A 
TITLE CARD. 

** DEFAULT = 0 ** 
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THc ABOVE COMPRISES A COMPLETE BASIC SET UF INPUT DATA IF 
RODE = 0 IN S.PICT. FOR KCOE = 1, 2, OR 3, ADDITIONAL SECTIONS OF 
THE BASIC DECK MUST BE REPEATED. THE DECK MUST END WITH 
NAMELIST SPICT HAVING KDDE = 0. 




DESCRIPTION CF GRAPHICS SUBROUTINES 


GRAPHICS 
SUBROUT IN E 

USEO BY FOLLOWING 
PLOTTING DEVICE 

CALCuMP 

CALCLMP 

LEROY 

CALCCMP 

PSEuUU 

VAR I AN 

CUC23 0 

CRT 

CAL PL T 

CALCCMP, VARI AN, CRT 

PNTPLT 

CALCOMP ,VARIAN,CKT 

NOTATE 

CALCCMP, VARIAN, CRT 

number 

CALCCMP, VARIAN, CRT 

nframe 

VARIAN 

NEXT 

CRT 

mESAGE 

CRT 

PARAMS 

CRT 

KFURMAT 

CRT 


CALLED BY FOLLOWING 
PROGRAM SUBROUTINES 

1NITAL 

INITAL 

INITAL 

INITAL 

MAIN, CCRT3,PST0P,TRIN, BOUND, LAY OUT 

POINTS 

MAIN 

TRIN, BOUND, LAYOUT 
MAIN, LAYOUT 

MAIN*CCRT I, CCRT 2,CCRT3 
MAIN,CCRT1,CCRT2,CCRT3 
C CRT I 
CCRT1 
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PURPGSc- 

USE 

COMMENT S 

PURPOSE 

USE . 
COMMENT S 

PURPOSE 

USt 

PURPOSE 

USE 


SUBROUTINE CALCUMP 

IMS IS THE NORMAL MODE PROCESSOR. THE NECESSARY 
PARAMETERS ANO LINKAGE ARE StT UP TO OUTPUT A TAPE 
FOR THE CALCUMP 78D/7fa3 0.010/0.005-INCH STEP PLOTTER. 

CALL CALC CMP 

THIS CALL MUST BE GIVEN BEFORE THE FIRST CALL TO A 
PLOTTING ROUTINE. 

SUBROUTINE LEROY 

THE PARAMETERS NECESSARY TO ACCOMODATE PLOTTING WITH THE 
LICUIO INK PEN ARE SET UP BY CALL LEROY. 

CALL LEROY 

THIS CALL SHOULD ONLY BE USED WITH THE CALCOMP PROCESSOR. 
IN ADDITION TC REDUCING THE SPEED OF THE PLOTTER FOR 
ALL PLOTTING MOVEMENTS* THE NUMBER OF PLOT VECTORS IN ANY 
ANNOTATION IS CONSIDERABLY INCREASED. 

SUDRCCTINc PSEUDD 

INITIALIZES PLOT VECTOR FILE FOR VAR i AN PLOTTER 
CALL PSEUDO 

SUBROUTINE CDC250 

INITIALIZES CATHODE RAY TUBE CONSOLE. 

CALL C DC 2 50 
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SUBROUTINE CALPLT 


TO WOVE ThE PLOTTER. PEN TO A NEW LOCATION WITH PEN UP 
OF CORN AND TO SIGNAL THE END OF A JOB SEGMENT BY 
INCREMENTING THE BLOCK ADDRESS NUMBER. 

CALL CALPLT (X»Y» IPEN ) 

WHERE 

X f V ARE THE FLOATING POINT VALUES FOR PEN MOVEMENT. 

IPEN = 2 PEN OOW.N 
= 3 PEN OP 

NEGATIVE OPEN WILL ASSIGN X=D, Y^O 
AS ThE LOCATION OF THE PEN AFTER MOVING THE 
X » Y (CREATE A NEW REFERENCE POINT! AND 
INCREASE THE BLOCK NUMBER BY ONE. 

ALL X AND Y COORDINATES MUST BE EXPRESSED AS FLOATING 
POINT INChES (ACTUAL PAGE DIMENSIONS) IN DEFLECTION FROM 
THE ORIGIN. 


SUBROUTINE PNTPLT 


TO CRAW NASA STANDARD PLOT SYMBOLS CENTERED ON A GIVEN 
COORDINATE VALUE. 

CALL PNTPLTl A,B»NU, IS) 

WHERE 

A IS THE X COORDINATE FOR THE CENTERED SYMBOL 

IN FLOATING POINT INCHES. 

B IS THE Y COORDINATE FUR THE CENTERED SYMBOL 

IN FLOATING POINT iNChtS. 

no is an integer specifying the symbol to be usedv 

IS IS AN INTEGER VALUE SPECIFYING THE SIZE SYMBOL 

TO BE USED. 


79 



C' r n c> 


APPENDIX C - Continued 


C 

L 

C 

c 

C PURPOSE 

c 

c 

c 

c ust 

c 

c 

c , . 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c , 

C 

c 

c 


SUBROUTINE NOTATE 


TO CRAW ALPHANUMERIC INFORMATION FOR ANNOTATION AND LABELIN 
ANU PRUV1UE SPECIAL CENTERED SYMBOLS FOR ANNOTATION OF 
GAIA POINTS. 

CALL NOTATE! XtY, HE I.GHT ,BCD,THtTA,Ni 
WHERE 

X , Y ARE THE FLOATING POINT PAGE COORDINATES OF THE 

FIRST CHARACTER. FOR ALPHANUMERIC CHARACTERS 
THE COORDINATES OF THE LOWER LEFT-HAND CORNER 
CF THE CHARACTERS ARE SPECIFIED. 


HEIGhT 


SPECIFIES CHARACTER SUE AND 
POINT INCHES FUR A FULL-SIZE 
WIDTH OF A CHARACTER WILL BE 
.SPACE BETWEEN CHARACTERS IS 


SPACING IN FLOATING 
CHARACTER. THE 
(4/71 *HE IGH T AND THE 
2/7i*HEIGHT. 


BCD IS THE STRLNG OF ALPHANUMERIC CHARACTERS TO BE 

DRAWN. 

THETA IS THE ANGLE IN FLOATING' POINT OEGREES AT WHICH 
THE INFORMATION IS TO BE DRAWN. 

N IS THE NUMBER OF CHARACTERS, INCLUDING BLANKS, IN 

THE LABEL* 
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comment s 


PURPOSE 

USE 


PURPUSE 


USE 


SUBROUTINE NUMBER 


TO CONVERT A FLOATING NUMBER TO BCD (EXPRESSED IN 
F FORMAT I , AND DRAW THE RESULTING APLHANUMERIC CHARACTERS. 

CALL NUMBER! X, Y, SI ZE.FPN, THETA.NI 


WHERE 

X,Y ARE THE COORDINATES IN FLOATING POINT INCHES 

CF THE LEFT LOWER CORNER OF THE FIRST DIGIT 
OF OUTPUT- 

SIZE IS THE HEIGHT OF THE PLOTTED NUMBER IN FLOATING 

POINT INCHES. 

FFN IS THE FLOATING POINT NUMBER TO BE DRAWN. 

THETA IS THE ANGLE IN FLOATING PUINT DEGREES AT WHICH 

THE NUMBER IS TO BE DRAwN. 

N IS THE NUMBER OF DECIMAL DIGITS TO THE RIGHT 

OF THE DECIMAL POINT FOR OUTPUT. 

N = -1 AND N = 0 BOTH SPECIFY NO DECIMAL PLACES, 
HOWEVERt -1 SUPPRESSES THE DECIMAL POINT. 

ThE NUMBER IS RESTRICTED TO A MAXIMUM OF 12 DIGITS. 

THE ROUTINE TRUNCATES THE FLOATING POINT NUMBER AT THE 
REQUIRED CECIMAL PLACE. 


SUBROUTINE NFRAME 


USED BY VARIAN PLOTTER TO ADVANCE PLOTTING FRAME. 

CALL NFRAME 

SUBROUTINE NEXT 

PROVIDES A BREAK POINT OR HALT CURING APPLICATION PROGRAM 
EXECUTION. OPERATOR MUST PRESS FUNCTION KEY TO RESUME, 
AND NUMBER OF KEY IS RETURNED IN CALLING PARAMETER- 

GALL NEXTINKEYJ 

WHERE 

NKEY IS NUM8EK OF FUNCTION KEY PRESSED. 
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APPENDIX C - Continued 

SUBROUTINE MESAGE 

PURPOSE PROVIDES THE CAPABILITY TO DISPLAY A MESSAGE ON THE CRT250. 

USE CALL MESAGE I I ,BCD,N> 

WHERE 

l INDICATES INTENSITY OF CHARACTER DISPLAY. 

BCD IS ADDRESS OF ARRAY CONTAINING THE MESSAGE 
IN HOLLERITH FORM. 

N IS THE NUMBER OF CHARACTERS IN THE MESSAGE ILESS THAN 50) 
SUBROUTINE PARAMS 

PURPOSE USED TO GENERATE A TABLE OF SYMBOLIC NAMES THAT CAN BE 

ACCESSED USING IHE ALPHANUMERIC KEYBOARD ON THE CRJ250. 

USE CALL P ARAMS ( BCD f VAR t 

WHERE 

BCC IS THE HOLLERITH REPRESENTATION OF SYMBOLIC NAME. 

VAR IS PROGRAM VARIABLE REFERRED TO BE SYMBOLIC NAME. 

COMMENTS UP TO i PAIRS CF VARIABLES MAY BE SPECIFIED IN A SINGLE 

CALL TU) PARAMS. TABLE HAS CAPACITY FOR 42 PAIRS. 

SUBROUTINE KFURMAT 

ALLOWS PROGRAMMER TO CHANGE FORMAT FOR KEYBOARD INPUT. 

CALL KFORMAT (NHBCDi 
WHERE 

N IS THE NUMBER OF CHARACTERS IN BCD. 

H IS REQUIRED. 

BCC IS THE RE CD I RED FORMAT (1 4 , F4. 2 , A 1 0 , ETC. ) 

RET URN 

END UF OOCMNT 


PURPOSc 

USE 
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SUBROUTINE CCRT1 

*** FUR CHANGING VALUES INPUT BY $PICT USING CRT. 

COMMON/CGNTRL/ 1 NFOR ,KGECM,KGATA* KSIGN.KPLQT , XORGN, YORGN, 

IP SC ALt, I SC ALE, WM AG, ICNTR , XLhT ,NXLAB , XLABl 10) , NYLAB , YLAB 11 0) 
CGMMUN/ SAVEV/ WMAGS, ICNTRS 
COMMON/ GK ODE/ KCCE.NPLOT 

COMMON/ SI 2EXY/ XM IMM , XMAX I M, YMI NiM , Y MAXI M, XSPACE , Y SPACE 

*** FOR INTEGER CONTROL VALUES 

CALL KFCKMAT (3HI Id 
CALL PARAMS 

CALL PARAMS(BLNPLCT,NPLCT,6LISCALE, I SC AL E , 6L I CNTRS , ICNTRS) 

CALL PARAMb I5LNXLAB, NXLA8.5LN YLAB , NYLAB i 

CALL MESAGEI 1 ,32PT0 CHANGE INTEGER CONTROL VALUES, 32 I 

CALL MESAGEI i , 13HVARIABLES ARttlii 

CALL MESAOEU .21PNPLQT , 1SCALE, ICNTRS, 21i 

CALL McSAGE 1 1 , 12HNXLAB , NYLAB, 12) 

CALL Me SAGE 1 1 , 17PANY KEY COM I NOES , 1 7 ) 

CALL NEXT(NKtY) 

*** PUR FLOAT ING POINT CONTROL VALUES 

CALL kFORMAT I5HF1G.3) 

CALL PARAMS 

CALL P ARAMS l 5LXURGN, XURGN, BL YORGN, YORGN, faLP SC ALE, PSCALEI 
CALL PARAMS (BLhMAGS , WMAGS ,4LXLHT,XLHT ) 

CALL PARAMS I 3L XL Adi , XL AS ID , 5L YLAB L , YLAB ( 1) ) 

LALL MESAGEI 1 , 39hT0 CHANGE FLOATING POINT CONTROL VALUES, 39) 
CALL MESauEI 1, 13FVARIABLES ARE, 13) 

CALL MESAGEI I .ZOHXORGN, YORGN, PSCALE,20) 

CALL MESAGEI I, ilPwMAGS, XLHT , i II 
CALL MESAGEI i» I2PXLAB1, YLAB1,12J 
CALL MESAGEII, 17PANY KEY CUNT INUES, L7 ) 

CALL NeXriNKtY) 

RETURN 

cND 0 F CCk T l 


SUBROUTINE OORT 2 

*** FUR SELECTING CCNTROL OPTIUN, KQUE, AT END UF JOB USING CRT. 
CUMMUN/ GKUut/ KGOE,NPLOT 


CALL MeSAGEIi, 

1 b h F N 

KEY 

34 

ENOS JOB 

,18) 



CALL MtSAGEIi, 

32HFN 

KEY 

3S 

TO 

ALTER 

EXIST 

ING 

PLOT, 32) 

CALL MtSAGEIi* 

4 2 h F N 

KEY 

36 

TO 

READ 

NEW SE 

T OF 

DISPLACEMENTS, 42) 

CALL MuSaGEI I , 
CCNT I NUE 

3 7HFN 

KEY 

37 

TO 

READ 

A COMP 

LETE 

NEW CASE, 37) 

CALL NcXTINKEY 

) 








IFINKEY.LT.34. 

LR.NKtY.G 

T.37J 

GO TO 

10 



1 F l NKtY . EQ.34) 

KODE 

= 0 







1FINKEY.EU.3B) 

KCCE 

= 1 







IFiNKl7.CU.36) 

KCCt 

= 2 







1FINKEY.EQ.37) 

KCOE 

= 3 








RETURN 

END UF CCRT2 
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SUBROUTINE CCRT3 

*** REMINDER TO PUT EOF ON PLOT FILE WHEN USING CRT. 

CALL ME SAGE (It 1H ,IJ 

CALL MESAGEU, 45FLAST REMINOER TO PUT EOF ON PLOT FILE, IF ANY, 45) 
CALL MESAGEU, 1H ,1) 

CALL MESAGEU, 4GF.C0 IT AT NEXT PLOT FILE COMPLETE MESSAGE, 40) 

CALL MESAGEU, IH ,1) 

CALL MESAGEU, 17FANY KEY CUNT INUES, 17) 

CALL NEXTINKEY) 

CALL CALPLT 112.0,0.0,-3) 

R t T URN 

END OF CCRT3 


SUBROUTINE PSTCP 
*** TO TERMINATE JOB 

CUMMuN/CUNTRL/ I NFOR »KGECM, KUATA, KSIGN.KPLUT , XORGN, YORGN, 
1PSCALE, I SCALE, WMAG, ICNTR , XLHT ,NXL AB ,XLAb(10) ,NYLAB, YLABI10) 
CALL CALPLT (0. 0,0.0,999) 

IF! KPL0T.EW.4) CALL CCRT3 
STOP 

end of pstop 
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SUBROUTINE INITAL 

*** lu SE1 VALUES FCR CONTROL PARAMETERS 

COMMON/ KU UN T / NiNCuE , NELMT .NNDEST, NELE ST 

COMMON/ CGNTRL/ INFGR ,KGECM, KCATA, KS IGN , KPLUT , XORGN, YGRGN, 
iPSCALE, l SCALE , KM AG , 1 CN TK , XLH T , NXL AB , XL ALi 1 10 I , NYL AB , YL AB ( 10 J 
COMMON/ GKUUE/ KCDE.NPLOT 
CuMMUn/CASEID/ I CCA SE 
COMMON/ SEONCE/ IRESEU 
COMMON/ VALUES/ NVALUS 
CGMMON/SAVEV/ WMAGS , ICNTRS 

COMMON/ S U EX Y/ XM INI M, XM AXI K , YM1N IM , YMAXI M, X SPACE , Y SPACE 
NAMELIST/ OPT ION/ NNDES T ,NE Lt ST ,KGEOM, KD AT A , NV ALUS , I RESEU, KPLGT , 
i 1 NFUR .XSPACE.KSIGN, l DC ASt 

*** UtSCRlPT ION OF VALUES IN ^OPTION GIVEN IN SUBROUTINE OGCMNT 


*** TO SET DEFAULT VALUES FGR iCPTIUN 

NNDES T - 200 
NELEST = 200 
KuECM ■= I 
kDATA - i 
NVALUS - 0 
IRESEw = 1 
KPLUT = 1 
INFOR = i 
XSPaCL = 10.0 
KSIGN = i 
ICCASE = 0 

*** TO SET DEFAULT VALUES FCR $PICT 

NPLUT = 4 
XORGN = 0.0 
YORGN = 0.0 
P SC ALE = 1.0 
ISCALE = 3 
to MAGS = 100.0 
ICNTRS = 10 

XLH T = 0.15 
NXL Ad = 0 

NYLAd - 0 
DO 10 1-1,10 
XL Ad ( I i = 0.0 
YL AtJI I ) = 0.0 
10 CONTINUE 
RODE = 0 
READI5 .OPTION) 

IFIKPL0T.LE.2l CALL CALCCMP 
IFIKPLUT . EQ .2 ) CALL LERCY 
IFIKPLOT .EG. 3) CALL PSEUDO 
IFIKPLUT. EG. 4) CALL CDC250 
wRl TEI6, OPTION) 

RETURN 

END OF INITAL 
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SOBRGU f INE GE0M1 (NUMPT, XPT , YPT ,WP T, RAOPT ,NUMEL , NODE 1 , N0DE2 , 

IN0DE3 ,NC0E4 ,CEMR) 

*** lu READ GEOMETRY OAT A FROM CARDS HAVING A GENERAL FORMAT 
C GMMGN/ KUUNT / NNLUE , NELKT « NNDEST , NELEST 

COMMON/ CGNTRL/ INFOR ,KGECM , K GAT At K.SI GN , KPLOT » XORGN , YORGN, 

IPSCaLE , ISCaLE, WMAG , I CNT R , XL hT , NXL A6 , XL AB < LO ) ,NYLAB, YLA8( 10) 
DIMENSION NOMPT! 1) ,XPT ( 1) ,YPT( 1) ,WPT! 1) ,RA0PT(1) , NOMEL I 1 J , 

INODE U 1 ) , NODE 2 ( 1 ) » NODE 3 ( l ) , NODE A! 1 J »CENTR( l i 
DIMENSION ABCDI8) 

DIMENSION FORM! t 7 ) 

TEST1 = 1GHENDGKID 
TEST2 = LOHENOGECM 
TEST3 = IGHFORMAT 

*** TO READ GRID INFORMATION 

100 CONTINUE 

HE AD I 5 » 10 ) ABC D 
10 FORMAT ( 8A 10) 

IF(AbCOill.EO.TESTi) GO TO 1000 
IFIABCU(L).BQ.IEST3T SO i 60 
30 CONTINUE 

DECODEIaOf55»ABCD) FORM T 
33 FORMAT! IOX, 7A10) 

GO TO 100 
o 0 CONTINUE 

NNODE - NNODE *• 1 
I FINNODE.GT. NNDEST) 550,531 

550 CONTINUE 
wRlTE(t>,535) 

555 FORMAT! IX ,///, 1CX,*XXXXX SORRY, THE ESTIMATE OF NUMBER OF NODE 

1 POINTS WAS EXCEEDED XXXXX*) 

STuP 

551 CONTINUE 

DECODE! BOrFURMT, AECDJ NUMPT i NNODE J , XP T ! NNODE ) , YP T ( NNODE ) 

GU TO 100 
1000 CONTINUE 

*** TO READ ELEMENT CONNECTION INFORMATION 

200 CONTINUE 

READ ( 5 , 10 J ABC 0 

IF! AbCOiU.Ew. TEST2) GO TO 2000 
IF! AoCUl U. EU.TEST3J 250,260 
2aO CONTINUE 

uECUDt ! BO , 55 , ABCQ ) FORM T 
GO To 200 
260 CONTINUE 

nelmt = NELMT+I 
lFiNELMT.GT.NELESTi 650,651 
650 dONTINUE 

WRI TE 16 ,536) 

556 FORMAT! IX,///, 10X,*XXXXX SORRY, THE ESTIMATE OF NOMBER OF ELEM 

ItNTS WAS EXCEEDED XXXXX*) 

STOP 


86 



APPENDIX C — Continued 


651 CONTINUE 

DECODE 80. FOKMT ,ABCO) NUMEU NELMT) . NODE i(NELMT) .NUDE2 (NELMT ) , 
INODE3 (NELMT ), NOD E4 ( NELMT ) 

GO TU 200 
2000 CONTINUE 
KETURIM 

END OF GEUM1 
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SUBROUTINE GEGM2 4 (vUMPT » XPT ,YPI»WPT*RADPT » NUMEL » NODE l » NU0E2 » 
1 NODES , N0DE4 , CEN TR ) 

*** TU REAU NASTRAN GEOMETRY DATA 

COMMON/ KUUNT/ NNCCE , NELHT , NN0EST-, NELE ST 

COMMON/CONTRL/ I NFOR , KGEGM, KOAIA, KS I ON, KPLGT , XORGN, YORGN, 
1PSCALE, 1 SCALE, *MAG, 1CNT R , XLFiT » NXLA6 , XLAb l 10) , NYLAB, YLAB110) 
D l MEN SI UN NUMPFl 1) , XPT! 1) ,YPTI D.WPT I 1 ) » R ADPT ( l ) ,NUMEL4 1) , 
1NUUE1 ( 1 1 ,NQUE2( 1 J , NODE 3 ( I J , NOOE44 ii »CENTR( ii 
DIMENSION ABCO I 8 ) 

DIMENSION NDPT44) 

DIMENSION . TR] A EL 4 9 J ,0UADEL(9) 
iO FORM AT I 8 A 10 J 
TESTG = IOH 
TEST I = 10HGR ID 
TESTS = 10HTRIAEL 
TEST4 = 10HOUADEL 
TESTS = 1UHENDGECM 
DU 13 1=1,9 
TRI AELt I ) = IOH 
UUADEL ID - IOH 
IS CONTINUE 
3UQ CONTINUE 

READ! 5, ID) AeSCU 
IF!£uF,5) 1111,750 
750 CONTINUE 

UECUUE (60*50, ABCO) *ORD 1 
50 FORMAI I A8) 

IF { wjKii l.EQ.TESTO) GO TC 800 
I F( iaOkDl.Eu.TEST l ) GO TC 200 
I F ( WORD 1. EG . T E S I 3 J 62,63 

62 DECODE i 80,101, A3C0J I TR I ALL II J , 1= 1 , 9i 
101 FORMAT I 8X,9A8) 

GO lu 600 

63 CONTINUE 

I F ( NURD 1 . EW . TE S I 4 ) 64,65 
o4 DECODE I 80, 101 , A BCD) (UUAOELt 11 , 1=1,91 
GO TU 300 
65 CONTINUE 
DO 70 I =1,9 

IFIwURDl.EQ.TRIAELl III GO TC 300 

IF! NUR01.EW.QUAGEL4 U I GO TO 400 
70 CONTINUE 

IPIWOROl. EU« TESTS) GO TC 20G0 
GO TO 300 

*** TO READ GRID CARCS 


2u0 CONTINUE 

DtCUDcIdO ,201, A8CCJ I0UM , XDDM, YOUM 
201 FORMAT! 6X ,A8 , 8X , 2F8 . 01 
CALL IRiTE(iDUM) 

NNODE = NNGDE+1 
iF(NNOUE.GT.NNUtSI) 550,351 
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^50 CONTINUE 

WRITE (8, 355) 

553 E uRMAT ( 1X,///,1CX,*XAXXX SURRY, THE ESTIMATE UE NUMBER OF NUDE 

1 POlNTo WAS EXCEEDED XXXXX*) 

STOP 

35i CONTINUE 

NUMP T 1 NNUDE ) — ICCM 
XPT (NNUDE) = XDLiM 
YPT (NNUDE) = YDUM 
GU TU 800 

*** Tu READ CAROS CCNTA INING ELEMENTS WITH 3 GRID POINTS 

3 00 CONTINUE 

0ECUUc(80*301»ABCC) IUUM,NDPT( 1 ) ■» NDP T ( 2 ) , NOPTI 3) 

301 FOKi'lAT (aX,A8,8X,3A8) 

CALL I R I T E ( IDUM) 

CALL I R I T E ( NDPT ( I ) ) 

CALL 1R I TE( NDPT ( 2 ) ) 

CALL IR1TE(NDPT (3) i 

NllMT ^ NELMT+1 . ... 

lF(NELMT.GT.NELEST) 650,651 

050 CONTINUE 
WR1TE(6,556) 

5oo FURMAK IX ,///, 10X ,*XXXXX SORRY, THE ESTIMATE OF NUMBER OF ELEM 

IENTS WAS EXCEEDED XXXXX*) 

STOP 

051 CONTINUE 
NUMEL(NELMT) = IDEM 

NUDEi (NELMT ) = A OF T I 1) 

N00E2 (N ELMT ) ■= NDPT ( 2) 

NODEdi NELMI ) = NOP 1(3) 

N0UE4 (NELMT ) = C 
GO TO 8 00 

*** TO READ CARDS CONTAINING ELEMENTS WITH 4 GRID POINTS 

400 CONTINUE 

D EC ODE ( 80, 40 1, ABC C) IDUM,NDPT( I), NDPT (2) , NOP T ( 3) , NDPT (4) 

401 F0RMAT(8X,A6,8X,4A8) 

CALL I R 1 T E ( 10UM) 

CALL 1 R I T E ( NDPT ( 1 ) ) 

CALL i R I T E ( NDPT ( 2 ) ) 

CALL I K I T E ( NOP T ( 3 ) ) 

CALL I R IT E ( NDPT ( 4 ) ) 

NElMT = N E L M T 4 1 

I F(NELMT.GT.NELEST) 761,751 
761 CUNT INUE 

WRITElo»55o) 

STOP 

751 CONTINUE 

NUMEL ( NElMT ) - I DOM 
NOD El (NELMT) ^ NDPT 1 1) 

NU0E2(NELMT) = NDPT(2) 

N0DE3 (NELMT ) = NDPTI3) 

N00c4lNELMT ) - NLPT(4) 

Gu TU dOO 
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ti.ii continue 
2000 CONTINUE 
RETURN 

END uF 0ECM2 


SUbRUuT 1NE IRlTEINUMi 
C 

TO RIuH T ADJUST INTEGERS IN A F IEiO WIDTH OF EIGHT 

DIMENSION N ( 8 ) 

LANK. = 1H 

DECODE 1 a , 1 , NUM ) NSAVE 

1 FORMAT 1 131 

DECODE! a, 2 , NUM) (MI ), 1 = 1 , 8 ) 

2 FuRM AT ( 3A 1 ) 

00 10 1 = 1,8 

II = 9-1 

IF INI I I ) .NE.LANK) GO TO 20 
10 CONTINUE 

20 NUM = NSAVE/1 10**18-1 1) ) 

RET URN 

END OF IRITE 


SUBROUTINE GE0M9 I RUMP T , X PT , V P T , WP I, R A DP T , NUMEL , NODE 1 , NODE 2 , 
iNuDEi,N0Dt4»CENTR) 

*** USER SUFPLIED GECMETRY INPUT ROUTINE 

COMMON / ROUN 1 / NNUC E , NEL M T , NNCt ST , NELtST 

COMMON/ CUNT RL/ l NFOR ,KGECM , KCAT A, KS IGN.KPLOT , XURGN, YORGN, 

IP SCALE, IaCALE,WMAG, ICNTR,XLHT,NXLA3,XLABI10) ,NYLA8, YLABI10) 
DIMENSION NUMPTU) ,XPTll) ,YPTI litWPTl 1) , RADP T 1 1 1 , NUMEL I 1) , 
iNOOEil I ) , NUDE 2 1 I) , NODE 3 I i) , N0DE4I 1) ,CENTK( 1) 

*** INSERT ROUTINE F.ERt 

RETURN 

END uF GE0M9 
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SUBROUTINE ZERORUUMPT, WPT,NUMtL ,CENTRi 
C 

*** to ZERO our VALUES of data to be plotted. 

0 UMMUN/ KUUtsT / NNCCE.NELMT .NNDEST , NELEST 

uUMMUN/ CONTRL/ l NFOK .KGEOM , KCA TA. KS IGN .KPLGT , XORGN , YORGN, 

IP SCALE, IjCALE, WMAG.ICNTR, XL hT.NXLAB »XLAB( 10I.NYLAB, YL A8 ( 10 ) 
DIMENSION NUMPTUJ ,wPT( 1 1 , NUMEL U I , CE NIR ( II 
GO TO { LOG, ZOO) 1NFCR 
100 CONTINUE 

DO 1.10 di, NNDEST 
WPT I 1 I = 0.0 
110 CONTINUE 
GO TO 300 
ZOO CONTINUE 

DO 210 1=1 , NELEST 
CENTKI1J = 0.0 
210 CONTINUE 
500 CONTINUE 
RET URN 
END OF lERQw 
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SUBROUTINE DAT AI ! NUKPT f WPT .TvUMEL, CENT R) 

*** TU RE All THE DATA TO BE PLOTTED 

COMMON/ KOUNT/ NNCDE , NELRT .NNCEST , NELEST 

COMMON/ CONTRL/ I NFOR ,KGECM , KCATA, KS I GN , KPLQT , XORGN , YORGN, 

IP SC ALE, I SCALE, WMAG, I CNT R , XLHT , NXL AB , XL AB { 10 ) , NYL AB, YLAB ( l 0) 
COMMON/ SEONCE/ IRESEG 

DIMENSION NUMPTil) ,WPT! 1} ,NUMEL{1 ) , CENT R{ 1) 

DIMENSION ABCJIBJ 
DIMENSION PORMT ( 7 ) 

DIMENSION ISAV! lC).DSAV(lO) 

10 FORMAT ! dAlO) 

TEST - oHFURMAT 
RE AD t 5 , 10) ABCO 
DECUJtiaO,45,ABCD) WORD, KVALU 
•+5 FORMAT ( Au, 1X,1 l) 

IFtKVALU.EQ.O) KVALU = 1 
l F t WORD • 2 U. TE S T ) 300*200 
200 WR1TE(6*20) 

20 FORMAT! 1H1 ,/// ,2CX,* SORRY, FORMAT FOR DATA NOT GIVEN*) 

S TUP 

300 LON T I NOE 

DECODt ( BO ,50, ABCO) FURMT 
50 FORMAT! 10X.7A10) 

100 CONTINUE 

READ ! 5 , FGRMT ) !1SAV!K),QSAV!K) ,K=1, KVALU) 

GO TU 1510,610) INFOR 

*** FOR DATA AT GRID POINTS 

510 CON T 1 NU E 

DU 560 K=l, KVALU 
I DOM - 1SAVIK) 

IFl IDUM.Ew.O) GO TU 100C 
WDUM - CSAViK) 

IF! IreSEw.LU . Li GC TO 520 
WPTUDUM) - WDUM 
GO TO 560 
520 CONTINUE 

UU 500 J=1»NN0UE 

IF ( NUMPT! J) .EG. iCUM) 501,500 

;>01 WPT ( J ) = WDUM 
GO TO 550 
500 CONTINUE 
550 CONTINUE 
5o0 CONTINUE 
GO TO 100 

*** FUR uaTA AT ELEMENT CENTROIDS 

o 1 0 CUN T 1 NU E 

DO o60 K= 1 , KVALU 
I CUM = 1SAV1K) 

IF! lDUM.EO.O) GG TO 1000 
WDUM = US A V ! K ) 

IF i IRESEw.EQ. 1 i GO TO 620 
CEN IRI I GUM) = wDUM 
uG TU oou 
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APPENDIX C - Continued 


CuN I I NU E 

DO 600 J= 1 , NEL M T 

iFINUMcU Ji .EQ.IOUMJ 601,600 

CfcNTK(J) - WOUM 

GO 10 o 5 0 

0 G N T 1 NU E 

CONTINUE 

CuNTl NO E 

GO TO 100 

6 ON 1 1 NO E 

RcT URN 

END OF 0ATA1 


6 Ub ROUT INE DAI A5 I NUMPT , InPT , NUMELi CENT R) 

TO READ DISPLACEMENT DATA FROM TAPE20 
COMMON/ ROUNT/ NNLOE ,NELMT ,NNOEST, NELEST 

CuMMON/CONT RL/ IN FOR ,KGECM, KDATA, KS 1GN , RPLOT , XORGN , YORGN , 

IP SC ALE, I SC ALE, wMAG, I CNTR , XL FT , NXL AO , XL AB I 10 ) , NYL A 8 , YL A6 ( 1 0 ) 
COMMON/ VALDES/ NVALUS 
COMMON/SEGNCE/ IRESEQ 

DIMENSION NUMPT I 1) * WPT I 1J » NUMEL 1 1 ) , CENTRI 1) 

Do 10 I=i,NVALOS 
R E AD l 20 J 1UUM , ROOM 
GO TO I510»6l0i 1NF0R 

Fur oaia at uRId points 

CON II NOE 

1FI IRESEQ. EQ. 1 I GG TO 520 
wPTIIDuM) = ROUM 
GO TO 10 
CONTINUE 

DU 500 J= 1 , NNOOE 

I FI NUMP T( J J .Eq. IOUM) 501,500 

WPT I J J = ROOM 

GO TO 10 

CONTI NO E 

GO TO 10 

FOR DATA AT ELEMENT CENTROIDS 
CON (I NOE 

IFl IRESEQ. EQ.li CO TO 620 
wPllluUM) — ROOM 
Gu TO 10 
CONTINUE 

DO 600 J=1,NELMT 

IFINOMEU J) .Eu.lCUMJ e01,60C 

CENTRI J) •= WOUM 

GO TO 10 

CONTI NOE 

CONTINUE 

RETURN 

END OF CAT A5 
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SUtiRUUT INE OAT A S I NUKPT , |»PT , NUMEL, CENTRJ 
C 

*** USER SUPPLIED w CATA INPUT ROUTINE 

COMMON/ KOUNT/ Ms ODE ,NELMT ,NNDEST , NELEST 

COMMON/ CUNTRL/ IN FOR ,KGEGM , KDA TA i KS l C N , KPLOT , XORGN , YQRGNt 
IP SCALE, I SCALE, *MAG, ICi\iTR,XLhT , NXLAB ,XLAB ( IU ) , NYL A6, YLABI 10) 
DIMENSION NOMPT l 1) ,WPTt 1) , NOMEL U ) , CE NTRl 1) 

*** INSERT ROUTINE HERE 

RETURN 

END uf DA T A 9 
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SUB ROUT IinIE PNT OU.T! lOUT.NUMPT.XPT* YP.T, WPT.RAOPT.NUMEL, 

INODE I , N00E2 , NODE3,NCUE4»CENTR) 

*** PGR PkINTED OUTPUT OF INFORMATION IN BLANK COMMON - ZZZ 

COMMON/ KOUNT/ NNOOE ,NELMT .NNOEST, NELEST 

COMMON/ CONTRL/ INFUR.KGECM.KDATA-, KS l&N*KPLOT , XQRGN, YORGN, 

1P3CALE, I SCALE, WMAG, ICNTR.XLHT , NXLA6 ,XLAB,( 101 , NYLA6, YLAB I 10) 
DIMENSION NUMPT(l) , XPT { 1 ) , YP H 1 1 , WPT I I) , RADP T (1) , NUMEL I 1) , 
lNOUci! I ) .NODE 2 I i 1 « N0DE3 ( 1 1 .N0DE4I 1) .CENTR(i) 

GO TO ( LOCO . 2000 1 ICUT 
1000 CONTINUE 

*** FOR OUTPUT UF GEOMETRY INFORMATION 
WRITE16.16) 

16 FORMAT!///. 5X,*GRIU POINT INFORMATION*,///) 
rtRlTE(6»17) 

17 FORMAT! 5X,*RESEQUENCED*,4X,*USER INPUT* ,35X ,*BOUNDARY*/ 

15X , *GR I D POINT*, 5X,*GRI0 POINT*, 7X,*I NPUT * , 10X,*INPUT*. 

2 tt X , *PUi NT*/ 

36 X, *NUMBER* »9X ,* NUMBER* , 1.3X , *X* , 1 4X , * Y* , 1 OX , * IND l CA TOR*// ) 

00 30 I = 1, NNOOE 

WRITEIB.IB) I, NO MPT ( 1) ,XPT( 1 ) , YP T ! I) .RADPT! I) 

18 FURMATi2X,l 10 , 5X , 1 1 0 , 3X , 3E1 5 .4 ) 

30 CONTINUE 

WHITE! o,l 4) 

19 FORMAT !///, 5X , *ELEMENT INFORMATION - WITH RESEQUENCED GRID POINTS 
1 *///) 

WRITE!6,21) 

21 FURMAT ! 6X,*RcStQCENCcD*,4X,*0SER I NPUT*, 19X,* GRID POINTS*/ 
J.5X,*tLtMENT*,8X, *ELEMENT*/ 

2i>X, *NUM8ER*»9X , *NUMBER* , 13X ,* 1*, 9.X , *2* , 9X , * 3* , 9X , *4*//) 

DO 38 l - l.NELMT 

WRI lc!o,22) I . NUMEL ! I) , NODE 1 ( I ) ,NGDE2( I ) .NODE 3 1 I) , N0DE4 ( I ) 

22 FORMAT ! 2 X , I 10, 5 X » I10.4X.4I1C) 

35 CONTINUE 

RETURN 

2000 CONTINUE 

*** FOR OUTPUT OF DATA TU BE PLOTTED 


WKITE(o,210) 

210 FORMAT! / / / , pX , *0 ATA TU BE PLG I TED* , / / / J 
I F ! I NFOR. Ew • 1 ) WRITE !6» 220) 

220 FORMA T l 5X , *RESEQUENCED GRID POINT NUMBERS AND DATA VALUES*,//) 
IF! INFUR.EQ.2) WRITE!6,230) 

230 FORMAT! 5X,*RESEGUENCEU ELEMENT NUMBERS AND DATA VALUES*,//) 

1 F ! I iIFOR. EQ . 1 ) W P I T E 16 » 240 ) i ( I , M PT ! I ) ) , 1 = 1 , NNODE ) 

IF l INFUR. EQ.2) WRITE16,240) l ( l , C EN TR ! I ) ) , I = 1 , NEL MT ) 

240 FORMAT! 5I2X, 18 ,E15.5i) 

RETURN 

END OF PNTGUT 
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SU6R0UT IN t XYSCAL1NUMPT,XPT,YPT,WPT ,RAOPT ,NUMEL , NODEi »NODE2 . 
INODE 3 »NL>DE4 » CENT R ) 

*** TO DETERMINE MAXIMUM DIMENSIONS QF THE DATUM SURFACE 
AND RENUMBER ELEMENT GRID POINTS. 

*** w DATA IS SCALED IN SUBROUTINE SCALEW 

COMMON/ KOUNT/ NNOUE ,NELMT .NNOEST-, NELL S! 

COMMON/ CONTRL/ INFQR ,KGECM , KDATA, K.S I GN * KPLOT , XORGN, YORGN, 

IP SCALE, ISCALE,WMAG,ICNTR,XLHI,NXLAB,XLABI 10) ,NYLAB, YLAB(IO) 
COMMON/ SECNCE/ IRESEQ 

COMMUN/SI ZEXY/ XM IN I M , XMAXIM » YMIN I M , Y MAXI M, XSPACE ,Y SPACE 
DIMENSION NUMPTI 1 ) » XPT ( 1 ) » YR II 1 ) * WP T (1 ) ,RADPT 1 1 ) , NUMELI 1 ) , 
INODE H i) t NUDE 21 1) ,NG0E3( l) ,NGDE4Il) »CENTR( L ) 

DIMENSION N0PTI4) 

*** TO SCALE DATA ANC DETERMINE X MINIMUM AND MAXIMUM 

XKINIM = XPTII) 

XMAXIM = XPTI I) 

YMIN1M = YPTIU 
YMAXIM = YPTl 1) 

DO tlO 1 = 1, NNODE 

IFIXPTI D.LT.XMIMM) XMINIM = XPTII) 

IF l XPTI I) .GT. XMAXIM) XMAXIM = XPTII) 

IFIYPTi I) .LT.YMIMM) YMINIM = YPTII) 

IPl YPTI 1) .GT. YMAXIM) YMAXIM = YPTII) 
iiO CONTINUE 

IF! IRESEO.NE.l) GC TO 7C0 

*** TO RENUMBER ELEMENT NODES 

DO oGO 1= l » NELMT 
NOPTli) = NODE l l I) 

NDPTC2) = NODE 2 I I) 

NOPTI 3) = NODE 3 I I) 

NDP T I 4) = NUDE4I l) 

DO t>U5 J-1,4 

I F ( N OP T l J ) • EU . 0 ) GO TO 605 
DU 610 K= 1 » NNGDE 
IF INUMP TiiO.Eu.NCPTIJ) ) E15.6IC 
615 IFtJ.Ew.l) NODE III) = K 

IFIJ.EU.2) NUDE 2 I I ) = K 
IF ( j.Ew.3) NODE 3 1 I ) = K 
IP t J. EW.4) NODE 4 l I ) = K 
Gu TO 601 
610 CONTINUE 

wRITE<6,555) NDP Tt J I .NUMEL 1 1 ) 

555 FORMAT I 5X ,* NUDE *,110,* IN ELEMENT *,110,* NOT IN USER OATA*) 
S fuP 

oOl CONTINUE 
605 CONTINUE 
oOO CONTI NUfc 

/oo continue 
re r URN 

END UF XYSCAL 
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SU6R0UT INE SCAL£W{wPT*CENTR) 

*** THIS SUBROUTINE SCALES W DATA 

COMMUN/.KuUNl / NN CUE t NELKT »NNDEST»NELEST 

COMMON/ CUNT RL/ INFOR ,KG£GM, KDATAi KS I ON, KPLOT ,XGRGN, YGRGN, 
iP SCALE, 1 SCALE , WMAG , i CNTR t XLHT , NXLA6 , XLABl 10 > , NYLAB, YLAB i 10) 

COMMON/ SAVEV/ WM ACS , 1CNTRS 
DIMENSION wPT( i) ,CENTRl 1) 
oG TO ( 1U0 , ECU ) INF OR 

*** FOR NODAL DATA 

100 CONTINUE 

ini M I N = WPTlii 
WMAX = WP T ( 1 ) 

DO 120 I- i » NNUDE 

IFIWPT { I ) .LT.WM1N) WMI N ^ WPT(l) 

IFIwPTt I) .GT.WMAX) WMAX = wPT(l) 

L20 CONTINUE 

wBIG - Ad $ l WMI N) 

IFIABS(wMAX).GT.WBIG) WBIG = ABSlwMAX) 
wRANG = A6S I WMAX-WM IN ) 

WKITEto, I 2 1 i WMAX.WMIN 

121 FORMAT I IX ,// / , 5X » *MAX I MUM VALUE QF DATA = * , E20. 8 , 10X , *M INI MUM VALU 
IE OF DATA =* » E2C. 8 ) 

GO TO I ICO, 160, 1BC) 1SCALE 
140 CONTINUE 

WMAG = WM AGS 
ICNIR = ICNTkS 
GO TO 190 
ISO CONTINUE 

wMAG = WM AGS 

LEXP = AL0G10IWBIG) 

I F< LEXP.LT.O j LEXF = LEXP-i 
IFILEXP.EQ.01 GG TO 161 
WMAG ^ WMAGS/l 1G.G**LEXP) 

161 CONTINUE 

WRANG = WRANG*WMAG 

CINT = WRANG/FLGAU ICNTRS) 

00 TO BOO 
1 60 CONTINUE 

IFIwMAGS.LT. 2.0) 333 , 33 A 
333 WRl lEIo ,336) 

333 FOR MAI I IX ,///» 1QX»* SORRY, VALUE OF WMAG MUST BE 2.0 OR GREATER*) 
STOP 

334 CONTINUE 

WMAG - WMAGS/WBIG 

WRANG = WKANG*WMAG 

CINT = WKANG/FLCATi ICNTRS) 

BOO CALL INTRVL (CINT , ICNTR ) 

190 CONTINUE 

WKl IE lo, 122 ) WMAG , l CNT K 

122 FORMAT! IX,///, 5X,*WMAG =* ,E20. 8,10X»*-ICNTR =*,16) 

GO T 0 300 

*** FOR ELEMENT DATA 
200 CONTINUE 


97 



APPENDIX C - Continued 


CMIN = CtNTRIlJ 
C HA X = CENTRIl) 

OU 220 l~ l » NELMT 
IFtCtNTRCli .LT.CMIN) CMIN = 
IFICENTRl I J .GT .CMAXJ CMAX = 
C 6 1 G = AB S I CM 1 N J 
IF! ABSI CMAX) .GT.CBIGJ CBlU 



CHANG = 

ABS I CM AX— CM I 

NJ 

220 

CON T l NO 

E 



HHi TE I o 

,121) CMAX, CM 

IN 


GU TO ( 

240,260,2801 

ISC 

240 

CUN T I NO 

E 



wMaG = 

WMAGS 



ICNTH - 

ICNTRS 



uU TO 290 

260 CONTINUE 
rtMAG = wMAGS 

lEXR - ALOG 10( CB IG ) 
IFILEaP.LI .OJ LEXP = LEXP-1 
IFILcXP.EO.OJ GU TO 261 
WMAG = *MAGS/l 10.C**LEXPJ 

261 CONTINUE 

CHANG = CHANG*WMAG 
CINT = CKANG/FLCATI ICNTHSJ 
GG TO 900 
260 CONTINUE 

to MAG = toMAGS/CBIG 
CHANG = CHANG* WMAG 
CINT = CKANG/FLOATI ICNTHSJ 
900 CALL INTKVHC1NT, ICNTRJ 
290 CONTINUE 

WRi TE I 6 » 122 ) IN MAG , I CNT H 
300 CONTINUE 
K t T U H N 

END OF SCaLEW 


CENTRI I J 
CENTH l I J 

= ABSI CM AX J 
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SUGRuuT INfc INTRVLIC INT , ICNTR) 

*** TO CAoCUL AT E CONTOUR INTERVAL. ICNTR 

DO 165 1=1 , 10 
CUP = 1.5*1 10.0**I)/10.0 
IP l C1NT .LT.CUP) 166,167 
l6o ICNTR = 1*110**1)710 
GO TO 140 
167 CONTINUE 

CUP = 3.5*1 1G. 0**11/10.0 
IFICINT.LI. CUP) 166 , 169 
166 ICNTR = 2*1 10** I ) / 10 
GU TO 140 
164 CONTINUE 

CUP = 7.5*( 10.C**I J /10.0 
IFICINr.LT.CUPi 170,171 
I/O ICNTR = 3*1 10**1 J /10 
GO TO 140 
171 CONTINUE 
lo5 CuNTl NOE 

wRITE(6»175) 

175 FORNAX 1a,///,5X,*SURRY, CONTOUR INTERVAL NOT FOUND, CHECK INPUT*) 
STOP 

140 CONTINUE 
RETURN 

END OF INTRVAL 
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SU8KGU T INE ENCLOSINUMPT , XPT »YPT ,WPT ,R-A0PT-,NUHEL,NQDEl t NQDE2, 
1NU0 c3 , N0DE4 , CENTR J 

*** TO CALCULATE ENCLOSED ANGLES AROUND NODE POINTS 
COMMON/ KGUNT/ NNOCE , NELP T , NNOE ST , NE LE ST 

DIMENSION NUMPTI 1 J t XPT 1 1 ) , YPT l 1 J » WPT III* R ADPT (1) , NUMELI 1) , 
LNCOEil 1 1 , NUDE/ I I i , N0UE3 1 ii , N00E4I 1 J »C ENTR ( l J 
DIMENSION XXI 3) »YY(3)»wl*C3) ,RRADl,3) 

DO 25 I^i,NNDEST 
25 KAOPTUJ = 0.0 
DO 500 1-= I » NELM T 
NDi = N00E1I I I 
ND2 = N0QE2 III 
NO 3 = NOD £3 1 1 J 
ND4 = NOD £4 I I ) 

XX (1) = XPTINOll 
YYI1I = YPTIND1I 
XX ( 2 1 = XPTIND2J 
YYI2I = YPTINU2) 

XXI 3 ) = XPTIND3I 
YYI3I ^ YPTIND3I 
CAlL AROUND ( XX t Y Y • RR AU ) 

RADPT (NDI ) = R ACP T ( NDI I +RRAU 1 1 i 
R AuP T I ND2 J - RADPIIN02) +KKAUI2I 
RADPT IN 031 = RADPTIND3 J+RRADI3I 
I E ( Nu4. Eg .0 J GO TO 505 
X X ( 2 i = XPTIND4) 

YY(2i = YPTIND4J 

CALL AROUND I XX * Y Y »RRAD ) 

RADPTINUI) = RADPTINDD+RRADtil 
RADPT IN04J = RADPT! ND4i +RRADI2I 
RADPTIND3I = RADPT ( ND3I +RRADI 3i 
5 D 5 CONTINUE 
500 CONTINUE 

PI = 3.14159265 
EPSIl 5 . 0*P I / 1 80 . 0 

TEST ■= 2. 0*PI— EPS IL 
DO 600 l — 1 » NNODE 
IPIRAUPTID.lt. TEST) 601,602 

601 CUNT iNuE 
RADPriiJ = 1.0 

GO TO oJO 

602 CONTINUE 

RADP TI I ) = 0.0 
600 CON T i NU E 
RETURN 

END uF ENCLUS 
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SUBROUTINE AROUNDIX , Y.RAD1 
C 

C *** DETERMINES A SPECIFIED CORNER ANGLE OF AN ELEMENT 
C 

DIMENSION XI 3) ,Y ( 31 ,RADt 31 

DONE = SQRTIIXI31-XI2) J**2+!YI31-Y(2) 1**2) 

OTwO = SORT! tXt3)-XI 11 I **2*(Y<3)-Y( 1) 1**2) 

OTRI = SwRT((Xl2)-X( m**2+lY(2)-Y( LI 1**21 
C OS I ■= lorw0**2*CTRI**2-00N£**21/I2.0*0TW0*0TRl 1 
C0S2 = 100Nt**2-»GTRI**2-UTWC**2l/ (2 . 0*Q0NE*0TRI 1 
COS 3 = tUUNE**2+CTWG**2-OTRI**2)/ I 2 . G*OONE*OTWQ I 
RAU( 1) = ACUSICCSI) 

RADI2! = ACOSICCS2) 

RADI3) = ACOSI COS31 
RETURN 

END OF ARUUND 


SUBROUTINE POINT St KIND, ISIZE.XPT, YPT) 

C 

*** TO PLOT SYMBOLS AT NODE POINTS 

COMMON/ KOUNT/ NNCOE , NEL MT , NNDEST, NELE ST 

COMMON/ CONTRL/ INFOR ,KGEOM , KOA TA, KS IGN , KPLOT , XORGN, YORGN, 

IP SC AC E, 1 SCALE, WMAG, I CN TR, XLHT ,NXLAB , XL ABI ID I , NYLAB , YL AB I 1 0 I 
DIMENSION XPTI II , YPT 1 1 1 
DO 1 DO I-ltNNi'DE 
X = XP T I 1 i / PSCALE 
Y - FLOAT! XSIGN1*YP II I I /PSCALE 
CALL PNTPLT IX, Y.KINO, ISIZE1 
I Du CONTINUE 
RETURN 

END OF POINTS 
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SUB ROUT INE WNOOE (NOMPT,XPT ,YP T , wP T , RADP T , NUMEL , NODE l , NQDE2 , 
lNODt3*NOOE4fCENTRJ 
C 

*** FOR PLOTTING NOCAL DATA 

COMMON/ KUUNT / NNDDE*NELMT f NNOEST* NELEST 

COMMON/ CONTRL/ I ISFCR .KG EQM , RCA TA , KS IGN.KPLOT , XORGN, YORGN, 
IPSCALt , ISCALt. wMAG, ICNTR , XLHT .NXL A6 , XLAB< 10) , NYL AB , YL AB ( 1 0) 
DIMENSION NUMPT-I 1 J , XPT ( 1 » , YP T C IJ j WPT I I J .RADPTUJ . NU MEL ( IJ, 
INuDEil lit NODE 2 ( 1 J .NOOE31 1) , N00E4I 1) .CENTRl IJ 
DIMENSION XX (3 J .YY13J , WW13J .RRAOl 3J 
DO 300 I- I * NELMT 
NO l = NODEUIJ 
NU2 = NC0E2II) 

NO 5 - NODE 3 I I J 
NU4 = ND0E4 ( I J 
X a I I ) ■= XPT ( NO 1 ) 

, YYI1J = YPTINOIJ 
WWl I) - WPT ( NO 1 J 
XX ( 2 J - XPTIND2J 
YY12J = YPTIND2J 
WW l 21 = WPTIND21 
XXI 3 J ^ XPT (NO 3) 

Y Y ( 3 J » YPT (ND3 ) 

WWl 3) - WPT (NO 3) 

CALL TRINlXX, YY.Wwl 
IFlND4.EW.0J GO TO 505 
XX l 2 J •= XPTIN04J 
YY ( 2 i - YPT ( ND4 J 
WW ( 2 J = WPTIND4J 
CALL TRINlXX, YY.WWJ 
505 CONTINUE 
300 CONTINUE 
RETURN 
END OF wNODt 
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SUBRUuT INE WELMT I NUMPT , XP T , YP T , WPT »RADP T , NUMEL , NODE 1 , NODE 2 , 
1 NUDe3,N0DE4,CENTR) 

*** FOR PLOTTING ELEMENT OATA 

COMMON / KUUNT / NNCDE,NELMT,NNOESP,NELEST 

COMMON/ CONTRL/ I NFOR ,KGEOM , KC A TA, KS ION, KPLOT , XORGN, YORGN. 

IP SCALE, I SC ALE, WM AG, I CNTR , XLHT , NXL AB , XL ADI 10) ,NYLAB, YLAB(IO) 
DIMENSION NUMPT ( i) ,XPT( 1) , YPT( 1) iWPTI 1) .RADPTI 1 ) , NUMEL ( 1) , 
INUOEit U,N0uE2( 1 ), NODES U ), NG0E4I U .CENTRM) 

DIMENSION NDSAVI4) , NEL ( 20 ) , XI 20) » Y( 20) ,R120) ,ANG( 20) 
DIMENSION XXI 3 ) , Y Y < 3 ) , rth ( 3 ) 

DIMENSION ISAV (4) , J S A V 1 4 ) 

PI = 3.1413426536 
DO 100 I— 1 , NNUUE 

*** TO DcTERMlNE ELEMENTS CONNECTED TO NODE l 

NKEEP = 0 
DO 300 J=1,NELMT 
NDSAV(i) = NODE 1 I J ) 

NDSAVI2) = N0DE2IJ) 

NQ5AVI3) - N0DE3IJ) 

NOS AV ( 4 ) = NUDE 4 ( J) 

DO 505 J J— 1 , 4 
IFINOsaVIjJ) . E u . I ) 510,305 
510 NKEEP = NnEEP+ 1 
NLLINKtEP) = J 
505 CONTINUE 
500 CONTINUE 

IF(NrEcP.EQ.O) GO TO 100 

*** CALCULATE INFORMATION FOR CENTROID OF EACH CONNECTED ELEMENT 

DO 530 K=1,NKEEF 
IDUM = NEL(K) 

NO 1 = NODEKIDUM) 

N02 = N0DE2 ( IDUM) 

NOj ■= NODES I IDUM) 

ND4 - NQDE4 I IDUM ) 

IFIND4.EQ.0J 560,561 

360 CONTINUE 

X ( K ) = (XPT IND1)+XPT(ND2)+XPT(ND3) ) /3.0 
Y { K ) = IYPT (ND1 ) ♦ YPT IND2I+YPT CNO J ) ) / 3 .0 
GO TO 565 

36 1 CONTINUE 

X(K) = (XPT (ND1 ) FXPT (ND2) +XPT IND3)+XP T (ND4) ) /4.0 
Y(tO = IYPT INDD+YPT (N02J+YPT 1 NU3.') -f YP T ( ND4) ) / 4.0 
5s5 CONTINUE 

XU 1ST = X1K)-XPTII) 

YU 1 S T = Y (K)-YPT ( 1 ) 

R ( N ) = SORT (XDIS7**2*-YDISJ**2) 

IFIXOISJ.EQ.O.O.AND. YD I S T . GE . 0 . 0) ANG(K) = Pl/2.0 
IF( XUIST.EQ.O. O.AND.YDl ST.LT.0.01 ANG(K) = 3.0*PI/2.0 
IFIXDISI.NE.O.O) ANGIK) = ATAN21 YDIST ,XDI ST ) 

IFIAnGIK) .lT.O.C) ANG(lO = ANG(K) +2.0*PI 
330 CONTINUE 
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*** to reorder infcrmat icn in oroer of increasing angle 

IFINKEEP.EG.il GO TO 60C 
DO 5 70 K= 1 , NKE E P 
KK = NKEEP-K 
DO 570 J= 1 , KK 

IFt ANGI JI-ANGI J+iJ ) 570,570,575 
575 TEMPI = ANG { J ) 

ITEMP2 - NEL(J) 

TEMP3 = XIJ) 

TEMP-+ = Y(J) 

TEMP5 = RIJI 
ANG I J ) = ANG (J + 1 ) 

N E L U ) = NELU + 1) 

XU) - XIJ+l) 

YU) = YU+1) 

K I J > = R I J + 1 ) 

ANGU + i) - TEMPI 
NEL(JFi) = ITEMP2 
XIJ+1) = TE MP3 
YU + i) = TEMP4 
R U 1 ) - TEMPS 
570 CONTINUE 
600 CON T I NUE 

*** TO CALCULATE WEIGHTED AVERAGE AT NODE POINT 


OENOM = 0.0 
DO o 1 0 K= 1 , NKE tP 
OENOM = OENOM+ l.Q/RIK) 

610 CuN TINGE 

wPTlIJ = 0.0 

00 615 K=1,NKEEP 
WEIGHT -= 1 1.0/RIK) ) /DENCM 
IUUM = NEUKI 

WPTII) = WPT( I )+WEIGHT*CENTR< IDUM) 
o 1 5 CONTINUE 

*** TO PLOT CONTOURS AROUND POINT 

IFINKEcP.EQ.l) GO TO 625 
XXII) ■= XPTIII 
YYI l) = YPT (I) 

WWU) = wPTll) 

DO 620 K^l.NKEEP 

1 l = K 
12 = K+l 

I F I K. EG .NKE EP i 12 = 1 
IFIRADPTID.EU.C.OI 660,650 
C *** TO CHECK CONNECTIVITY OF CENTROIDS ABOUT A BOUNDARY POINT 
650 CONTINUE 

IDUM = NELt I 1 J 
UDUM * NElI 12) 

ISAVt 1) = NODE 1 ( ICUM1 
I S AVI 2 ) = NODE 2 1 ICUM) 

ISAVI3) = NODE 3 I ICUM) 

ISAV14) = NODE 4 1 IDUM) 

JSAV(l) = NUDE 1 1 JOUM J 
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JSAVI2) ■= N0ufc2( JUUMJ 
J S AV ( 5 ) = NUOEJt JDUMi 
JSAV(A) = NUDE4U0UM) 

00 o3 L i I l- 1 * A 

I FI ISAV » I m .fcu.O.OK.i SAV( I lit .Eg. IJ GO TO 651 
L)U 652 JJ J=1 *4 

It( JiAV { J Jj ) .Ew.O.OR. JSAVI J JJ) .EU. i 1 GU TO 652 
lEUSAVUIi 1 .EQ.JSAVUJJ)) 660,652 
652 CONT INUE 
o51 CONTINUE 
GO TU 620 
66 0 CON T 1 NuE 

XX ( 2 1 = XdiJ 
YY12) = Yl I 1J 
IUUM = NEU ill 
WM{ 2) = CENTRl 10UMJ 
XXI3J = X ( I 2 J 
Y Y ( 3 J = YI12J 
1DUM = NEL112) 

3) CENTRC ICuMJ 
CALL TR INi XX, YYtlnWi 
620 CONTINUE 
o25 CONTINUE 
100 CONTINUE 
RETURN 

END OF *ELMT 



c. n r c, n r> o r» o o r r> 


APPENDIX C - Continued 


SUBROUTINE TR1N(XG»YG,WG) 

plots contour lues within 4 given triangle 

COMMON/ CuNTRL/ UFGR ,KGECM, KOATA.XS IGN.XPLOT , XORGN, YORGN, 
1 P,SCAll, ISCALE.WiMAG, I GNTR , XLHT ,NXLAB , XLAB ( 101 , NYL AB , YLAB ( 1 0 ) 
OIMENSICN XGI3 ) ,YG( 3) ,WGI3 J 
DIMENSION X I 3 I , Y(3),W<3J 
DIMENSION XPNTI 3) ,YPNT(3I 

*** TO SoAlE « DATA BEFORE PLOTTING 

DO 10 1=1,3 
XIII = XG ( I I / PSC ALE 
Yll) = FLOAT (KSIGNJ*YG( IJ/PSCALE 
will = WG I I ) *W M AG 
10 CONTINUE 

*=** TO FIND MIN AND MAX VALUES OF W 

WMIN = will 
W M A X = will 
DO BO 1=2,3 

IFUUJ .LT.WM1NJ WMIN = Will 
1F( Wi U .GT.WMAX) WM AX = will 
20 CONTINUE 

*** TO DETERMINE CONSTANT VALUES OF M, WITHIN TRIANGLE 

IDUM = WMIN/ICNTR 
ODOM = IDUM* 1 ICNTR 
1F( JDUM.GT.O) ISTRT = JCUM> ICNTR 
IFIJDUM.Lc.OJ ISTRT = JOOM 
IFIFLOATI jOUMJ .EC. WMIN) ISTRT = JDUM 
IDUM = WM AX / ICNTR 
JDUM = IDUM* ICNTR 
1F( JDUM.GE.OJ I STOP = JDUM 
IF( JDUM.LT.O) ISTOP = JOUM-ICNTR 
IF (FLUA.ri JOUM) .EG.WMAXJ ISTCP = JDUM 
IF( ISTOP. LT. IS TRTJ GO TO 10C0 
C 

c *** DETERMINES INTERSECTIONS OF CONSTANT M LINES AND TRIANGLE SIDES 
0 

ALLOW = 1 .OE-O S 

EPSILN = l.OE-lu 

ELUw = O.O-EPSILN 

EUP = 1.0+EPSILN 

VuAVE = 1.0E+20 

DO 100 1=1STRT , ISTOP, ICNTR 

WCONST = FLOAT! I) 

I ON T = 0 
DEN = Mil) — M(3i 
1FI OtN.Eu.O.OJ GO TO 120 
XL I = IwCGNST— W 13J J/UEN 
IFl XLl.GE.fcLOW.AND.XLl.LE.ECPl LI 0, 120 
110 I ON T = ICNT + 1 

XPNTIIONTJ = XL I*X( 1J + I 1.0-XLI )*X(31 
YPnTUOnTJ = XLI*Y( l) + i 1.0-XLI 1*Y.I3J 
120 CONTINUE 
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UEN = w(3)-Wt2) 

IFI UEN. tg.U.O) GC TU 140 
XLM = IwCGNST— Wt2) ) / UEN 
IFIXLM.GE.cLUW.ANC.XLM.LE.EUP) 130. 140 
130 I CN T = ICNTD 

XPNl(lCNT) - I 1.0-XLM)*X12)*XLM*XI3) 

YPNTdCNT) = l l.G-XLM) *Y(2 ) *XLM*Y.I3) 

140 continue 

DEN ^ « ( 2 ) - W ( 1 j 
IF! UEN. Ew.O.O) GG TO 16C 
XL J = { WCUNST-*( l))/OfcN 
1FIXLJ.GE.ELQW . AND. XL J.LE.EUP) 150.160 
150 ICiMT = iCNT+l 

XPNTliCNT) = I i.O-XLJ)*Xl 1}+XLJ*X(2) 

YPN T l IC N T ) = l 1.0-XLJ)*YI1) +XLJ*Yti2 ) 

160 CONTINUE 

1 F l ICNT.LE. 1) GO TO 100 
IFl ICNT.EQ.2) GO TO IbG 

iH AO 3 4 XPNT ( 1 ) -XPNT l 2) ) .L T . ALLOW- AND. 
1A6S(YPnT 11)-YPNT(2)) .LT .ALLOW ) IZQi 1B0 
170 C UN T l NUE 

XPNT! 2 ) = X PNT I 3 ) 

YPNTI2) = Y PNT ( 3 ) 

1 60 CONTINUE 

Call Calpltixpnku .ypnti l) ,3i 

CALL CALPLT I X P NT I 2 i .YPNT 12) » 2 ) 

CALL CALPlT UPNT (2) .YPNTI2) ,3) 

IFINXLAB.EQ.G) GC TO 300 

IFI AbSIXPNT ( 1J-XPNTI 2) ) .LT. ALLOW): GO TO 300 
GO Jiu K=1,NXLA8 
XLAoEL = XLAo(KJ/PSCALE 

lF(XPNIlL).L£.XLAEEL.AND.XPNTl2).Gt.XLABELi GO TO 350 
IFlXPNn2).LE.XLABEi.AAI0.XPM<lDGE.XLABEU GO TO 350 
GO TU 310 
350 CONTINUE 

XSIKT = XLA8EL 

YSTKT ^ YPNT I 1 )*{YPNT(2)-YPNT ( 1))*1XLABEL~XPNTI 1) J/ 

1 1 XPNT ( 2 J-XPNT ( IJ J 
YSTKT = YSTKT— XLFT/2.0 

I F l A3 S I VSAVE— YSTRT ) . LT . XLHT I GO TO 310 
VSAVE = YsTRT 

CALL NUMBEk IXSTRT, YSTKT, XLH.T.WCUNST ,0.0,- D 
310 CONTINUE 
30 U CONTINUE 

IF INYLAtJ. EU.U) GC TO 40C 

IFIABSI YPNT ( li-YPNT I2j ) .LT. ALLOW) GO TO 400 
DO 410 K = 1 * NYLAB 

YlAoEL = FLOAT I K S IGN ) # YLAB ( K ) /PSC AL E 

IFlYPNTtl).Lt.YLAUEL.ANC.YPNT(2)*GE.YLABEL) GC TO 450 
IFI YPNT(2).LE.YLABEL.AN0.YPNT(1).GE.YLABEL) GO TO 450 
GO TU 410 
430 CONTINUE 

YSTRT - YLAbEL 

XSTKT = XPNT! 1 ) MXPNT(2)-XPNT( ID *1 YLABEL-YPNTl 1) )/ 

1 ( YPNT (2 )— YPNT ID) 

XSIKT = XSTKT+XLFT/2.0 

IFIABSI VSAVE-XSTBT) .LT.XLHT) GO TO 410 
VSAVE = XSTRT 

CALL NUMBER C XSTKT .YSTRT, XLHT.WCONST, 90. 0,-1) 

410 CONTINUE 
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*00 

continue 

too 

CONT INUE 

1000 

CUNT I NUc 


UO 500 1=1,3 


wm = wm/wHAG 

500 

CON T I NUE 


RETURN 


END OF TRIN 


Continued 


108 



o r r n r* 


APPENDIX C - Continued 


SUBROUTINE BUUNUINUMPT, XPT , YPT , WPT, RaUP T , NUMEL , NODE 1 ,N00E2, 
i NuD E 3 , NCOE4 > CE N TP 1 
C 

C *** PLOTS bOUNDARY CP STRUCTURE AND LABELS CONTOUR LINES 
C 

CCMMUN/KOUNT/ NNCCE , NE LPT , NN BEST * NELEST 

COMMON/ CONT RL/ I NFOR ,KGECM , KDATA-, KS 1 GN , KPLOT , XOKGN, YORGN, 

IP SC ALE, i SCALE, WMAG, ICNT R , XLPT , NXL AB , X L AB ( 101 , NYLAB, YLAB( IOJ 
DIMENSION NUMP T C U ,XPT(1) ,YPT( U , wPTI 1J ,R AOPT (II , NUMEL I 1 1 , 
INODE 1 m , N00E2 ( l ) »N00E 3 ( 1 1 * NGUE4I 1J,CENTR(L) 

DIMENSION ND( 3 i , NOCONI 30, El 
DIMENSION XXI3) ,YY(3),WW(3) 

*** EXAMINE ALL- ELEMENT CONNECTIONS AND ADD 1.0 TO RADPTS WHEN TWO 
*** NODES ARE CONNECTED WHICH HAVE RADP T . NE .0 . 0. ALSO FLAG 
*** ELEMENTS CONTAINING SUCH CONNECTIONS BY MAKING NUMEL NEGATIVE. 

DO 4 0u 1 = 1 , NELMT 
NDI li ■= NODE 1 1 IJ 
NO! 2 1 = N0DE2I 1 J 
N 0 1 3 1 = N0DE3I I) 

ND(4l ■= NUDE4II) 

NDI 5 I = NODE 1(11 
DO 410 J= 1 , 4 
K = NDI J ) 

IFlK.Ew.DI GO TO 40C 
L = NOIU+11 

IF(L.Ew.O) L = N01J+21 

IF(RAUPTIK) .NE.O.O.AND.RADPTIL J.NE.O.O) 411,410 
411 CONTINUE 

RAUPTIN) = RADPTIK1+1.0 
RADP T l L 1 = RADPTIU + l.O 
IF! NUMELl 1 J .GT .0 J NUMEL < 1) = -NUMEL I I * 

*10 continue 
40U CONTINUE 
C ' 

C *** EXAMINE ALL ELEMENTS AND SAVE ALL CONNECTIONS WITH BOTH 
C *** ENDS HAVING RACPT.GE.4.0 IN NOCON. 

C 

NPDSS = 0 

DO 420 1=1, NELMT 

IFINUMELl I ) .GT .0 I GO TO 420 

Nul 11 = NODE 1 ( li 
NDI 21 = NODE 2 I 11 
NDI 31 = N0DE3 I 1 1 
NDI 4! = NUDE4I 11 
NDI 31 = NuUEll II 
DU 430 J= 1 » 4 
K = NUlJl 

IF(K.EQ.O) GU TO 42G 
L = NDIJ+I) 

IFIL.EU.01 L = NDIJ+21 

I F ( RADP I I K) .UE.4.G.AND.RA0PTIL1 . GE.4. 01 431,430 
4 3 1 CUN T I NU E 

iFlNPOSS.EO.SO) GC TO 430 
NPDSS = NPGSS+1 
NOCCNI NPOsS , 1 } = K 
NuCUNInPUSS ,2) = L 
430 CONTINUE 
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420 C UN i 1 N'J E 

*** CHECK POSSIBLE NC CONNECTS AND SAVE ANY MATCHED PAIRS 
NBAD - 0 

IFlNPOSS.LE.il 00 10 4/0 
I ENU = NPOSS-i 
DO 430 1 = 1 , I END 
JOU = I +1 , 

SEND = NPGSS 
DU 460 J=JGG»JENC 

IF! NOCUN ( 1,1). EQ.NUCGN! J, 1 ) . AND. NOCON ( I , 2 ) . EQ .NOCON ( J , 2 ) ) 
aGO TO *61 

l FI NOCON! I , 1) . EQ .NOCON I J , 2 ) .AND. NOCON (1,2) . EQ .NOCON I J , 1 ) ) 
A GO T 0 46 1 
Gu Tu 460 
461 CONTINUE 

NBAD = NBAD+1 
NOCONINBAU, 1) = NOCON! 1,1) 

NOCUN! NBAD»21 = NCCCNi 1,2) 

460 CONTINUE 
430 CONTINUE 
470 CONTINUE 

DO 100 1= 1 ,NELMT 
I F I NOME L I I ) . GT . 0 ) GO TO 100 
Null) = NUOElii) 

NO! 2) = NCDE2! i ) 

NO! 3) = NODE 3 ( I ) 

NO l 4) = N0UE41I) 

NO 13) = NUDE 1 ! I ) 

DO 30 J = 1 , 4 
K = NU(J) 

IFIN.Eu.O) GO TC 100 
L = NO! 0+1) 

1FIL.EQ.0) L = ND1J+2) 

IF! RAOPTIn) .GE. 3.C.AN0.RAUPT (D.GE.3.0) 111,30 
ill CONTINUE 
NO 1 = K 
N02 = L 

IF! J.G..1) 1DUM = J-l 
iFlJ.tw.l) I OUM = 3 
NAUX = NO! IOUM j 

*** TC PLOT BOUNDARY LINES 

IF! NBAD.EQ.Oi GO TO 40 
00 30 JJ = 1 ,NBAC 

1F1ND1. Ew.NUCUNl J J , 1 ) . AND .NC2 . EQ- NOCON 1 J J , 2 ) ) GO TO 50 
IFiNUi.EQ.NUCGN(JJ,2).AN0.N02.EQ.N0CQN(JJ,i)) GO TO 50 
30 CONTINUE 
40 CUNT INUE 

XI = XPT(N01)/PSCALE 
Y 1 -= FL0AT!kSIGN)*YPT!ND1)/PSCALE 
X2 = XPT IN02J/PSCALE 
Y 2 = FlUAT!kSIGN)*YPT1ND2)/PSCALE 
CALL CALPLT 1 XI ,Y 1 ,3) 

CALL CALPLT 1X2 ,Y2, 2) 

CALL CALPLT 1 X2 ,Y2 ,3 ) 

IF! INFUR. EQ. 2) 886,889 
8BB CONTINUE 
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*** TO PLOT CONTOUR LINES TC EDGE OF ELEMENT DATA PLOTS 

AXll) = XPTIND1) 

Y Y ( i i = YPTINUl) 

WXI il - WP T I imD I i 
XXI 21 = XPTIND2) 

Y Y l 2 1 = YPTIND21 
WW(21 - WPT l ND2 1 

N 1 - NO 1 1 J 
N2 - NUI2) 

Ni = NU131 
N*» ~ NUI'tl 
IFINA.EG.G1 910,920 
91 J lONTINUE 

XCENT = lXPTlNl)+XPTlN21+XPT{N31i/3.Q 
YCENT = ( YPTl Nl J+YPT (N2)*YPTIN3JJ / 3 . 0 

Go TO 950 

920 continue 

XCcNT = (XFTIN1 1 +XP T l N2 1*XP T l N3H-XPT I N4 1 j / A . 0 
YCENT = l yPTlNiJ*YPUN2l+YPT(N3H-YPT(N4) 1/4.0 
9 50 CON T i NU E 

XXI 3 J - XCENT 

Y Y l 3 1 = YCENT 
**131 = CENTklil 
CALL IkiNTXXfYY , W W 1 

ddv oun T i N<J c 
C 

C *#*■ TU LABEL CONTOUR LINES 
C 

wPTINOlJ ■= WPT INCH *WMAG 
WPTIND2) = WPT I NC2) *WMAG 
I F I WP T (NUli .EU.WPHND2) i GO TO 066 
wMI N = WP T l ND1 1 
«(MAX = WPTIN021 

IFInPTINUI) .GT. fcPTlND2) ) WMAX = WPTINOl) 

IFIwPT (NU2) .LT .WPTIND1) 1 wMIN = »*PT I ND2J 

I0UM-WMIN/1CNTR 

J GUM = I UUM* ICNTR 

IFI JUUM.GT.O j ISTRT = jDUM+ICNTR 

IFl JUUM.LE.O 1 ISTRT = JCUM 

IFIFLOAll JDUM) .EC.wMIN) ISTRT = JUUM 

IUUM = WMAX/ ICNTR , 

JUUM = IUUM* ICNTR 

I FI JUUM.GE.01 I STOP = JDUM 

IFI JUUM.LT.U1 ISTOP = JDUM— ICNTR 

IFiFLUATl JDUM) .EG. WMAX) ISTGP = JDUM 

IFI ISToP.LT. ISTRT) GO TO 666 

Y LINE - l.QE+20 

XlINE = I.GE+20 

IF IXPT I NU2) — XPT IND11 .NE.0.0) YLINE = YPTlNDll+lYPTlNU 21— YPT I NU 1 ) 1 
I* IXP TINAUXJ-XP UNO! 1 1/ I XPT I ND21-XPT INOi) ) 

I F l YPTl ND21-YPT INUl) .NE.0.01 XLINE = XP T ( ND 1 ) + (XP T ( ND21-XPT { ND l ) ) 
l* IYPT iNAUXl— YPT IND11 1/ l YPTI N0 2 l-YPT I NO 1 ) 1 
XU 1ST - XPT iNAUXl-XLINE 
YU I ST = YPT iNAUX)— YLINE 
THETA = 0.0 

IFIABSl YUIST/XUIST) .LT.1.0) THETA = 90.0 
XSAVE = I.0E+20 
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YSAVE = l.GE+20 

DO 600 JJ = ISTRT,ISTGP,ICNTR 

n CONST = FLOAT < J J J 

i F i wCGNST .LT. wMIN. UR.WCCNST.GI . wMAX J GO TO 600 
X = XPT IND1J+IXPT I ND2J— XP 1 1 ND 1 ) ) $ l WCONST— WP T I NU 1) } 

1/ ( MPT I ND2 J— WP T (NCI) i 
X = X/PSCALE 

Y = YPTINDDMYPII ND2)— YP T ( ND1 ) ).* I WCONST— WPT ( NQ1 ) J 
1/ ( WPT I ND2 J-WPT (AO 1 J i 

Y = FLOAT (ksign)*y/pscale 
DO 650 IG = 1,10 

I LUM = 10** IG 

IF l I ABS ( JJ) .LT. IDUM) 055, 65C 
o55 ND1G = IG 

IF(JJ.LT.01 NOIG = NUi GFl 
GD TO 660 
650 CONTINUE 
660 CONTINUE 

IFl THETA. EU. 0.01 GO TO 670 
IFITHETA.EG.90.0J GO TO 680 
670 CONTINUE 

IF l AbSl YSAVE-YJ .LE.XLHT J GO TO 600 
YNUM = Y-XLHT/2.C 

IF(XDIST.LE.Q.Q) XNUM = X+I 6 . 0/7. 0) *XLHT 
1 F ( XD IST.GT .0.0) XNUM = X-( t. 0/7.0) *XLHT*FLOAT (ND1G+1 1 
YSAVE = Y 

GO TO o 90 I 

680 CONTINUE / 

IF(AaSIxSAVE-X) .LE.XLHT J GO TO 600 
XNUM = X+-XLHT /2 .0 

1 F ( YD13T.LE.0.Q) YNUM — Y+ ( 6 . 0/7.0) *XLHT 

IF! YDIST.GT .0.0) YNUM - Y- ( b . 0/7. 0) *XLHT* FLOAT! ND IG* 11 

X S A V E = X / 

690 CONTINUE 

CALL NUMBER I XNUM , YNUM, XLHT , WCCNST »THETA,— 1) 

600 CGNT INUE 

6o6 CONTINUE - 

WPT(NOl) = WP T (NCI) /WMAG 
MPT i N 02 ) = WPT ( NC2 ) /WMAG 
50 continue 
100 CGnT INUE 

*** TO RESTORE RAUPT OF BOUNDARY POINTS TO 1.0 
*** AND NUMEL TO POSITIVE VALUE. 

Du 310 1=1, NNUDE 

lF(RAUPTll) .GT.O.CJ RADPTII) = 1.0 
3 10 CuNTINUE 

DO 315 1=1, NELMT 
NUMElH) = IABSINUMELI l) ) 

313 CONTINUE 
RETURN 

■ END OF BOUND- 
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SUBRUUT INE LAYOUT (NPLUT , NUMPT , XPT , Y PT , V.PT , RAuPT , NUMEL , NUD El , NQUE2, 
lNGDc3»NCDE4,CENTf<l 

*** TO PLUT LAYOUT CF ELEMENTS 

COMMON/ KGUNT/ NNGOE , NELMT , NN DE ST , NtL E ST 

COMMON/CUNTRL/ INFOK ,KGEGM , KCATA , KSIGN»KPLUT ,XGRGN, YORGN, 

IPSCALE* lSCACE»wMAG* ICNT R , XL h T , NXL A6 , XLAB I 10 1 ,NYLAB, YLA6 ( 10) 
DIMENSION NUMPTI 11 ,XPT ( 11 , YPTI 11 iWPIl 11 ,RADPT( 11 » NO MEL I 1) , 

INGuEl ill , NODE 2 ( 1) ,N00E3C 1) , N0DE41 11 ,CENTR( 1) 

C *** TO MAKE ALL GrIU POINT NUMBERS NEGATIVE 
UO 44 0 l-l»NNUUE 
NUMPTI 1 1 = -NUMPTI I 1 
440 CONTINUE 

UO 10U l- 1 » NtLMT 
NO 1 = NODE 1 I 1 1 
N02 ■= N0UE2II1 
N03 - NO LIE 3 1 1 1 
NU4 = N0UE4II) 

C *** TO MAKE ALL GRID POINT NUMBERS CONNECTED BY ELEMENTS POSITIVE 
iMUMPTINOll = I Afl SI NUMPTI Nil l 11 
NUMPTIN021 = IABSINUMPT (NG2 1 1 
NUMPTI NU31 = i ABSI NUMPTI N031 1 
NUMPTIN041 = IAfiSlNUMPT(ND411 
XI = XPTINO.il/PSCALE 
X2 - XP T (NU21/PSCALE 
X 3 = XP T ( N031/PSCALE 
Y 1 - FLUAT(KS1GN.)*YPTIN011/PSCAL£ 

Y 2 = FLOAT! KS l GM *YPT I ND2 1/ PSC ALE 

Y 3 = FLOAT IkS1GM*YPTINU31/PSCALE 
CALL CALPLT 1X1 ,Y1 ,31 

CALL CALPLT IX2,Y2,21 
CALL CALPLT(X3,Y3,21 
IFIN04.EQ.0J 200,300 
200 CONTINUE 

CALL CALPLTIX1 ,Y 1,21 
CALL CALPLT <X1 ,Y1,3) 

I F ( NPLUT . NE . 31 GC TO 25C 
C *** TO NUMBER TRIANGLES 

XCENT ■= IX1+X2+X 31/3.0 
YCENT = I Y1+Y2 + Y 31/3*0 
X - XCENT— (6. 0/7. Cl* XL FIT 

Y - YCENT-XLHT/2.C 
A - NUMELIIJ 

CALL NUMoERIX,Y,XLHT,A,0.0,-ll 
250 CONTINUE 
GU TU LOO 
300 CONTINUE 

X4 ■= XPTINU41/PSCALE 

Y 4 = FLOAT IKSIGM*YPT I NC4 1/ PSCALE 
CALL CALPLT 1X4, Y4, 21 

CALL CALPLT 1X1 , Y 1 ,2J 
CALL CALPLT(X1,Y1,31 
lFINPLOT.NE.il GC TO 350 
C *** TO NUMBER QUADRILATERALS 
XCENT = IXl+X2*X3+X4J/4.0 
YCENT = I Y1+Y2+Y3+Y41/4.Q 
X = XCENT-I 6.0/7 .Cl * XL FIT 
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y ~ rcewr-xLHr/a.c 

A = NUrtEL ( i } 

CALL NUMBER!*, Y,XLHJ,A,0.0,-1J 
3 30 CuiMl T NUb 
LUO CONTINUE 

L FLNFLOT.NE.i J GC TO 430 
C *** Tu NUMBER Gk ID POINTS 
UO 400 I=1»NNUDE 
IFINUMPTt li .Lt.Q) GQ TO 400 
X = XPf ( H/PSCALE+0.5*XLHT 
y = ELQAr(KSlGN)*yPHn/PSC/'LE«-0.5*XCHT 
A = NUMPT < l I 

CALL NUMBERLX, Y.XCHf, A, 0.0, -Li 
400 CONTLNUt 
430 CONTINUE 
RETURN 

tNU UF LAVOUr 
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APPENDIX D 


INPUT DATA FOR OBLIQUE ORTHOGRAPHIC PROJECTION OF EXAMPLE PROBLEM 


FLAT PLATfc UNOEk UN i FORM PRESSURE - GBLIGUE ORTHOGRAPHIC PROJECTIONS 
SUPTION 

NNOEST = j6i NUDISP = 1 , NVlISP = It NwDI SP = It 
KGEUM = 2 1 KD AT A = I, KPLGT = 3, IDCASE = 1 
$ 

QUADEL CwUAOi 

IU GILES, TEST 

APP DISPLACEMENT 

SOL i,0 

TIME oO 

CEND 

TITLE = FLAT PLATE UNDER UNIFORM PRESSURE 
OLOAD = ALL 
SPCFOftCES = ALL 

DISPLACEMENTIPRINT , PUNCH] = ALL 
STRESSIPRINT , PUNCH) = ALL 
LOAD = 100 
SPC = 100 
BEGIN BULK 


PUUAuI 

1 

1 

0.01 

1 

8.333E-91 

0.01 

o 

• 

o 

i-GARYl 

+GAR y 1 

u. GOB 

-0.005 







MAT 1 

1 

30. 

♦ 6 


0.3 0.3 




SPC1 

100 

123 

10 

20 

30 40 

50 

60 

+ S1 

+ SI 

70 

130 

1 90 

250 

310 320 

330 

340 

+S2 

*S2 

330 

3o 0 

120 

160 

240 300 




PL0AD2 

100 

100. 0 

fa 

THRU 

10 




PLUAD2 

100 

100. 0 

16 

THRU 

20 




PLUAD2 

100 

100.0 

20 

THRU 

30 




PLGAD2 

100 

100.0 

36 

THRU 

40 




PLQA02 

. 100 

100.0 

46 

THRU 

50 




GROSET 


0 



0 

6 



GRID 

10 


0.0 

0 . 0 

0. 0 




GRID 

20 


2.0 

0.0 

0.0 




GRID 

30 


4.0 

0.0 

0.0 




GR I D 

40 


6.0 

0.0 

0.0 




GR I D 

30 


8 . 0 

0.0 

0.0 




GRID 

GO 


10.0 

0.0 

0.0 




GRID 

70 


0.0 

2.0 

0.0 




GRID 

80 


2.0 

2.0 

0.0 




GRID 

DO 


4. 0 

2.0 

0.0 




GRID 

100 


6.0 

2.0 

0.0 




GkI u 

110 


8.0 

2.0 

0 . 0 




GRID 

120 


10 .0 

2.0 

0.0 




GRID 

130 


0.0 

4.0 

0 . 0 




GRID 

140 


2.0 

4.0 

0 . 0 




GR I D 

150 


4.0 

4.0 

0.0 




GRID 

160 


6. 0 

4.0 

0.0 




GRID 

170 


8.0 

4.0 

0.0 




GR ID 

180 


10.0 

4.0 

0 . 0 




GRID 

IDO 


0.0 

6.0 

0.0 




GRID 

200 


2.0 

6.0 

0.0 




GRIU 

210 


4.0 

6.0 

0.0 




GRID 

220 


o. 0 

6. C 

0 . 0 




GR I D 

230 


8.0 

6.0 

0.0 







115 



APPENDIX D - Continued 


GRIu 

240 


10.0 

6. 0 

0. 0 



GRID 

250 


0.0 

8.0 

0.0 



GRID 

260 


2.0 

6.0 

0. 0 



GRID 

270 


9.0 

8.0 

0.0 



GRID 

260 


6.0 

6.0 

0. 0 



GRID 

290 


8. 0 

8.0 

0.0 



GRID 

300 


10.0 

8.0 

0.0 



GRID 

310 


0. 0 

10.0 

0. C 



GRID 

32 0 


2.0 

10.0 

0.0 



GRID 

330 


9. 0 

10. 0 

0. 0 



GRID 

390 


6.0 

10.0 

0.0 



GRID 

350 


a. o 

10.0 

0.0 



GRID 

360 


10.0 

10.0 

0.0 



CQUADI 

. 6 

1 

10 

20 

80 70 



CQUADi 

7 

1 

20 

30 

90 80 



CuUADl 

8 

1 

30 

40 

100 90 



CQUADI 

9 

1 

90 

50 

110 100 



CQUADI 

10 

1 

50 

60 

120 110 



CQUADL 

16 

1 

70 

80 

140 130 



CQUADi 

17 

1 

80 

90 

150 140 



CQUADI 

18 

1 

90 

100 

160 150 



CQUADi 

19 . 

1 

LOO 

110 

170 160 



CQUADI 

20 

1 

UO 

120 

180 170 



CQUADI 

26 

l 

130 

140 

200 190 



CQUADI 

27 

l 

190 

150 

210 200 



CQUADI 

28 

1 

150 

160 

220 210 



CQUADI 

29 

1 

160 

170 

230 220 



CQUADI 

30 

l 

170 

180 

240 230 



CQUADi 

io 

1 

190 

200 

260 250 



CwUADl 

37 

1 

200 

210 

270 260 



CQUADI 

36 

1 

210 

220 

280 270 



CQUAui 

39 

1 

220 

230 

290 260 



CQUADI 

0 

1 

230 

240 

300 290 



CQUADi 

4o 

1 

250 

260 

32 0 310 



CQUADi 

97 

1 

260 

270 

330 320 



CQUADi 

98 

1 

270 

280 

340 330 



CQUADi 

99 

1 

260 

290 

350 340 



CQUADi 

50 

1 

290 

300 

360 350 



ENUOATA 








ENDGtUM 








Z 

- Dl SPLACEMENT S 






FURMAT 

1 (110, 

8X » 3E18 . 

6/) 






10 

u 

0. 


0. 

0. 

10 

-CUNT- 



-5. 7916 1GE+0 1 

5. 741610E+01 

0. 

11 


20 

G 

0. 


0. 

0. 

12 

-CUNT- 



2. 83183QE+03 

4.045395E+01 

0. 

13 


30 

G 

0. 


0. 

0. 

14 

-COi'-JT- 



4.962056E +03 

2 . 75 1 1 37E + 0 1 

0. 

15 


90 

G 

0. 


0. 

0. 

16 

-CUNI- 



4.462056E+03 

-2.75 11.37E + 0 1 

0. 

17 


50 

G 

0. 


0. 

0. 

18 

-CONF- 



2.631830E+03 

-4.045395E+01 

0. 

L9 


60 

G 

0. 


0. 

0. 

20 

-CUfiT- 



-5.7416 

10E+01 

-5. 741610E+01 

0. 

21 


70 

G 

0. 


0. 

0. 

22 

— CUD T— 



— 4. 0453 

55E+01 

—2. 83 1 830 E +03 

0. 

23 


80 

G 

0. 


0. 

5. 2842 50E + 03 

24 

-cun r- 



2. 1929 

13E+03 

-2. 192913E+03 

0. 

25 
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90 

G 

0. 

0. 

8. 351675E+03 

26 

-CON ^- 



3 . 472 7 62E + 0 3 

-7. 8002 94E+02 

0. 

27 


100 

G 

0. 

0. 

8. 351675E+03 

28 

-CONT- 



3.472762E+03 

7. 800294E+02 

0. 

29 


1 10 

G 

0. 

0. 

5. 2842 50E+03 

30 

-CONT- 



2. 19 29 1 3E +0 3 

2. 1929.1 3E + 03 

0. 

31 


120 

G 

0. 

0. 

0. 

32 

-CONT- 



-4.045 395E+01 

2. 8 31 8-30 E+ 03 

0. 

33 


130 

G 

0. 

0. 

0. 

34 

-CONT- 



-2. 751 137E+01 

-4.462056E+03 

0. 

35 


140 

G 

0. 

0. 

3.351675E+03 

36 

— CON T- 



7. 800294E+02 

-3.472762E+03 

0. 

37 


150 

G 

0. 

0. 

1. 324540E+04 

38 

-CGNT- 



1.25C892E+03 

'- 1 . 2 508 92E + 03 

0. 

39 


ioO 

G 

0. 

0. 

1. 32454GE+04 

40 

-CUN1- 



1.250892E+03 

1.2508 92 E+0 3 

0. 

41 


170 

G 

0. 

0. 

8. 351675E+03 

42 

-CONT- 



7. 800294E+02 

3.4727 62 E + 03 

0. 

43 


1 80 

G 

0. 

0. 

0. 

44 

-CONT- 



-2. 751 137E + 01 

4. 4620 56E + 03 

0. 

45 


190 

G 

0. 

0. 

0. 

46 

— CGiNT- 



2. 751 1 3 7E + 0 l 

— 4.462056E+03 

0. 

47 


200 

G 

0. 

0. 

8.351675E+03 

48 

-CGNT- 



-7. 800294E+02 

— 3.472 762E + 03 

0. 

49 


210 

G 

0. 

C. 

1. 324540E+04 

50 

-CGisi T— 



- 1. 250892E+C3 

-1.250892E+03 

0. 

51 


220 

G 

0. 

0. 

1 . 324540E + 04 

52 

-CCIM- 



- 1. 250892E+03 

1.250892E+03 

0. 

53 


230 

G 

0. 

0. 

8.351675E+03 

54 

-CONT- 



-7.800294E+02 

3.472762E+03 

0. 

55 


240 

G 

0. 

0. 

0. 

56 

-CONT- 



2.751 137E+01 

4 . 4620 56 E + 03 

0. 

5 7 


250 

G 

0. 

0. 

0. 

58 

-CONT- 



4. 045 395E+0 l 

-2.831830E+03 

0. 

59 


260 

G 

0. 

0. 

5. 2842 50E+03 

60 

-CONT- 



-2. 1929 1 3E + 03 

-2.1929.13E + 03 

0. 

61 


270 

G 

0. 

0. 

8.35 16 75E + 03 

62 

-CQNT- 



-3.472762E+03 

- 7 .8002 9<t E + 02 

0. 

63 


260 

G 

0. 

0. 

8. 351675E+03 

64 

-CON'I - 



— 3.472 7 6 2 E + 0 3 

7. 600294E+02 

0. 

65 


290 

G 

0. 

0. 

5. 2842 50E+03 

66 

-CUNT- 



-2. 19291 3E+ 0 3 

2. 1929.13E+03 

0. 

67 


300 

G 

0. 

0. 

0. 

68 

-CONT- 



4. 045 39 5E + 0 1 

2 . 831 8 30E + 03 

0. 

69 


310 

G 

0. 

0. 

0. 

70 

-CONT- 



5. 741610E+0 1 

5. 741610E+01 

0. 

71 


320 

G 

0. 

0. 

0. 

72 

-CONT- 



-2. 831830E+03 

4.045395E+01 

0. 

73 


330 

G 

0. 

0. 

0. 

74 

-CONT- 



-4.462056E+03 

2.751137E+01 

0. 

75 


340 

G 

0. 

0. 

0. 

76 

-CONT- 



-4. 4620 56 E + 0 3 

-2.751L37E+01 

0. 

77 


350 

G 

0. 

0. 

0. 

78 

-CON T — 



-2.83183CE+03 

-4 . 045 3 95E + 0 1 

0. 

79 


360 

/■» 

O 

0. 

0. 

0. 

80 

-CCINT- 



5. 741610E+01 

-5.741610E+01 

0. 

81 


(BLANK CARD) 
(BLANK CARD) 
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$P LOT 
KHOR2 = 

I SC AL t = 

KDlSP = 

$ 

SP I CT 
KD ISP = 

$ 


L, KVERT = 
2, XUHGN = 
i , I DMAG = 


3, PHI = 20.0, THETA = 0.0, 
1.0, YORGN = 1.0, PSCALE 
2, DMAGS = 2.0, KODE = It 


PS l 

2 . 0 , 


3, KGOb - 0 , 


10. Ot 
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APPENDIX E 


INPUT DATA FOR CONTOUR PLOT OF EXAMPLE PROBLEM 
FLAT Plate under uniform PRESSURE - CONTOUR PLOT 


$UPT I UN 

nnuest 
kplot = 
$ 

FORMAT 
GR ID 

- 36* NELtST = 2 3* 
it IDCASE = I* 

(6X, 13, i:>X,2Fa.l 1 
ID 0.0 

0.0 

0.0 


GRID 

20 

2.0 

0.0 

0.0 


GRID 

30 

5.0 

0.0 

0.0 


GRID 

50 

6.0 

c.o 

0.0 


GRID 

50 

6. 0 

c.o 

U.C 


GRID 

60 

10.0 

c.o 

o.c 


Grid 

70 

0. 0 

2.0 

o.c 


GRID 

60 

2.0 

2.0 

o.c 


GRID 

90 

5.0 

2.0 

o.c 


GRID 

100 

6.0 

2.0 

0.0 


GRID 

HO 

6.0 

2.0 

0.0 


GRID 

120 

10.0 

2.0 

0.0 


GRID 

130 

0.0 

5.0 

0 . 0 


GRID 

1 90 

2. 0 

5.0 

0.0 


GRID 

150 

5.0 

5.0 

o.c 


GRID 

loO 

6.0 

5.0 

o.c 


GRID 

170 

6.0 

5.0 

o.c 


GRID 

160 

10. 0 

5.0 

o.c 


GRID 

150 

0.0 

6.0 

o.c 


GRID 

200 

2. 0 

6.0 

o.c 


GRIO 

210 

5.0 

6.0 

0.0 


GRID 

22 0 

6.0 

6.0 

0.0 


GRID 

230 

6.0 

6.0 

0.0 


GR I D 

250 

10.0 

6.0 

o.c 


GRID 

230 

0.0 

6.0 

0.0 


GRID 

260 

2.0 

6.0 

0 . c 


GRID 

27 0 

5.0 

6.0 

0.0 


GRID 

260 

6.0 

6.0 

0.0 


GRID 

290 

6.0 

6.0 

0.0 


GRIu 

300 

C 

• 

o 

8.0 

0 . c 


GRID 

310 

0.0 

10.0 

o.c 


GRID 

32 0 

2.0 

10.0 

o.c 


GRID 

330 

5.0 

10.0 

o.c 


GRID 

350 

6.0 

10.0 

c. c 


GRID 

33 0 

6. 0 

10. c 

0.0 


GRID 

560 

10. 0 

10.0 

o.c 


ENDGKID 

FORMAT 

CQUADI 

l 6X, 12 

6 

, i 5X * 9 1 1 3 * 5X ) ) 
1 10 

20 

60 

70 

CQUADI 

7 

1 20 

30 

9C 

60 

CQUADI 

6 

1 30 

50 

100 

90 

CUUADI 

9 

1 50 

50 

1 10 

100 

CQUADI 

10 

1 50 

60 

120 

no 

CQUADI 

lo 

1 70 

60 

150 

130 

CQUADI 

17 

1 60 

90 

1 5 0 

150 

CUUADI 

16 

1 90 

100 

L6C 

150 
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APPENDIX E - Concluded 


CQUADI 

19 

1 

100 

110 

CQUA01 

20 

1 

110 

120 

CQUAOi 

26 

1 

130 

140 

CQUAUL 

27 

1 

140 

150 

CqUAOI 

28 

1 

150 

160 

CQUADI 

29 

1 

160 

170 

CQUADI 

30 

1 

170 

180 

CQUAOi 

36 

1 

190 

200 

CQUAOI 

37 

1 

200 

210 

CQUADI 

36 

1 

210 

220 

cquadi 

39 

1 

220 

230 

CQUAOi 

40 

I 

230 

240 

CQUADI 

46 

1 

250 

260 

CQUAOI 

47 

1 

260 

270 

CQUADI 

48 

1 

270 

280 

CQUADI 

49 

1 

280 

290 

CQUADI 

50 

I 

290 

300 

ENDGEUH 





l 

- DI iPLACEMc 

NTS 


FORMAT 

4 ( 4 i 15, 

El !>• 

611 


10 

0. 


20 0 . 


50 

0. 


8 0 0 . 


90 

8. 35 16 75t «• 03 

100 6. 

351675E+Q3 

130 

0. 


140 8 . 

35167 5E * C3 

170 

8. 35 1 6 75t *■ C3 

LuQ 0 . 


210 

1.324540E404 

22 0 1 . 

324540E+C4 

250 

0 . 


260 5. 

284250E+03 

290 

5« 284250E+03 

300 0. 


030 

0 . 


340 0. 



(BLANK CARD) 
spier 

pscAie - 2.0, 

NXi.AH ^ I, AlABlil = 5.0, NYCAb = i, 
$ 


170 

160 


180 

170 


200 

190 


210 

200 


22C 

210 


230 

220 


240 

230 


260 

250 


270 

260 


280 

270 


290 

280 


300 

290 


320 

310 


330 

320 


340 

330 


35C 

34 0 


360 

350 



30 

0 . 

40 

0 . 

7 C 

0 . 

80 

5.2842 50E+03 

110 

5. 264250E403 

120 

0 . 

150 

1.324540E+04 

160 

1. 32 45 4 OE * 04 

190 

0 . 

200 

8.35 1675E+03 

230 

6. 35 1675E* 03 

240 

0 . 

270 

8.351675E+03 

280 

8.351675E+03 

310 

0 . 

320 

0 . 

350 

0 . 

360 

0 . 


YLABll) = 5.0 
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APPENDIX F 


LISTING OF CONTOUR PROGRAM OUTPUT FOR EXAMPLE PROBLEM 

FLAT PLATt UNDER UNIFORM PRtSSURE - CONTOUR PLOT 


$OPT I ON 



nnuest 

- 

30 * 

NELEST 

= 

15, 

KGEOM 

- 

i t 

KDATA 

= 

1. 

NVALUS 

= 

0, 

IREStQ 

- 

i t 

KPLUT 

= 

3 , 

INF OR 

= 

i . 

XSPACE 

.•= 

0 . IE+02 » 

RS1GN 

- 

1# 

iOCASE 

- 

l t 

iEND 




BLANK. COMMON STGRAGE ill REQUIRES AT LEAST 330 LOCATIONS FOR THIS CASE 


GRID POINT INFORMATION 


RESEQUENCED 

USER INPUT 


GRID POINT 

oRIU POINT 

INPUT 

number 

NUMBER 

X 


BOUNDARY 

INPUT POINT 

Y INDICATOR 


1 

l 

3 

A 

5 

6 


iu 

0. 

0. 

20 

2 .OOCOE+OO 

0. 

30 

A.OOOCE+OO 

0. 

AO 

O.0000E+00 

0. 

50 

8.00C0E+CG 

0. 

60 

l.OOOOE+CI 

0. 


i.OOOOE+OO 
1. QOOOE+QO 
i.OOOOE+OO 
l.OOOOE+OO 
I. OOOOE+OO 
l.OOOOEFOQ 
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APPENDIX F - Continued 


7 

70 

0. 

2. OOOOE + OO 

1 . OOQOEtOO 

a 

80 

2.00CGE+C0 

2. OOOOE + OO 

0. 

5 

90 

4.0000E+CG 

2.0QQGE+00 

0. 

lu 

100 

6.GUCGE+C0 

2 . QQQQE+OQ 

0. 

it 

110 

8. OOOOE+OO 

2.0G00E+Q0 

0. 

12 

120 

i.OOCOE+Ci 

2. QQQOE+OQ 

1 . OOOOE+OO 

13 

130 

0. 

4.QG00E+Q0 

l . OOOOE+OO 

14 

140 

2 .OOQCE+OO 

4. OOOOE+OO 

0. 

15 

150 

4.0UCGE+00 

4.0Q00E+00 

0. 

16 

160 

6.00CCE+00 

4. OOOOE+OO 

0. 

1 7 

1 76 

8.00CCE+CC 

4. OOOOE+OO 

0. 

id 

180 

i.GOGOE+Cl 

4 . OQOOE+OQ 

1. OOOOE+OO 

19 

190 

0. 

6. OOOOE + OO 

1. OOOOE+OO 

20 

200 

2 .OOGOE + 00 

6. OOOOE+OO 

0. 

21 

210 

4.00C0E+C0 

6. OOOOE+OO 

0. 

Zz 

220 

6 .OOGOE + 00 

6 . OQOQE+OQ 

0. 

23 

2 30 

3.000QE+CO 

6.0QG0E+00 

0. 

24 

240 

l.OOOOE+Ol 

6 . OOOOE <-00 

1. OOOOE+OO 

2 5 

280 

0. 

d.QOOOE+OO 

1. OOOOE+OO 

26 

2o0 

2 . OOCOE + CO 

8.0QU0E+00 

0. 

27 

270 

4.0000E+GG 

8. OOOOE *-00 

0. 

2d 

280 

6.0QCOE+GO 

8 . QOOOE+OO 

0. 

29 

290 

d.OUCOE+OO 

8. OOOOE+OO 

0. 

30 

300 

l.GQCGt+01 

8. OOOOE+OO 

1 . OOOOE+OO 

31 

310 

0- 

1 .0000E>01 

1. OOOOE+OO 

32 

320 

2. OOCOE + CO 

1 • OOOOE +0 1 

1. OOOOE+OO 

33 

330 

4.0000E+00 

1. OOOOE +01 

1. OOOOE+OO 

34 

340 

6.00CCE+CG 

l.OOOOE+Ol 

1. OOOOE+OO 

35 

350 

d . OGGOE + OO 

1 .0G00E+01 

1. OOOOE + OO 

36 

360 

1 . QOCGE + O 1 

I.OOOOE+Ol 

1. OOOOE+OO 

ELEMENT 

INEGRM AT10N - WITH 

RESEQUENCEO GR 10 

PO IN TS 



RESEOUENCEU 

USER INPUT 


GRID POINTS 


tL E(1 ENT 

ELEMENT 




NUMdER 

NUMdER 

1 

2 3 

4 


1 

6 

1 

2 8 

7 

2 

7 

2 

3 9 

8 

3 

8 

3 

4 10 

9 

4 

9 

4 

5 11 

10 

5 

10 

5 

6 12 

11 

6 

16 

7 

8 14 

13 

7 

17 

8 

9 15 

14 

8 

18 

9 

10 16 

15 

9 

19 

10 

11 17 

16 

10 

20 

11 

12 18 

17 

11 

2o 

13 

14 20 

19 

12 

27 

14 

15 21 

20 

13 

28 

15 

16 22 

21 

14 

29 

16 

17 23 

22 

15 

30 

17 

18 24 

23 

1 6 

36 

19 

20 26 

25 

1 7 

37 

20 

21 27 

26 
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APPENDIX F - Continued 


18 

36 

21 

22 

28 

27 

19 

39 

22 

23 

29 

26 

20 

40 

23 

24 

30 

29 

21 

46 

25 

26 

32 

31 

22 

47 

26 

27 

33 

32 

23 

46 

27 

28 

34 

33 

24 

49 

26 

29 

35 

34 

25 

50 

29 

30 

36 

35 


2 - UI SPLACEMENTS 


Li AT A TO BE PLGTltO 


RESEQU ENDED GR 10 POINT NUMBERS AND DATA VALUES 


1 

0* 

2 

0. 

3 

0. 

o 

Q • 

7 

0. 

6 

5 • 28425E+ 03 

11 

5 .2 6425C+03 

12 

c. 

13 

0. 

1 6 

1 . 324 54E+ 04 

17 

6.35167E+03 

18 

0. 

21 

1.32454E*04 

22 

1 « 32454E+Q4 

23 

8.35167 E*-03 

26 

5. 26425E+03 

27 

6.35 167E + Q3 

28 

8.35167E+03 

31 

0. 

32 

C. 

33 

0. 

36 

0. 

1 

f 4 

0 . 

5 

0 . 

These items appear to the 

9 

14 

6. 35167E+03 
8.35167E+03 

10 

15 

8.35167E+03 

1.32454E+04 

right of the above items in 2 

19 

0. 

20 

8.35167E+03 

the actual printout 

24 

29 

0. 

5.28425E+03 

25 

30 

0 . 

0. 



34 

V. 

0. 

35 

0. 


$P I C T 



NPLOT 

- 

4, 

XGRGN 

= 

o 

• 

o 

YGRGN 

- 

0 • 0 9 

PSCALE 

- 

0 . 2 E *• 0 1 t 

ISCALc 

= 

3 9 

RMAGS 


0. 1E+03, 

ICNTRS 


10, 
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APPENDIX F - Concluded 


XLH r 

0.15E+00, 






NXLAB 

It 






XLAB 

0.5t«-01 , 

0 • 0 » 0 • 0 , 0 « 0 f 

o 

• 

o 

.0.0. 

o 

• 

o 

0.0, 0.0, 0.0, 

NYLAB 

It 






YLAB 

0.5E+01 t 

J.Of O.C-9 0-0, 

0.0, 

o 

• 

o 

— 

o 

• 

o 

0.0, 0.0, 0.0, 

KOUt 

Ot 






$ENU 







MAXIMUM 

VALUE UP 

DATA = 1.3245400CE+04 

MINIMUM 

VALUE 

OF DATA = 0. 

wMAO = 

7. 549790B7E-G3 

ICNTR = 

10 
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(b) Element numbers. 

Figure 4.- Exploded plots of elements at a selected fuselage frame-wing spar location of figure 1. 
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(b) View from side. (c) View from rear. 

Figure 5.- Contour plot and oblique orthographic projections of vibration mode shpaes of aircraft in figure 1. 
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Figure 7.- Diagram showing subareas of a triangle required 
to locate a generic point p by using area coordinates. 



REMAINDER OF INPUT DECI 
REPEATS PARTS OF BASIC 
DECK DEPENDING ON VALU 


° £ 
ll II O 

i— i 

UJ I— 
£- is) 
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§ 2 g 

i— i «=- 
< 



—I Q_ 
< O 


</> ^ ^ 
5 5 q 

■Z ^ _ j ^ 

UJ Z 3 < 

Q- O O 

Q M i, Ll 


O LU I 

ZOuj 

U. ^ ^ 


Q O 
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Figure 8.- A proper sequence for input cards. 







(a) Deformed shape. 



Figure 9.- Oblique orthographic projections of a simply supported, 
square, flat plate under uniform pressure. 




k < x k.yk. z k> 



l(x i.yi.Zi) C j (xj.yj, zj) 

A,- = area of Ajkp 

Aj = area of A kip 

A k = area of Aijp 

A = area of Aijk 

Figure 11.- Diagram showing subareas of a triangle required to 
locate its incenter p by using area coordinates. 


136 




Call GE0M1. GEOM2, or GEOM9 
to read geometry input data deck 


I^ZEROD to initialize 
displacements to zero 


Call DATA1, DATA2, or DATA9 to 
read deck of data to be plotted 


Read NAMELIST P1CT> 


Call DSCALE to determine 
scale factor for displacements 


Set pointers to allocate 
storage in blank COMMON 



PNTOUT prints data \ 
to be plotted / 


Call CCRT1 to change 
values in NAMELIST PICT 



Call BOUND to determine maximum 
dimensions of model to be plotted 


Call ROTAT to calculate 
coefficients of rotation matrix 


Call PLOT which uses XYSCAL, XPLOD, \ 
and GARROW to generate plots / 


Call NDLET to annotate 
grid points on plot 
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( Start program MAIN ) 


/ Call INITAL to set default values \ 
\ and read NAME LIST OPTION / 


Call GEOM1, GEOM2, or GEOM9 \ 
to read geometry input data deck / 


Call XYSCAL to determine bounds' 
of datum surface and renumber 
^element grid points 


/ Call ENCLOS which uses AROUND \ 
\ to determine boundary points • / 




Call ZEROW to initialize all 
data to be plotted to zero 


> 


Call DATA1, DATA2, or DATA9 . 
to read deck of data to be plotted 


> 


( Read NAMELIST PICT) 


Call SCALEW which uses INTRVL' 
to scale data to be plotted and 
^determine contour interval 


Call POINTS, WNODE, or WELAAT, 
and BOUND which use TRIN to 
^generate contour plots 


KODE = 0 


< 


DOCMNT contains 
user documentation 


) 


- ( Read title card for run ) 


< Set pointers to allocate \ 
storage in blank COMMON / 


— ( PNTOUT prints geometry data ) 

< Read case identification \ 
card if IDCASE ^ 0 / 


_/ PNTOUT prints data 


to be plotted 


) 


Y Call CCRT1 to change \ 
A values in NAMELIST PICT/ 



< Call CCRT2 to change KODE ) 


( Call P.STOP o terminat7^i> — — Xlfpiot vector ^ E0F ) 


( End program MAIN ) 

Figure 13.- Flow chart for program which generates contour plots. 
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